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Abstract 


The  flap-lag  equations  of  motion  of  an  isolated  rotor 
blade  and  those  for  a  rigid  helicopter  containing  four  blades 
free  to  flap  and  lag  are  derived.  Control  techniques  are 
developed  which  stabilize  both  systems  for  a  variety  of 
flight  conditions. 

Floquet  theory  is  used  to  investigate  the  stability  of  a 
rotor  blade’s  flap- lag  motion.  A  modal  control  technique, 
based  on  Floquet  theory,  is  used  to  eliminate  the  blade's 
instabilities  using  existing  collective  and  cyclic  pitch 
control  mechanisms.  The  technique  shifts  the  unstable  roots 
to  desired  locations  while  leaving  the  other  roots  unaltered. 
The  control,  developed  for  a  single  design  point,  is  shown  to 
significantly  reduce  or  eliminate  regions  of  flap-lag 
instabilities  for  a  variety  of  off-design  conditions.  Both 
scalar  and  vector  control  are  successfully  used  to  stabilize 
the  blade’s  motion. 

Coupling  the  flap-lag  equations  of  motion  of  four  rotor 
blades  to  a  rigid  airframe  alters  the  flap,  lag,  and  airframe 
roots.  The  airframe  roots  are  stabilized  using  a  combination 
of  the  body’s  pitch  attitude  and  pitch  rate  feedback  to  the 
main  rotor's  longitudinal  cyclic  pitch.  The  modal  control 
technique  is  used  to  eliminate  multiple  blade  instabilities 
by  first  controlling  a  pair  of  unstable  roots  at  a  specific 
design  point.  The  resulting  closed  loop  system  is  a  new 
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linear  system  with  periodic  coefficients.  Another  modal 
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TIME  PEBIODIC  CONTROL  OF  A  MULTI -BLADE  HELICOPTER 

I .  Introduction 

A  helicopter  is  a  mechanically  complex  aircraft  whose 
stability  and  control  characteristics  are  often  marginal 
unless  a  reliable  automatic  control  system  is  used. 
Effective  modeling  of  dynamic  effects  is  crucial  to  the  safe 
design  of  a  successful  helicopter.  It  is,  therefore, 
essential  to  consider  the  basic  dynamic  behavior  of 
helicopters  so  that  potential  instabilities  can  be  simulated 
and  eliminated  in  the  design  process.  However,  the  dynamics 
of  a  helicopter  in  hover  and  in  forward  flight  involve  the 
coupled  motions  of  the  rotor  blades  and  the  fuselage.  It  is 
convenient  to  analyze  and  understand  the  helicopter  rotor 
blade  and  fuselage  dynamics  separately  before  their  coupling 
is  examined. 

Rotor  Blade  Dynamics 


Many  researchers  have  examined  the  stability  of  an 
isolated  rotor  blade  in  hover  and  in  forward  flight.  The 
helicopter  blade  can  be  modeled  as  a  rigid  body  rotation 
about  the  blade  root,  which  is  attached  to  the  rotor  hub 
(1;2).  The  blade  has  three  degrees  of  freedom:  flap;  lag; 
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and  pitch  (or  feather)  rotations  about  hinges  at  the  blade 
root  (see  Figure  1)  .  Vote  that,  in  the  figure,  Ci  is  the 
rotor’s  rotational  speed  while  ft  is  the  blade  flap  angle  and 
is  positive  for  upward  notion  of  the  blade.  The  blade  lag 
angle,  C.  is  positive  when  it  is  opposite  the  direction  of 
the  rotor  blade’s  rotation,  and  6,  the  blade  pitch  angle,  is 
the  feathering  motion  produced  by  rotating  the  blade  about  a 
hinge  at  the  root  and  is  positive  for  nose-up  blade  rotation. 
The  blades  on  a  helicopter  are  maintained  in  uniform 
rotational  motion,  where  the  rotation  direction  is  assumed  to 
be  counterclockwise  as  seen  from  above. 

Two  fundamental  types  of  mathemat ical  models  are  used  to 
describe  the  flap  and  lag  motion  of  a  blade.  The  first  model 
consists  of  a  hinged,  spring-restrained,  rigid  blade.  Hinges 
at  the  blade  root  allow  free  motion  of  the  blade  both 
perpendicular  to  and  in  the  plane  of  rotation.  A  schematic 


of  the  hinge  arrangement  is  shown  in  Figure  2. 

A  flapping  hinge  alleviates  the  root  stresses  and  hub 
moments  by  permitting  blade  motion  out  of  the  disk  plane. 
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Figure  2.  Rotor  Hinge  Arrangement  (2:7) 


However,  this  flapping  motion  introduces  aerodynamic  and 
inertial  forces  in  the  disk  plane,  so  a  lag  hinge  is 
introduced  to  reduce  chordwise  root  loads  by  allowing 
in-plane  motion.  A  pitch,  or  feathering,  hinge  is  also 
required  to  control  the  rotor  by  changing  the  blade's  angle 
of  attack.  On  a  hinged  blade,  the  pitch  hinge  is  usually 
outboard  of  the  flap  and  lag  hinges  (Figure  2).  Motion  about 
the  flap  and  lag  hinges  is  restrained  by  centrifugal  forces, 
while  the  control  system  restrains  motion  about  the  pitch 
hinge  (2)  . 

The  second  model  describing  blade  motion,  applied  to 
hingeless  rotors,  treats  the  blade  as  an  elastic  beam.  The 
blade  is  attached  to  the  rotor  hub  without  flap  or  lag 
hinges,  though  there  are  often  hinges  for  the  feathering 
motion.  The  blade  is  attached  to  the  rotor  hub  by  a 


cantilever  root  restraint,  and  blade  motion  occurs  through 
bending  at  the  root.  The  structural  stiffness  is  small 
compared  to  the  centrifugal  stiffening  of  the  blade;  hence, 
the  fundamental  structural  mode  shape  is  not  too  different 
from  the  flapping  motion  of  a  hinged,  rigid  blade. 
Therefore,  the  remainder  of  this  study  deals  with  the  motion 
of  a  hinged,  rigid  blade. 

The  basic  derivation  of  a  helicopter  blade's  equations 
of  motion  has  been  well  documented  (1-5).  The  simplest  model 
considers  an  isolated  rotor  blade’s  flap  dynamics  in  hover 
and  in  forward  flight.  Flapping  instabilities  at  moderate  to 
high  advance  ratios  have  been  extensively  examined  using  this 
model  (6-10).  Sissingh  and  Kuczyncki  (11),  as  well  as 
Hohenemser  and  Prelewicz  (12),  extended  their  investigations 
of  flapping  instabilities  to  include  the  effect  of  the  blade 
torsional  degree  of  freedom  on  the  rotor’s  stability 
boundaries  and  response  characteristics . 

Various  procedures  have  been  incorporated  to  analyze  an 
isolated  rotor  blade’s  equations  of  motion  and  obtain  the 
rotor's  flapping  stability  limits.  Some  of  these  methods 
have  included:  generating  time  history  responses  by  analog 
or  digital  Integration  of  the  equations  of  motion  (8;9;11); 
analytically  solving  the  system's  equations  of  motion  by  the 
rectangular  ripple  method  (10);  and  modeling  the  helicopter 
blade  as  a  linear  system  with  periodically  varying 
coefficients  and  using  Floquet  theory  to  determine  the 
system's  characteristic  exponents  (7).  Hohenemser 


and 


Prelewicz  (12)  discovered  that  a  simple  system  identification 
method  applying  a  linear  sequential  estimator  could  be  used 
to  solve  the  transient  responses  in  the  flap-bending  and 
torsion  of  a  rotor  blade.  Even  with  these  varying  methods  of 
solution,  all  of  the  studies  indicated  flap  instabilities  at 
high  advance  ratios  and  therefore  considered  the  effects  of 
reversed  flow. 

To  examine  a  rotor  blade’s  stability  at  more  realistic 
advance  ratios,  investigations  have  also  been  extended  to 
include  analyzing  the  blade’s  coupled  flap-lag  motion 
(13-22).  Friedmann  and  Silverthorn  (13)  examined  the 
flap- lag  motion  of  a  cantilevered  rotor  blade  at  arbitrary 
advance  ratios.  They  concluded  that  instabilities  predicted 
by  flap  only  models  could  be  inaccurate  due  to  the  neglected 
lag  degree  of  freedom  and  that  coupled  flap-lag  instabilities 
occurred  at  lower  values  of  the  advance  ratios  than  when  only 
the  flap  degree  of  freedom  was  considered.  Peters  (15) 
investigated  the  flap-lag  stability  of  hinged  rigid  blades  in 
forward  flight  and  concluded  that  flap-lag  stability  was  very 
sensitive  to  the  moment  trim  and  propulsive  trim  conditions 
placed  on  the  rotor.  Further,  for  advance  ratios  between 
0.10  and  0.50,  excluding  periodic  coefficients  and  using 
constant  coefficient  equations  of  motion  resulted  in 
erroneous  flap-lag  stability  predictions.  Kaza  and  Kvaternik 
(16)  examined  the  stability  of  a  hinged  blade  in  hover  and  in 
forward  flight  for  both  flap-lag  and  lag-flap  hinge 
sequences.  These  authors  discovered  that  the  form  of  the 


coupling  terms  in  the  helicopter  blade’s  equations  of  motion 
depended  on  the  hinge  sequence  used  which,  consequently. 


influenced  the  blade’s  stability.  Finally,  studies  (2; 15) 
reported  that  the  effect  of  the  reversed  flow  region  was 
negligible  for  advance  ratios  less  than  0.50. 

An  improved  model  of  a  helicopter  blade  is  the 
introduction  of  a  third  degree  of  freedom,  torsional  motion 
(which  is  also  known  as  feathering  rotation).  Again,  studies 
have  examined  the  stability  of  coupled  f lap-lag- torsion 
motion  on  hinged  rotors  (23-26) .  Panda  and  Chopra  (26) 
investigated  the  effects  of  flap  bending,  lag  bending,  and 
torsion  on  the  stability  of  an  isolated  blade  in  forward 
flight.  They  found  that  raising  the  blade’s  torsional 
stiffness  increased  lag  mode  damping  at  lower  forward  speeds 
(advance  ratios  less  than  0.30)  and  decreased  lag  mode 
damping  at  higher  forward  speeds.  Consequently,  stability 
margins  predicted  from  flap-lag  studies  may  be  lower  than 
f lap- lag-torsional  analyses. 

Calculating  the  stability  regions  of  both  the  flap-lag 
motion  and  the  f lap-lag- torsion  motion  of  a  helicopter  blade 
have  been  accomplished  with  various  procedures.  Several 
studies  have  linearized  the  equations  of  motion  about  a 
periodic  equilibrium  position  to  obtain  a  system  of  linear 
ordinary  differential  equations  with  periodic  coefficients 
and  used  Floquet  theory  to  determine  the  system’s  stability 
( 13 ; 15 ; 16 ; 25 ; 27) .  Xaza  and  Kvaternik  (16)  analyzed  rotor 
flap-lag  instabilities  by  transforming  the  blade’s  equations 
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of  motion  to  a  nonrotating,  apace-fixed  coordinate  system  and 
time  averaging  the  equations'  periodic  coefficients  to  obtain 
an  approximate  system  of  linear  equations  with  constant 
coefficients.  Crimi  (21)  extended  Hill’s  method  to  solve  a 
system  of  second  order  linear  differential  equations  with 
periodic  coefficients.  This  method  yielded  a  set  of  complex 
exponents,  in  Floquet’s  form,  as  an  initial  value  solution, 
and  provided  a  quantitative  measure  of  the  system's 
stability.  Quasi-linearization  was  utilized  by  Panda  and 
Chopra  (26)  to  solve  the  nonlinear  equations  of  motion  with 
an  iterative  procedure  based  on  Floquet  theory.  This  type  of 
solution  contained  all  the  harmonics  for  flap,  lag  and 
torsion  response  amplitudes.  Finally,  harmonic  balancing  was 
used  (15)  to  calculate  a  blade's  steady-state  response  by 
assuming  the  response  was  periodic  and  consisted  of  a  sum  of 
a  finite  number  of  harmonics.  This  procedure,  however,  could 
become  quite  involved  for  coupled  systems  with 
nonlinearities.  Frequently,  though,  accurate  solutions  were 
obtained  when  the  harmonic  balance  method  was  restricted  to 
only  zero  and  first  order  harmonics. 

Friedmann  (27)  provided  an  excellent  review  of  the 
aeroelastic  stability  and  response  problems  associated  with 
isolated  rotor  blades  and  coupled  rotor/f uselage  systems  in 
hover  and  in  forward  flight.  He  summarized  the  stability 
analyses  of  coupled  flap- lag  and  coupled  flap-lag-torsion 
motions  of  blades,  and  discussed  the  importance  of  trim  and 
nonlinear  terms  on  blade  stability.  Further,  he  reviewed  the 


coupled  rotor/f uselage  aeroelastic  problem  and  noted  that 


helicopter  stability  analyses  have  basically  been  restricted 
to  single  blade  or  isolated  blade  aeroelastic  problems. 
Friedmann  noted  that  interblade  mechanical  couplings  or 
couplings  between  the  rotor  and  the  fuselage  could  have 
significant  effects  on  a  helicopter’s  dynamic  stability. 


Coupled  Bo tor /Fusel age  Dynamics 


To  study  the  stability  characteristics  of  a  coupled 
rotor/fuselage  system  a  simple,  yet  accurate,  mathematical 
model  is  required.  Indeed,  this  area  of  study  has  been 
addressed;  several  authors  (1;2;5)  provided  a  brief 


introduction  to  developing  the  equations  of  motion 

of 

the 

entire  helicopter, 

to 

include 

both  the  rotors 

and 

the 

a i r  f  rame .  Bousman 

(28) 

obtained 

experimental  data 

on 

the 

stability  of  a  hingeless  rotor  mounted  on  a  special  gimbaled 
support  which  simulated  body  pitch  and  roll  degrees  of 
freedom.  He  then  compared  this  data  with  theoretical  results 
he  had  obtained  from  an  analytical  model.  Other  studies 
(20-38)  have  introduced  analytical  models  to  examine  coupled 
rotor/f uselage  dynamics.  Further,  a  detailed  description  of 
the  procedure  followed  to  obtain  a  set  of  rotor/fuselage 
equations  of  motion  was  accomplished  by  Venkatesan  and 
Friedmann  (39) .  These  authors  presented  a  set  of  governing 
coupled  differential  equations  for  a  model  representing  a 
Hybrid  Heavy  Lift  Helicopter,  which  used  a  multiple  rotor 
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system  for  its  operation.  They  subsequently  reduced  these 
equations  to  a  single  rotor  model  and  obtained  a  simplified 
system  of  coupled  rotor/fuselage  equations  which  produced 
results  agreeing  favorably  with  experimental  data  (40) .  The 
model  included  a  rotor  with  three  or  more  hinged,  rigid 
blades  and,  in  deriving  the  system's  equations  of  motion, 
neglected  terms  second  order  and  higher. 

Solving  the  equations  of  motion  to  determine  the 
stability  regions  of  a  coupled  rotor/f uselage  system  can  be 
accomplished  with  many  of  the  methods  used  for  the  isolated 
blade  cases.  Because  the  equations  are  so  complex,  industry 
has  developed  a  number  of  complex  analyses  and  have 
implemented  them  in  sophisticated  computer  programs 
(27 ; 35 ; 36) .  It  is  possible,  however,  to  simplify  the  coupled 
rotor/fuselage  equations  of  motion  so  that  relatively 
uncomplicated  procedures  can  be  used  to  determine  the  sytem’s 
stability  (3 1 ; 32 ; 34 ; 37) .  For  example,  Straub  and  Warmbrodt 
(31)  and  Straub  (37)  converted  linearised  perturbation 
equations  with  periodic  coefficients  into  a  constant 
coefficient  system  using  a  Fourier  coordinate  transf  or  station 
so  the  system’s  stability  could  be  evaluated  by  performing  an 
eigenvalue  analysis. 

Helicopter  Control 


Rather  than  Just  documenting  where  the  instabilities  of 
an  isolated  rotor  blade  or  a  coupled  rotor/fuselage  system 


will  occur,  it  would  be  desirable  to  actually  reduce  or 
eliminate  those  instabilities. 

For  example,  Peters  and  Hohenemser  (7)  used  Floquet 
theory  to  establish  the  stability  of  a  rotor  blade's  flapping 
equations  of  motion  with  <5g  (pitch-flap  coupling  angle;  it  is 
positive  if  the  pitch  is  decreased  when  the  blade  flaps  up) 
and  tilt-moment  feedback.  They  found  that  adding  the 
feedback  to  the  linearized  equations  delayed  the  onset  of 
instabilities.  After  linearizing  an  isolated  rotor  blade's 
flap-lag  equations  of  motion,  Johnson  and  Hohenemser  (22) 
described  a  rotor  thrust  or  tilting  moment  feedback  control 
system  which  essentially  eliminated  the  blade's  first  order 
flapping  harmonics.  Even  though  the  system  did  not  eliminate 
higher  blade  flap  bending  harmonics,  it  did  minimize  the 
chordwise  and  flapwise  blade  motions.  Similarly,  Gaffey  (17) 
concluded  that  negative  <5g  feedback  eliminated  flap-lag  blade 
motion  instabilities  since  it  separated  the  flap  and  lag 
natural  frequencies.  However,  this  also  introduced  the 
possibility  of  a  flapping  divergence.  Further,  Gaffey 
discovered  that  the  effect  of  negative  <5g  on  pitch-lag 
instabilities  was  stabilizing. 

Bo tor  blade  control  has  also  been  accomplished  by 
designing  the  blade  as  an  el 1 iptical ly-shaped  circulation 
controlled  airfoil  with  a  simple  cyclic  control  mechanism 
based  on  blade  pressure  variation  rather  than  blade  pitch 
variation  (24) .  A  thin  jet  of  air  blown  from  a  spanwise  slot 
along  the  blade's  rounded  trailing  edge  controlled  blade  lift 


(23).  This  model  delayed  the  onset  of  f lap-lag-torsion 
instabilities  without  any  moving  parts  other  than  the 
rotating  blades. 

Quinn *8  (41)  rotor  control  design  placed  the  control 
actuators,  power  supplies  and  computers  in  the  rotating 
control  system.  Bef erred  to  as  an  ‘Individual  Blade  Control 
Independent  of  a  Swashplate*  (IBIS) ,  this  system  consisted  of 
four  single  actuators  per  blade,  each  controlled  by  a 
different  power  supply  located  at  the  rotor  hub.  Four 
computers,  also  located  at  the  hub,  processed  the  control 
data.  Quinn  claimed  advantages  in  reliability,  drag,  weight, 
and  cost  over  conventional  swashplate  control  systems. 

As  for  controlling  coupled  rotor/f uselage  motion, 
Miyajimi  (20)  developed  a  stability  and  control  augmentation 
system  whereby  the  helicopter  was  represented  by  a  six-degree 
of  freedom  rigid  body  (excluding  the  rotors)  with  constant 
factors  used  as  weightings  for  the  state  and  control 
variables.  A  least  squares  design  method  was  applied  to 
determine  the  control  augmentation  system.  Then,  blade 
flapping  motion  was  included  using  the  method  of  multiblade 
coordinates,  and  linear  optimal  control  theory  was  applied  to 
determine  the  appropriate  feedback  gains  for  the  stability 
augmentation  system. 

Straub  and  Warmbrodt  (31)  ,  after  approximating  the 
coupled  rotor/f uselage  equations  of  motion  with  constant 
coefficients,  used  state  variable  feedback  with  appropriate 
closed- loop  feedback  phase  and  gain  margins.  They  then 
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flapping  stability  throughout  the  flight  regime. 


Their 


controller  did  indeed  eliminate  the  flapping  instabilities  of 
one  and  two  helicopter  blades  at  high  advance  ratios. 
Although  both  studies  (42; 43)  dealt  only  with  blade  flapping 
instabilities  at  high  advance  ratios,  the  technique  utilized 
is  directly  applicable  to  controlling  the  instabilities  of 
more  complex  helicopter  models. 

Implestenting  the  modal  control  technique  has  several 
advantages  over  previous  attempts  to  reduce  or  eliminate 
blade  (17;22;23)  and  rotor/fuselage  (31;37;44)  instabilities. 
For  exasiple,  thrust,  6g,  or  tilt  feedback  control  is  no 
longer  required,  and  conventional  swashplates  can  be  used  in 
the  control  system.  The  linear  periodic  system  is  not 
approximated  by  a  constant  coefficient  system  in  order  to 
apply  state  feedback  control.  Consequently,  active  control 
of  blade  and  rotor/fuselage  motion  using  collective  and 
cyclic  pitch  mechanisms  to  eliminate  the  system’s  unstable 
modes  could  possibly  enhance  a  helicopter’s  stability. 


Problem  Statement 

This  research  effort  is  intended  to  demonstrate  a 
control  technique  which  does  indeed  improve  a  helicopter’s 
stability.  As  previously  mentioned,  the  modal  control 
technique  effectively  eliminated  an  isolated  rotor  blade’s 
flapping  instabilities  using  existing  control  mechanissu. 
Successfully  demonstrating  this  technique  on  more 


1 

Is 


sophisticated  helicopter  models  would  therefore  serve  two 
purposes:  a  helicopter’s  stability  could  indeed  be  improved 
using  modal  control;  and  the  modal  control  technique  would  be 
validated  on  more  realistic  systems. 

The  first  step  in  pursuing  this  course  of  research  is  to 
extend  the  isolated  helicopter  blade  model  from  considering 
only  flapping  motions  (42;43)  to  considering  both  flapping 
and  lagging  motions.  Including  flap-lag  motion  introduces 
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not  only  another  degree  of  freedom  to  the  equations  of  motion 
but  also  adds  coupling  terms  to  those  equations.  This  more 


realistic  model  has  been  shown  to  demonstrate  a  rotor  blade’s 
instability  at  lower  advance  ratios  than  the  simple  flap 
model.  These  lower  advance  ratios  may  easily  fall  within  the 
performance  envelope  of  a  helicopter.  Finally,  extending  the 
modal  control  technique  to  controlling  blade  flap-lag  motion 
will  demonstrate  this  theory  on  a  more  sophisticated  model. 

The  control  of  a  helicopter’s  coupled  rotor  and 
longitudinal  fuselage  motion  will  then  be  considered.  The 
blade  will  again  be  modeled  as  rigid  and  will  have  both  flap 
and  lag  degrees  of  freedom.  Demonstrating  that  the  modal 


control  technique  eliminates  instabilities  in  the  coupled 
rotor/fuselage  dynamics  will  indeed  improve  a  helicopter’s 
stability  in  a  reliable  and  efficient  manner. 


In  summary ,  the  purpose  of  this  research  program  is 
twofold.  First,  the  equations  of  action  for  a  helicopter 
blade's  flap-lag  dynamics  will  be  developed.  The  blade’s 
stability  in  hover  and  in  forward  flight  will  be  examined 
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using  Floquet  theory,  and  the  results  compared  to  previous 
studies.  The  pole  placement  technique  mill  then  be  used  to 
actively  control  the  stability  of  the  helicopter  blade  in 
forward  flight  using  existing  collective  and  cyclic  pitch 
control  mechanisms. 

Once  an  isolated  rotor  blade’s  stability  has  been 
investigated,  the  equations  of  motion  for  the  coupled 
rotor/fuselage  dynamics  mill  be  derived.  Using  Floquet 
theory,  the  stability  of  the  coupled  rotor/fuselage  system 
mill  be  examined.  Finally,  control  of  the  system  mill  be 
accomplished  mith  the  pole  placement  technique,  as  me 1 1  as 
standard  feedback  control  mechanisms. 

The  development  of  the  blade  flap-lag  equations  of 
motion  is  discussed  in  Chapter  II  mhile  Chapter  III 
highlights  the  derivation  of  the  coupled  rotor/fuselage 
equations  of  motion.  The  control  technique  used  to  eliminate 
the  helicopter's  unstable  regions  is  reviewed  in  Chapter  IV. 
The  results  of  controlling  the  blade  flap-lag  motion  is 
presented  in  Chapter  V,  mhile  Chapter  VI  contains  the  results 
of  controlling  the  helicopter’s  coupled  rotor/fuselage 
motion.  Finally,  a  number  of  conclusions  concerning  the 
results  obtained  are  presented  in  Chapter  VII. 
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I I .  Flap-Lag  Stabl 1 1 ty  of  Hel ioopter  Rotor  Blades 

Deriving  a  helicopter's  equations  of  notion  and 
conducting  a  stability  analysis  on  the  system  is  a 
challenging  venture.  The  complexity  of  the  helicopter,  as 
well  as  the  rotor,  requires  using  assumptions  and 
approximations  which  simplify  the  equations  of  motion  and 
make  them  more  manageable.  However,  care  must  be  taken  to 
avoid  oversimplifying  the  problem  to  the  exclusion  of 
important  effects. 

A  logical  starting  point,  in  a  stability  analysis,  is  to 
derive  the  equations  of  motion  of  a  single  rotor  blade  in 
forward  flight  and  then  examine  the  rotor's  stability 
characteristics.  The  more  complex  task  of  deriving  the 
equations  of  motion  of  the  entire  helicopter  stay  then  be 
undertaken  once  this  is  accomplished. 

The  examination  of  the  stability  of  an  isolated 
helicopter  rotor  blade  is  based  on  a  system  of  equations  of 
motion  for  a  rigid  blade.  The  blade  is  centrally  hinged  and 
free  to  both  flap  and  lag.  The  assumptions  and  approach  used 
are  summarized  below. 

Blade  Dynamics 

An  analysis  of  the  various  forces  and  moments  acting  on 
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Figure  3.  Blade  Reference  Frames 


Figure  4.  Hinge  Geometry 


a  helicopter  rotor  in  forward  flight  may  be  conducted  in 
several  different  coordinate  systems.  Three  of  these 
reference  frames  are  quite  useful  in  relating  blade  motion  to 


Y  form  a  dextral  set  with  Z.  The  X  axis  is  along  the 

helicopter’s  velocity  vector  when  it  is  in  forward  flight. 

In  deriving  the  equations  of  motion  for  an  isolated  rotor 

blade,  it  is  assumed  that  this  frame  at  most  translates  with 

uniform  velocity  relative  to  a  fixed  space.  Hence,  can  be 

considered  to  be  an  inertial  system. 

Frame  F.  is  also  attached  to  the  hub  but  rotates  with 
n 

the  blade  azimuth  rate.  The  z  axis  is  coincident  with  the  Z 
axis,  and  the  x  axis,  which  lies  in  the  X-Y  plane,  is  in  the 
direction  of  the  rotor  blade.  The  y  axis  is  defined  so  as  to 
yield  a  dextral  set. 

The  final  reference  frame,  F.  ,  is  a  blade-fixed  frame 

b 

with  the  x  axis  coincident  with  the  blade  feathering  axis. 

When  considering  a  blade’s  flapping  and  lagging  motions  it  is 

necessary  to  specify  the  rotation  sequence.  A  flap-lag 

rotation  is  used  when  the  lag  hinge  is  outboard  of  the  flap 

hinge  [Figure  4(a)]  while  a  lag-flap  rotation  is  used  when 

the  lag  hinge  is  inboard  of  the  flap  hinge  [Figure  4(b)] 

(16:877).  In  Figure  4(a),  frame  F.  is  found  from  F.  by  a  2-3 

b  b 

rotation  through  the  angles  -ft  and  C-  On  the  other  hand. 
Figure  4(b)  indicates  that  Fb  is  found  from  Ffc  by  a  3-2 
rotation  through  the  angles  C  and  -ft.  In  both  situations, 
the  flap  and  lag  hinges  are  assumed  to  be  coincident  and 
located  at  the  rotor  shaft.  The  y  and  y  axes  form  a  dextral 
set  with  %. 

The  equations  of  motion  used  in  this  study  are  developed 


Figure  5.  Free  Stream  Velocity 


assuming  a  flap-lag  hinge  Sequence  (16:877).  In  addition, 
the  derivation  of  the  aerodynamic  forces  is  based  on  linear, 
quasi-steady  strip  theory.  The  induced  inflow  is  uniform, 
and  stall  and  compress ibility  are  not  considered. 

Figure  5  shows  the  relationship  between  the  free  stream 
velocity  vector  and  the  inertial  coordinate  system.  The 
helicopter’s  forward  velocity  has  magnitude  V  at  an  angle  of 
attack  a  to  the  X-T,  or  disk,  plane.  The  rotor  induced 
velocity,  v is  normal  to  the  disk  plane,  while  the  advance 
ratio  m  and  inflow  ratio  X  are  dimensionless  velocity 
components  parallel  to  and  normal  to  the  disk  plane, 
respectively.  These  terms  are  defined  as  follows  (2:25,28): 

V  =  Vcosa/  (BO)  (2.1) 

A.  =  (Vsin  a  -  tO  /  (HO)  (2.2) 

where  B  is  the  blade’s  length. 

Forward  flight  introduces  a  phenomena  known  as  reversed 
flow.  The  region  of  reversed  flow  is  an  area  on  the 
retreating  side  of  the  hub  plane,  near  the  blade  root,  where 


R.ir".if».»  »jt  ■>  vjvwt  wv  wvvi»  11 » t"t**r*jiHj’Rr**rwmnwwiiwwu»ir« 


the  velocity  relative  to  the  blade  is  directed  from  the 
trailing  edge  to  the  leading  edge  of  the  blade.  This  means 
that  in  this  region  the  rotational  velocity  is  smaller  in 
magnitude  than  the  component  of  the  forward  speed,  and  the 
flow  is  reversed  (2:152).  Sissingh  (8:57)  documented  three 
different  types  of  flow  regions  a  blade  encounters:  normal; 
reversed;  and  mixed  flow.  However,  Johnson  (2:152)  reported 
that,  for  low  advance  ratios,  the  reversed  flow  region 
occupies  only  a  small  portion  of  the  rotor  disk,  and  it  is 
characterized  by  low  dynamic  pressure  until  the  advance  ratio 
becomes  large.  Since  the  root  cutout  typically  extends  from 
15  to  30  percent  of  the  rotor  radius,  it  covers  much  of  the 
reversed  flow  region.  Hence,  the  effects  of  the  reversed 
flow  region  are  negligible  up  to  an  advance  ratio  of  -  0.50 
(2:152).  Because  most  helicopter  forward  speeds  represent 
values  of  /j  between  0.30  and  0.40,  reversed  flow  is  not 
considered  when  deriving  the  isolated  rotor  blade’s  equations 
of  motion  (2 : 152 ; 16 : 878) . 

Since  the  rotation  sequence  for  the  rotor  blade  is  flap 
followed  by  lag,  the  coordinate  transformation  relating  the 
blade-fixed  axes  x,  y,  y  to  the  rotating,  hub-fixed  axis 
system  x,  y,  z  is 

'x 
V 

where 


X 

X 

=  LblLlh 

y 

Lbh 

y 

K 

z 

■  m 

z 
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<-  /  /  <" 


where  aft  =  sin/?,  eft  =  cos/?,  s £  =  sinC  ,  and  c£  -  cosC  • 


How,  if  a  and  b  are  any  two  vectors  in  the  same 
frame,  the  components  of  the  vector  product  a  x  b 


found  from 


m  • 


where 


dextral 
may  be 


(2.9) 


In  order  to  determine  the  system's  kinetic  energy,  the 
derivative  of  the  position  vector  as  seen  by  an  observer  in 
the  inertial  frame  is  required.  This  derivative,  expressed 
in  the  blade-fixed  axis  system,  can  be  written  as(45:50): 


-  P)b  *  [“b/I]bBb 

However ,  n  ,  the  derivative  of  the  position  vector 

seen  by  an  observer  fixed  in  the  blade-fixed  reference  frame, 
is  zero.  Hence, 


0 

II 

M 

r( C  *  O cos/?) 

r(/?cosC  “  Osin/?sinC) 

(2.11) 


The  kinetic  energy  of  a  single  blade  is  given  by 


t  =  y  (pv  {pi}„Hdr 


(2.12) 


where  is  the  density  of  the  point  mass. 
The  system’s  potential  energy  is 


(2.13) 


ic 


A 

h 


!■ 


i 


r. 


A 

V 


t  ‘ 


V  = 


=  k^3  /2  ♦  kcC  /2 


with  and  k,.  being  the  spring  rates  about  the  flap  and  lag 


hinges,  respectively.  Vote  that  the  gravity  potential  is 
ignored . 

Since  the  Lagrangian,  L,  is  defined  as  L  =  T  -  V  , 
solving  the  expressions  for  the  system’s  kinetic  energy,  and 
combining  this  result  with  the  potential  energy  of  the  system 
yields  the  following: 


=  i-I  [c2^2cos2C+02  [cos2^sin2/)sin2c] 


^2fty?sin^sinCcosC+2QCeo8^j  ~  ^k^?2  ~ 


(2.14) 


iriiere,  assuming  a  thin  blade,  the  blade’s  mass  moment  of 
inertia  for  both  flap  and  lag  motions  are  equal  and  can  be 
approximated  as 

B 


2 

r  p  dr 


(2.15) 


Meirovitch  (45:72-76,88-91)  derived  Lagrange’s  equations 
for  a  system  with  both  conservative  and  nonconservative 
forces : 

f*l 


d_ 

dt 


*Ii 


-  *k  =  Q 


k  =  1,2 . n 


(2.16) 


where  qfc  are  generalized  coordinates  and  are  generalized. 


non-conservative  forces. 

The  generalized  forces  include  both  the  aerodynamic 
forces  acting  on  the  rotor  blade  and  the  forces  resulting 

23 


c 


from  the  torsional  dampers  situated  along  the  hinge  line. 
Thus , 

Qk  =  “0^9  -  ccC  ♦  Qk  (2.17) 

where  c^  and  are  the  viscous  damping  coefficients  in  flap 

# 

and  lag,  respectively,  and  are  the  generalized  aerodynamic 
f  orces . 

If  the  angles  ft  and  C  are  assumed  to  be  small  and  terms 
through  second  order  are  retained  (16:878),  Lagrange’s 
equations  of  motion  can  be  obtained: 

ft  ♦  2£X ft  +  c^ft/1  +  ^O2  +  k^/lj/?  =  Q^/I  (2.18) 

C  -  2Qftft  ♦  e^C/I  +  k^C/I  ■  Q^/I  (2.10) 

#  # 

where  and  Q^.  are  generalized  aerodynamic  forces. 

The  azimuth  angle,  y/,  is  given  by  y  =  fit  Hence,  the 

equations  of  motion  can  be  differentiated  with  respect  to  the 
azimuth,  resulting  in 

ft"  ♦  2C '  ft  *  2pC ijft'  *  P2ft  =  Q^/(02I)  (2.20) 

c”  -  2ft'  ft  *  2CC^C  +  w2C  =  Q^/(n2I)  (2.21) 

where:  p,  the  dimensionless  rotating  flapping  natural 

“2  —  2  — 2 

frequency,  is  1  ♦  and  ,  the  nondimensional 

nonrotating  flapping  and  lagging  natural  frequencies,  are 
k^/(02I)  and  k^./(f)2I)  ,  respectively;  and  end  C^- .  the 

flap  and  lag  viscous  damping  ratios,  are  c^/(2Iftp)  and 
c^./(2IOu>^.)  ,  respectively. 
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Figure  6.  Aerodynamic  Velocity  Components 


Botor  Aerodynamic  Forces 


To  completely  derive  the  blade  flap-lag  equations  of 


motion,  the  aerodynamic  forces  must  be  determined.  Consider 


the  velocity  components  seen  by  a  rotating  blade  while  the 


helicopter  is  in  forward  flight.  The  helicopter  has  forward 


velocity  V  and  a  blade  section  angle  of  attack  cx.  From 


Figure  6  it  is  apparent  that  the  component  of  the 


helicopter's  velocity  in  the  hub  plane  is  Vcosa  or  ijOR . 


Consequently,  the  flow  relative  to  the  blade  can  be 


expressed  as 


VA  = 


Vcosacosv' 
-Vcosasiny' 
Vsina-v , 


pORcosv' 

-/jftBsiny/ 

\0B 


(2.22) 


where  the  point  A  is  fixed  to  the  bub  but  rotates  through  an 


azimuth  angle  y.  The  induced  velocity,  v.,  is 


2 


c 


V.  =  COR/  [2  [m2*X2]‘"2] 


(2.23) 


and  the  inflow  ratio  can  be  expressed  as 


X  =  ptana  - 


_  2  .  2  1 1/2 

2  m  +X  I 


(2.24) 


with  C  being  the  thrust  coefficient. 

T 


The  expression  for  the  total  velocity  of  a  blade  mass 


element  is  the  sum  of  the  aerodynamic  and  dynamic  components 


(16:878) ,  or 


V  =  V.  ♦  V. 

A  dm 


(2.25) 


where  V 


■  PL  • 


It  is  customary  to  express  this  velocity  in 


dimensionless  radial,  tangential,  and  perpendicular 


components  as  Ur ,  UT ,  and  0p ,  respectively.  Hence,  in  the 


blade-fixed  frame  [following  the  orientation  used  by  Ka2:a  and 


Kvaternik  (16:879)), 


AJORc^cCcw  -  J-ifiBsCsv  +  XflBsfJcC 


Ol  =  rC  +  rO eft  ♦  /jOBc/?sC cf  ♦  pOBcCsv'  +  XOBs^sC  I  (2.26) 


|0pJb  |_r/7cC  ♦  rOs^sC  +  -  XfiRcft 


where  cy  =  cosv'  and  SV'  =  siny/. 


The  total  velocity  of  the  air  mass  is  now  written  in 


terms  of  the  blade-fixed  axis  system.  Further,  the  relative 


angle  of  attack  can  be  easily  defined,  and  blade  element 


theory  can  be  used  to  calculate  the  section  aerodynamic 


forces  (2 : 45-51 ; 3 : 1 ; 46) .  This  theory  assumes  each  blade 


Figure  7.  Blade  Section  Aerodynamics 


K 


section  acts  as  a  two-dimensional  airfoil,  and  the  induced 
velocity  at  the  section  accounts  for  the  influence  of  the 
rotor  wake.  Thus,  two-dimensional  airfoil  characteristics 
can  be  used  to  evaluate  the  section  loads  in  terms  of  the 
blade  motion  and  the  aerodynamic  forces. 

From  Figure  7  the  resultant  velocity,  U,  and  the  inflow 

(2  < ✓? 

0T+r  and  4>  - 

tan  1  fup/U  1  The  blade  section  pitch  6  is  measured  from 


<  V 


K 


4 


the  reference  plane  to  the  zero-lift  line;  it  includes  the 

collective,  cyclic  cosine,  and  cyclic  sine  pitch  angles,  , 

Q  ,  6  ,  as  well  as  the  blade’s  built-in  twist  (2:168-169). 
c  s 

Thus,  the  section  angle  of  attack  can  now  be  expressed  in 
terms  of  the  pitch  and  inflow  angles  as 


a  -  6  -  <p  -  6  -6  cosy/  -  6  siny/  -  tan 

O  C  8 


(2.27) 


For  an  elemental  section  of  length,  dr,  the  lift  and 
drag  forces  can  be  expressed  in  terms  of  elemental  components 
as  (5:130-135:16:879) 
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dL  =  ^pU  cC ^dr 


(2.28) 


dD  =  ipU2cCddr  (2.29) 

where  c  is  the  section  chord  and  p  is  the  air  mass  density. 

The  lift  and  drag  coefficients,  and  ,  are,  in 

general ,  functions  of  the  section  angle  of  attack  a . 
Defining  the  two-dimensional  lift  curve  slope  as  a  =  C^/a 
Eqs  (2.28)  and  (2.29)  can  be  rewritten  as 


dL  =  ^pactJ2adr  (2.30) 

dD  =  ipacD2  ^Cd/ajdr  (2.31) 

The  components  of  these  aerodynamic  forces  resolved  into 

the  blade- fixed  reference  frame  are  F  ,  F  ,  and  F  (radial, 

x  y  9 

in-plane,  and  normal  forces,  respectively)  (2:170): 


-dLsin#  -  dDcosP 
dLcos0  -  dDsin<£> 


(2.32) 


where  the  radial  component  can  be  neglected  (16:879). 

Combining  Eqs  (2.30),  (2.31),  and  (2.32)  results  in  the 
following  equations  describing  the  elemental  components  of 
the  aerodynamic  forces,  expressed  in  the  blade  axis  system: 


fdF  1 

r-aUU  - 

(C  /a)  UTJ  1 

V 

1 

* 

p 

d  T 

dF 

<6 

aUTJ  - 

(C  ./a)  UU 

L  9s 

L  T 

d  p  J 

(2.33) 


Expressing  these  forces  in  the  inertial  reference  system 


yields 


-coa^ainrdF  -  aintfdF 
V  9 

coaCdF^ 

-ain/?ainCdF^  +  coa^3dF^ 


(2.34) 


Uaing  the  principle  of  virtual  work,  the  generalized 
forced  in  Eq  (2.17)  that  are  aaaociated  with  the  above 
aerodynamic  forced  have  the  form  ( 16 : 879 ; 48 : 60 ,72-79) 

-v =  [  ^UpV  *Fv-Uplk  iF>UpV 


<  -  [  dF*^{pI}b*  dF.-ypIV  dV^pI}J 


(2.36) 


By  taking  the  partial  derivative  of  the  velocity  vector 
p1  with  redpect  to  ft  and  C .  and  integrating  the  elemental 
aerodynamic  forced  over  the  length  of  the  rotor  blade,  the 

t  * 

generalized  aerodynamic  force8  and  may  be  obtained: 


w 


rcoaCdF 


Qc  =  |  rdF. 


(2.37) 


(2.38) 


Thua ,  the  blade  flap-lag  equationa  of  motion  can  be 


written  aa 


ft"*2C  ‘ft*2pt:^  *p2ft  =  Jr^o.UUT- [cd/aJuOpJcoaCdr  (2.39) 


C  -2ft  ft*2K 


u 

cVnc  ■  ft!  |'H^(cd/*)’™T]' 

-  A-l  o 


(2.40) 


Blade  Trim  and  Perturbation  Equations 


A  helicopter  blade's  stability  can  be  evaluated  by 
establishing  its  steady-state  equilibrium,  or  trim,  solution 
and  then  perturbing  the  blade’s  nonlinear  equations  of  motion 
about  the  trim  values  (15:5-8).  The  steady-state  flap  and 
lag  motions  are  described  by  a  Fourier  series  as: 


ft(yj)  =  ft  *  ft.  cosy/  +  ft .  8iny/  ♦  H.O.T. 
O  lc  18 


C  (y/)  =  C  ♦  C,  cosy/  ♦  c,  siny<  ♦  H.O.T. 
O  lc  18 


(2.41) 


(2.42) 


The  trim  solution  sought  is  not  a  constant  but  a 

periodic  solution.  Eqs  (2.41)  and  (2.42)  are  used  to 

construct  this  'trim'  solution.  Usually,  harmonics  above  the 

first  order  are  very  small  and  may  be  neglected  (2:157).  The 

ft  term  is  a  coning  motion  about  the  hub,  while  the  harmonics 
o 

and  ft^g  generate  once-per-revolution  variations  of  the 
flap  angle.  CQ  is  the  blade’s  mean  lag  angle  relative  to  the 
rotor  hub.  The  harmonic  C^c  produces  a  lateral  shift  of  the 
blade  in  the  plane  of  rotation  and  a  lateral  shift  in  the 
rotor’s  center  of  gravity.  Similarly,  C^8  produces  a 
longitudinal  shift  of  the  blade  and  a  longitudinal  shift  in 
the  rotor’s  center  of  gravity  (2:156-158). 

The  blade *8  pitch  motion  can  also  be  represented  as  a 


Fourier  aeries: 
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s 

r. 

K 

< 

•r. 

If 
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\r 


V 


«  ■ 


6(vy)  =  ©  -6  cosv  -  0  siny/ 

O  C  8 


(2 . 43) 

where  the  harmonics  higher  than  first  order  are  again 


neglected.  The  zeroth  harmonic,  ©q ,  is  the  average  blade 


pitch  and  is  called  the  collective  pitch.  The  once-per- 


revolution  variations  in  the  pitch  angle,  0  and  0  ,  are 

C  8 


called  the  cyclic  pitch  angles.  Collective  pitch  controls 
the  rotor  thrust  magnitude  while  cyclic  pitch  controls  the 


thrust  vector  orientation  ( 0  controls  the  lateral 

c 


orientation  while  9  controls  the  longitudinal  orientation) 

s 


(2: 158- 159 ; 3) 


Denoting  the  steady-state  values  of  ft  and  C  as  ftQ  and  CQ 


and  the  perturbations  from  equilibrium  as  Aft  and  AC,  then 

(2.44) 


ft  =  ft  +  Aft 

o 


C  =  co  *  AC 


If  the  above  equations  are  substituted  into  the 
nonlinear  equations  of  motion  [Eqs  (2.20)  and  (2.21)],  and  if 
the  perturbations  are  assumed  to  be  small  enough  to  neglect 
terms  that  are  second  order  and  higher,  two  sets  of  equations 
result : 


p\  =  %  ' 
o 

(*) 

(2.45) 

w*C  =  Qr  / 
C  °  Co 

(2.40) 

99 

Aft 

+  2ft  At;' 
o 

♦  2pC  pAft 

2.  . 

♦  P  A/? 

*  AV 

[o20 

(2.47) 

( 

09 

AC  - 

2ft  Aft' 

o 

*  2CcwcAC# 

-2  . 

♦  u^A C 

=  AS./ 

[^0 

(2.48) 
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Eqs  (2.45)  and  (2.46)  are  the  flap  and  lag  equations 


defining  the  trim  conditions  ft q  and  CQ  while  Eqs  (2.47)  and 


(2.48)  are  the  linear  perturbation  equations  involving  A/7  and 
AC  (16:879-080). 

Now,  trim  requires  force  and  moment  equilibrium  on  the 
blade.  The  rotor  is  maintained  at  a  fixed  value  of  thrust 


coefficient,  C  .  This  coefficient  can  be  defined  as 

T 


CT  =  2  4 

pnCi  E 


(2.49) 


to  be  zero,  C  Is  neglected,  and 
o 


much  less  than  one 


to  be 

small  , 

the 

and  C,  ) 

IS 

are  assumed 

is  assumed  to 

V 

be 

=  [DP  /UT 
V  O  i 

J  J' 

The 

in  Eq 

(2.41)  , 

are 

Still  considered. 

With  a  constant  thrust  coefficient,  vertical  and 
horizontal  force  equilibrium  (refer  to  Figure  6)  can  be 
expressed  as 

D  ♦  H  cosot  +  T  sinot  =  0  (2.50) 

p  o  o 


T  cosot  =  mg 
o 


(2.51) 


where  the  aerodynamic  forces  are  as  follows: 


T  ,  the  thrust,  is  normal  to  the  hub-plane  and  is  given  by 
o 


in  R 


h  l  j 


dFz  dv> 


(2.52) 


0  0 


Hq  ,  the  rotor  drag  force  in  the  hub-plane,  is 

2  n  R 

Ho  =  h  {  J  [dFX  coav'  “  dFY  «inV']d*' 


(2.53) 


0  0 


Dp ,  the  helicopter  parasite  drag,  is 

D,  "  ±PV*f 


(2.54) 


where  f  is  the  parasite  drag  area,  N  is  the  number  of 

helicopter  blades,  m  is  the  rotor  blade's  mass,  and  dF^  , 

o 

dFy  ,  and  dF^  are  the  steady-state  values  obtained  from  Eq 
o  o 

(2.34) . 

The  steady-state  values  for  the  thrust  and  rotor  drag 
forces  can  now  be  calculated: 


To  ■  “-^[*o(!  *  *’)  *)] 

_  MpacR3n2  fCJ  .  ^  1-  -  1 

o  4  La  oj  ?  O  c  2  sj 


(2.55) 


(2.56) 


Defining  the  blade  Lock  number  as  y  =  pacB  /I  and  the 

rotor  solidity  to  be  a  =  Mc/(rcB)  ,  the  thrust  coefficient, 

C  ,  and  H  become 
t  o 


-  T  «.(§  *  "*)  -  *  x[‘  *  H] 

ho  =  r-Hpl  pgd  -  xeo)  .  ^oec  ♦  ±xe 


^d  -  \e 


]  ♦  e  *  i\0 

>J  3  o  c  2  s 


(2.57) 


(2.58) 


Trim  also  requires  rotor  pitching  and  rolling  moment 


e 


equilibrium;  the  rolling  moment.  If  ,  and  pitching 

xo 

M  ,  on  the  rotor  hub  are: 
yo 

2n  B 


>ment , 


M  =  I  I  reinwdF  d 

xo  2n  J  J  K 


0  0 


2n  H 


V  =  '  s?  J  I 


(2.59) 


(2.60) 


0  0 


where  dF^  is  the  steady-state  normal  blade  axis  component  of 

the  aerodynamic  forces  [Eq  (2.32)].  Expanding  these  two 
moment  equations  yields 


\  ■  ^[§*v  •.&  *  Wl  *  H 


(2.61) 


=  +  |^2j  ♦  [i  ♦  |d][|^Vl8  ♦  !^o]j  (2.1 


However,  in  equilibrium.  If  =  M  =0.  Consequently, 

Xo  yo 


expressions  for 

the  cyclic  pitch. 

6  and 
c 

V 

can  be 

derived 

from  Eqs 

(2.61) 

and 

(2.62)  such 

that 

the 

f  irst 

harmonic 

f lapping 

coefficients 

'’ic  »n<1  Pl. 

are 

suppressed 

(16:880) . 

These  are 

1  - 

e  = 

c 

i 

“L- 

0 

(x  . 

r) 

(2.63) 

e  = 

a 

5-. 

*  5^(1  *  I4) 

/&* 

3  2^1 

8^  J 

(2.64) 

The 

trim 

equations  [Eqs  (2 

.45) 

and 

(2.46)  ] 

require 

computing  the  equilibrium  generalized  aerodynamic  forces  Q 
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(2.37)  and 


(2.38)  . 


these  equilibrium 
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and  Q  .  Using  Eqs 
^  o 

forces  are 

B 


r*o  =  |  rcosCdFio=  -  5*'®.* 


(2.65) 


R 

'  I  ’ 


rdF 


rQ  i 


4a 


©  -4x1  -  |x^e  -  ©  -  4 

o  *  o  3  J  4  H  s  6™  o  c  2 

[l+p2]  +  sinV'{-^X©o  -  ix©s  +  g^Xec  ♦Iv1} 

*  =o.^V.  -  5X®c  *  W’ohK)}] 


(2.66) 


The  complete  set  of  trim  conditions  have  now  been 

defined.  Bemembering  that  C  •  C,  .  and  C,  are  all  assumed 

O  lc  18 

to  be  zero  and  and  are  suppressed,  the  various  trim 

equations  can  be  used  to  solve  for  the  remaining  unknowns , 
fiQ,  ©Q ,  ©c ,  ©a ,  and  X.  The  equations  required  to  determine 
these  unknowns  can  be  summarized  as  follows  (noting  that,  due 
to  the  assumptions  used,  several  of  the  equations  derived 
earlier  are  either  redundant  or  identically  satisfied) : 
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(2.69) 
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+  0^  J 


(2.70) 


(2.71) 


independent:  /j  ;  X;  o ;  C  ;  and  C. 

t  d 


Mote  that  the  following  equilibrium  parameters  are 

For  this  analysis  all  of 
these  variables,  except  fj,  are  constrained  to  specific  values 
so  that  the  trim  conditions  at  various  forward  velocities  may 
be  calculated. 

The  steady-state  equilibrium  terms  appearing  in  the 
perturbation  equations  of  motion,  Eqs  (2.47)  and  (2.48),  have 
now  been  determined.  However,  before  these  equations  can  be 

/  9 

solved,  the  perturbed  aerodynamic  forces,  AQ^  and  AQ^.  ,  must 
be  evaluated. 

Perturbing  the  generalized  aerodynamic  forces  [Eqs 
(2.37)  and  (2.38)]  by  AQ^  and  AQ^.  yields 

B 


I.-  It 


A«„  = 


rcosC  A(dF  )  -  rsinC  A£dF 
o  9  ° 


B 


iQc  ■  I 


rA(dFy) 


(2.72) 


(2.73) 


where  the  perturbed  values  of  the  elemental  forces  dF  and 

4F 


dF^  are  obtained  by 


i  Taylor  series  expansion  about  the 
steady-state  equilibrium  condition  (16:880): 


dF  =  dF  +A(dF  ) 

V  vQ  v 
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(2.74) 


dF  =  dF  +  A(dF  ) 

»  *Q  * 

d(dF  )  d(dF  ) 

=  dF  +  9  AU  +  ^  9  AO  ♦  H.O.T. 

3  do  p  du  t 

O  P  O  TO 


(2.75) 


Neglecting  the  higher  order  terms,  assuming  (C^/a)  «  1 

£or  j^C^/a  -  1 J  =  - 1  j  ,  and  expanding  the  partial  derivatives 
in  the  above  equations  leads  to  the  following  expressions  for 
the  perturbed  values  of  the  elemental  aerodynamic  forces 
(16:880) 

A(dFy)  =  -ipac^0UT  -  2Up  }A0p+{<90p  +  2[|d]uT  ^Aojdr  (2.76) 

A  (dF  )  =  ipac[-UT  AUp+{20OT  -  0p  ^Aoldr  (2.77) 

*  *-  O  O  O'  J 


where  6  is  defined  by  Eq  (2.43). 

Notice  that  Eqs  (2.76)  and  (2.77)  contain  steady-state 
and  perturbed  values  for  the  velocity  components  Up  and  UT. 
Evaluating  the  expressions  in  Eq  (2.26)  at  the  trim  condition 
yields  the  steady-state  equilibrium  values 


0  =  Rnf-X  +  fjft  coSV'l  +  rOC  ft 

p0  l  °  J  °  ° 


(2.78) 


0  =  rfi  +  B^ofsinV'  +  C  coSV'l  ♦  BXOft  C  (2.79) 
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and  perturbing  the  expressions  about  the  equilibrium  yields 
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ADp  =  r ClAfi'  *  rn[f?oAC  +  C Q ♦  RflA/?  [\f?o  +  JJcosy'j  (2.00) 

AU  =  rOAC*  -  r  Qf3  Af3  +  RXO  [ft  AC  +  C  A/?l 
T  O  (.  °  °  J 


♦  R/jOIaCcosv'  ~  C  ACainv'  -  ft  C  A/?cosv'] 
C  o  o  o  j 


(2.81) 


With  the  expre88ions  for  the  perturbed  values  of  the 
aerodynamic  forces  now  defined,  the  perturbed  blade  flap-lag 
equations  of  motion  [Eqs  (2.47)  and  (2.48)1  have  been  formed. 
These  equations  are  ordinary  linear  differential  equations 
with  periodic  coefficients. 
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Veri f ication  of  the  Blade  Flap-Lag  Equations  of  Motion 


For  convenience,  the  equations  of  motion  can  be  written 
in  the  form 

x  (v>)  =  ACwlxCy)  (2.82) 


where  A(v/)  is  periodic  with  period  2n  and  the  state  vector 
x(w)  is  defined  to  be: 


x  (y) 


AC. 


(2.83) 


Hence,  the  periodic  matrix  A(v')  of  Eq  (2.82)  has  the 

f  orm: 


0  0  10 

0  0  0  1 

a31  a32  a33  *34 

*41  a42  *43  a44 


(2.84) 
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How  that  the  periodic  state  matrix  A(yO  is  known,  blade 
stability  regions  can  be  determined.  The  stability 

boundaries  generated  are  based  on  the  assumption  that  the 
blade  has  damping  along  the  hinge  line.  The  validity  of  the 
current  model  to  accurately  predict  stability  regions  is 
verified  by  comparing  the  results  with  those  obtained  from 
published  studies  (IS; 16).  In  generating  numerical 

solutions,  Kaza  and  Xvaternik  (16)  made  several  assumptions ; 
these  assumptions  are  also  used  for  this  study.  For  example, 
the  trim  angle  of  attack,  a,  the  parasite  drag  coefficient, 
C  ,  and  the  horizontal  force  coefficient,  C  ,  are  all  set  to 

Q  « 

p 

zero.  In  addition,  the  following  terms  are  set  to  fixed 
values : 


Stability  boundaries  are  obtained  using  Floquet  theory. 
The  state  transition  matrix  is  determined  numerically  and 
evaluated  at  the  end  of  one  period.  The  eigenvalues  of  this 
monodromy  matrix  are  the  system's  characteristic  multipliers. 
System  stability  is  assured  if  all  of  these  characteristic 
multipliers  have  a  magnitude  of  less  than  one. 

Several  different  methods  can  be  used  to  numerically 
integrate  an  isolated  rotor  blade’s  equations  of  motion. 
Qaonkar ,  Simha  Prasad,  and  Sastry  (47)  performed  a  stability 
analysis  of  a  rotor  modeled  as  a  system  with  periodic 
coefficients  and  presented  computer-generated  data  on  the 
comparative  efficiency  of  four  classes  of  solution 
procedures:  Bunge-Kutta  one-step;  Hamming's  predictor 
corrector;  Bui irsch-Stoer  extrapolation;  and  a  hybrid  or 
variable-order  method  which  embodied  the  special  features  of 
one-step  and  multistep  methods,  such  as  the  Gear  type  and  the 
Shampine-Gordon  type.  It  was  determined  that  Hamming’s 
fourth-order  predictor  corrector  method  was  the  most 
economical  with  respect  to  three-digit  accuracy.  This  method 
evaluated  the  equations  of  motion  twice  per  integration  step, 
as  compared  to  four  or  more  times  when  other  methods,  such  as 
Bunge-Kutta,  were  used  (7 : 27 ; 48 : 100- 101 ) .  Consequently,  for 
this  study,  Hamming’s  modified  predictor  corrector  srethod  is 
used  to  numerically  integrate  the  equations  of  motion  for  one 


period . 


The  curves  presented  in  Figures  8  and  8  are  stability 
boundaries  as  a  function  of  p,  the  rotating  flapping  natural 


(. 


frequency,  and  ,  the  nonrotating  lagging  natural  frequency. 

Figure  8,  which  illustrates  blade  flap-lag  stability  in 

hover,  compares  quite  well  with  the  results  presented  in  Kaza 

and  Kvaternik  (16:881).  The  general  size  and  shape  of  the 

curves  are  similar,  with  small  differences  resulting  from  the 

fact  that  Kaza  and  Kvaternik  set  the  value  of  Q  to  be  the 

o 

critical  collective  pitch  angle,  above  which  the  blade 

becomes  unstable.  They  determined  this  value  to  be  0.40.  On 

the  other  hand,  the  present  study  used  Eq  (2.69)  to  calculate 

the  trim  value  for  6  ,  which  turned  out  to  be  0.29722.  In 

o 

addition,  Kaza  and  Kvaternik  did  not  specifically  state  that 

the  steady-state  mean  lag  angle,  C  .  was  zero  when  they 

o 

generated  the  p  versus  curve  for  hover.  Because  C  was 

(.  o 

assumed  to  be  zero  for  the  remainder  of  the  curves  they 

presented,  the  current  research  effort  used  the  same 

assumption ,  even  for  the  hover  case.  If  C  were  indeed  a 

o 

value  other  than  zero,  the  two  curves  would  differ.  Still, 
both  sets  of  data  appear  to  be  in  close  enough  agreement  to 
assume  that  the  present  model's  results  are  valid  in  hover. 

Although  Kaza  and  Kvaternik  (16)  presented  a  stability 
curve  for  a  flap-lag  hinge  sequence  at  M  =  0.40,  this  curve 
was  constructed  using  an  approximate  method  of  solution. 
Even  so,  the  authors  concluded  that  there  was  close  agreement 
between  the  results  using  both  the  approximate  method  and 
Floquet  theory  for  (J  =  0.40  for  the  lag-flap  hinge  sequence. 
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BLADE  INPLANE  FREQUENCY 


BLADE  INPLANE  FF 


Figure  9.  Blade  Flap-Lag  Stability  Boundaries  at  p  =  0 

44 


The  results  presented  in  Figure  9  are  in  very  close  agreement 
with  Kaza  and  Kvaternik's  results  using  the  approximate 
method.  Hence,  it  can  be  concluded  that  the  present  study's 
results  for  p  =  0.40  accurately  compare  with  previous 
studies . 

Because  the  solution  to  the  blade  flap-lag  equations  of 
motion  was  found  to  be  well-behaved  throughout  the  region  of 
0.00  £  p  <  0.50  it  was  assumed  that  the  derivation  of  the 
equations,  as  well  as  the  assumptions  used,  were  consistent 
and  compared  well  with  previous  efforts.  Therefore,  it  was 
concluded  that  the  solution  to  the  rotor  blade's  equations  of 
motion  were  verified  and  considered  accurate  enough  to  be 
used  throughout  the  remainder  of  this  study. 
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III.  Stability  of  a  Coupled  Bo tor /Fusel age  System 

Chapter  II  summarized  the  development  of  the  equations 
of  motion  for  an  isolated  blade  experiencing  flap-lag  motion. 
With  this  information  available,  consider  now  the  equations 
of  motion  for  the  entire  helicopter.  These  equations  are 
derived  in  order  to  investigate  the  stability  of  the  coupled 
rotor/rigid  body  equations  of  motion.  The  assumptions  and 
approximations  used  are  summarized  below. 


Dynami cs 

In  free  flight,  the  helicopter  has  six  rigid  body 
degrees  of  freedom;  namely,  three  rotation  (yaw,  pitch,  and 
roll),  and  three  position  (vertical,  lateral,  and 
longitudinal)  components .  It  is  customary  to  split  the  six 
degrees  of  freedom  of  a  helicopter  into  two  groups  of  three 
each  and  assume  their  dynamics  can  be  analyzed  separately 
(2:822;49).  The  first  group  contains  the  lateral 
dynamics--lateral  velocity,  roll  attitude,  and  yaw 
rate--while  the  other  consists  of  the  longitudinal 
dynamics--longi tudinal  velocity,  pitch  attitude,  and  vertical 
velocity . 

Since  the  rotor  is  a  major  factor  in  the  helicopter’s 
stability  analysis,  the  rotor’s  motion  must  also  be  modeled. 
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For  a  rigid  hinged  rotor,  this  requires  three  additional 


degrees  of  freedom  for  each  blade;  namely,  flap,  lag,  and 
pitch  components.  For  this  study,  pitch  is  used  for  control 
while  flap  and  lag  are  free  variables. 

As  with  analyzing  the  isolated  blade’s  equations  of 
motion,  an  analysis  of  the  helicopter’s  equations  of  motion 
requires  the  establishment  of  several  sets  of  coordinate 
frames.  The  inertial  reference  frame,  X,  T,  Z,  is  fixed  to  a 
flat,  non-rotating  earth.  Further,  the  wind  axes  are  defined 
with  the  x^  axis  being  along  the  free  stream  velocity,  the 
axis  in  the  plane  of  symmetry  when  one  exists  and  nominally 
pointing  down,  and  y^  forming  a  dextral  set.  The  free  stream 
velocity  is  then  expressed  in  the  wind  axes  as 

H,- 

Another  reference  frame  used,  the  vehicle  carried  frame, 
Fy,  is  the  same  as  the  inertial  frame  when  the  earth  is 
assumed  to  be  flat  and  non-rotating.  This  axis  system  is 
denoted  as  xy,  y^,  and  z^.  The  helicopter’s  body  axes  are 
related  to  those  of  F^  via  a  set  of  Euler  angle  rotations 
through  the  angles  p  (yaw)  ,  $  (pitch)  ,  and  p  (roll)  .  This 
set  of  three  independent  rotations  involves  first  rotating  an 
amount  p  about  z^  to  produce  an  intermediate  set  of  axes  x^, 
yj,  Zj,  followed  by  a  rotation  through  an  angle  S  about  y^, 
carrying  the  axes  to  a  second  axis  set  x^ ,  y,  ,  z^,  and 
finally  a  rotation  about  x 2  through  an  angle  p  to  yield  the 
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body  axes  x  ,  y  ,  z  . 
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(3.3) 


where  c p  =  cos??,  ap  -  sin??,  c#  =  cos#,  s#  =  sin#,  cp  =  cosy, 
and  ap  =  siny. 

The  wind  axis  system  may  be  oriented  with  respect  to  the 
body  axes  by  means  of  the  aerodynamic  angles  ot  and  ft: 


_B/W  *  —  ~  _ 

where  co  =  -  ftk ■  and 
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coaftaincn 


-ainftaina 

coaft 

-sinftatna 


-sina 

0 

cosa 


(3.4) 


(3.5) 


Mote  that  a,  the  helicopter’s  angle  of  attack,  and  ft, 
the  sideslip  angle,  can  be  written  as  (50:10-11): 

a  =  tan  1  |w/u^  ;  ft  =  sin_1|v/vj  (3.6) 

where  u,  v,  and  w  are  the  components  of  the  velocity  along 
the  body  axes. 

To  accurately  model  a  helicopter,  its  rotors  must  be 
analyzed.  The  helicopter  considered  possesses  a  single  main 
rotor,  as  well  as  a  tail  rotor.  The  stain  rotor  is  assumed  to 
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be  directly  over  the  helicopter’s  center  of  mass  and  provides 
lift,  propulsive  force,  and  roll,  pitch,  and  vertical  control 
for  the  helicopter  (2:9).  The  tail  rotor  balances  the  sain 
rotor  torque  as  well  as  provides  the  helicopter’s  yaw  control 
(2 : 9 , 264 ; 5 : 93) .  When  the  main  rotor  is  attached  to  the 
helicopter’s  body,  the  hub-fixed  rotating  axis  system,  used 
in  deriving  an  isolated  rotor  blade’s  equations  of  motion, 
can  be  oriented  with  respect  to  the  helicopter  body  axes  by  a 

shaft  angle  of  attack  a  (positive  for  a  rearward  tilt  from 

not 

the  horizontal)  and  a  shaft  cant  angle  <p  (positive  to  the 

flft 

right)  (51:167).  Denoting  the  subscript  m  for  the  stain 
rotor,  this  orientation  can  be  expressed  with  the  following 
matrix : 


-ca 

0 

sa 

m 

m 

8a  8 <p 

c  <t> 

ca  s  <p 

in  in 

m 

in  in 

-aa  c <p 

m  m 

B<p 

m 

-ca  c <p 

m  nH 

(3.7) 


In  a  similar  aianner ,  the  orientation  of  a  tail  rotor  can 
be  specified  by  a  shaft  cant  angle  <p  (positive  for  an  upward 

v 

tilt,  as  sieasured  from  the  vertical)  and  a  shaft  angle  of 
attack  a  (positive  rearward) .  Note  that  the  subscript  t 

v 

denotes  the  tail  rotor.  In  addition,  the  definition  of  the 
tail  rotor  hub  axes  depends  on  the  main  rotor’s  rotation 
direction.  A  counterclockwise  rotating  main  rotor  requires 
the  tail  rotor’s  thrust  to  be  in  the  direction  of  y^  while 
clockwise  rotation  requires  the  thrust  to  be  in  the  opposite 
direction  (51:187).  If  A^=  1  represents  a  counterclockwise 
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rotating  main  rotor  and  A  =  -1  denotes  a  clockwise  rotating 

B 

main  rotor,  the  tail  rotor’s  hub  axes  can  be  related  to  the 
body  axes  by 


(3.8) 
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Assuming  neither  the  main  nor  the  tail  rotors  have  shaft 

angles  of  attack  and  cant  angles,  and  the  main  rotor  rotates 

in  a  counterclockwise  direction,  L  and  L  are  as  follows: 

Bm  Bt 
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(3.8) 


Since  the  earth  is  assumed  to  be  flat  and  nonrotating, 

the  angular  velocity  of  the  body  axis  System  with  respect  to 

— B/I  — B/V  — B/V 

the  inertial  frame  is  u  =  u>  ,  where  to  is  the  angular 
velocity  of  the  body  frame  with  respect  to  the  vehicle 
carried  frame.  Denoting  the  body  axis  components  of  the 
body’s  angular  velocity  with  respect  to  inertial  space  as  p, 
q,  and  r,  Eq  (3.2)  can  be  written  as 
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Eq  (3.11)  expresses  the  kinematics  in  the  body  axis 
system. 

To  aid  in  deriving  the  helicopter's  equations  of  motion, 
the  following  terms  need  to  be  defined:  B  denotes  the  entire 
helicopter,  to  include  both  the  main  and  tail  rotor; 
denotes  only  the  helicopter's  airframe  and  horizontal  tail 
(the  main  and  tail  rotors  are  excluded) .  It  should  be  noted 
that  the  airframe  includes  not  only  the  helicopter’s  body  but 
also  the  vertical  tail  assembly,  to  which  the  tail  rotor  and 
horizontal  tail  are  attached. 

The  equations  of  motion  for  the  helicopter  are 


=  (M  ) 

cm  ext 


(3.12) 


where  F  the  total  external  force  acting  on  the  body,  is 


Fext  Faero  *  FMR  +  FTR  + 


(3.13) 


F  represents  the  aeronautical  forces  acting  on  B. ,  F._ 

aero  1  MR 

and  F_  include  the  external  forces  acting  on  the  main  and 
rn 

tail  rotors,  respectively,  and  MgA^  is  the  weight  force.  In 
addition,  M  is  the  helicopter’s  total  mass,  including  the 

main  and  tail  rotors,  and  a  ,  the  inertial  acceleration  of 

cm 

the  helicopter’s  mass  center,  can  be  expressed  as 


cm  dt  dt 


(3.14) 


.  -W/I  -W/B  -B/I 

where  to  -  to  ♦  to 


is  the  inertial  rate  of  change  of  the  helicopter’s 


angular  ■omentum  about  its  mass  center,  and  includes  the 


9 


angular  momentum  of  B 


,  the  main  rotor 


and  the 


tail  rotor  (*JJ.  Finally.  <“cm>ext.  the  total  external 
moment  acting  on  the  helicopter,  includes  the  terms  **aoro  + 
“mb  *  Sr  ,  where  **aero  •  “mb  ,  and  “tr  are  the  moments  acting 
on  B^,  the  main  rotor,  and  the  tail  rotor,  respectively. 

Flap  and  lag  contributions  to  the  equations  of  motion 
for  B j  should  be  considered  in  this  analysis.  The 
contributions  to  the  force  equations  are  of  the  order  of  m/M; 
hence,  they  can  be  neglected.  Each  blade's  flap  and  lag 
motions  do,  however,  produce  coupling  effects  on  the 
airframe’s  moment  equations.  For  the  present,  the  moment 
equations  will  be  derived  ignoring  any  contributions  from  the 
flap  and  lag  motions.  Later  in  the  chapter,  the  coupling 
between  the  body  and  the  blades  will  be  addressed. 

Figure  10  summarizes  the  forces  and  moments  acting  on 
the  helicopter  (2:777).  The  notation  used  in  this  figure  is 


as  follows:  x  ,  y  ,  z  are  the  body  axes;  t  ,  t  ,  t  are  the 
Jd  b  d  x  y  z 

distances  from  the  center  of  mass  (c.m.)  to  the  tail  rotor’s 

hub;  l  is  the  vertical  distance  from  the  c.m.  to  the  main 
z 

rotor’s  hub;  T.,  T  are  the  distances  from  the  c.m.  to  the 
x  z 

horizontal  tail;  T  ,  H  ,  Y  ,  D  are  the  main  rotor’s  thrust, 

m  m  m  p 

longitudinal  (drag)  force,  side  force,  and  parasite  drag 
force,  respectively;  T. ,  H. ,  Y.  are  the  tail  rotor’s  thrust, 
longitudinal  force,  and  side  force,  respectively;  M^,  M  , 


M  are  the  B  ,  main  rotor,  and  tail  rotor  pitching  moments, 
yt  1 
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Figure  10.  Helicopter  Forces  and  Moments 


respectively;  ,  DT  are  the  horizontal  tail’s  lift  and 

forces;  and  L  .  D  ,  Y  are  the  lift,  drag,  and  side 
B  B  B 

acting  on  B ^  through  the  helicopter’s  center  of  mass. 
How,  in  matrix  form,  Eq  (3.14)  can  be  written  as 


Expressed  in  body  axes,  Eq  (3.15)  becomes 


From  Figure  10,  the  external  forces  are 


■D  -D 
B  T 


{'*•"},'  Ib  :  {f“}b= 


H  +D  cosa 
m  p 


T  +D  sina 

■  m  p 


H. 


drag 

forces 


(3. 15) 


(3.16) 


(3.17) 


Collecting  terms,  Eq  (3.12)  yields  the  force  equation 


LB*  ° 


lbv°  -[“W/I]blbw  0 

V  L®  J  W 


(3.18) 


Developing  the  moment  equations  for  the  entire 
helicopter  requires  calculating  the  contributions  from  as 

well  as  the  main  and  tail  rotors.  The  inertial  rate  of 
change  of  the  helicopter  body's  (Bj)  angular  momentum  about 
the  mass  center  is  expressed  as 

■  *5  *  P'IbWb 

where  |®B^  =  J  How,  the  moment  of  inertia, 

I  .is  about  the  helicopter's  center  of  gravity  and 

L  CJB 

includes  not  only  the  body  terms  but  also  the  terms  for  the 
main  and  tail  rotors: 


Mb  ■  Mb  *  ML  *  Ml 


(3.20) 


Assuming  the  helicopter’s  airframe,  excluding  the 
rotors,  is  symmetrical  in  the  x  s  plane,  the  moment  of 

D  O 


inertia  for  B^  is  given  by 


rwriUTwn 


e 


r 


f 


( 


Here,  B  is  the  main  rotor’s  radius,  r  is  the  radial 

m  m 

location  on  the  main  rotor,  and  p  is  the  main  rotor  blade’s 

n 

mass  density. 

The  term  in  Eq  (3.23)  must  now  be  expressed  in  the  body 
axis  system.  First,  Eq  (3.23)  needs  to  be  oriented  in  the 
rotating  hub-fixed  reference  frame  using  the  matrices  found 
in  Eq  (2.4).  In  this  situation,  L.  ,  remains  the  same ,  and 

D  1 

L,.  can  be  redefined  as  L,  The  moment  of  inertia  for  an 
lh  lm 

isolated  blade  can  now  be  expressed  in  the  hub-fixed  rotating 
system: 


IvL  =  hmlhlblvl 


L.  ,L. 
bl  li 


(3.25) 


where 


ml 


■  h-f 


and 


‘lb 


■  M 


Finally,  Eq  (3.25)  can  be  expressed  in  the  body-fixed 
reference  frame: 


(3.20) 


Similarly,  the  expression  for  the  tail  rotor’s  moment  of 
inertia,  in  terms  of  the  body  axis  system,  may  be  derived: 


[^Jb  =  ¥t  [jtb] 


B 


+  *t"t 


(t*  +  t2) 

y  * 

-t  t 
x  y 

-t  t 
X  z 


-t  t 

X  V 

(t2+t*) 
X  z 

-t  t 

y  * 


-t  t 
X  z 

-t  t 

2y  f 

(t  +t  ) 
x  y 


(3.27) 


JB 


where  m^.  is  the  mass  of  a  single  blade  on  the  tail  rotor, 
is  the  number  of  blades  on  the  tail  rotor,  and 
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(3.26) 


-  M  ^  V  ^  V  ^  V  ¥  — -V 


[!te] 


B 


=  L  L 
Btt 


Mote  that  L_.  =  L,  and 
It  In 


lLlb[ITEj 

[jte] 


L  L  L 
blit  tB 


i 8  identical  to  Eq 


(3.23),  except  that  tail  rotor  terms  are  used.  Likewise,  L 


mirrors  L  .  However,  since  the  tail  rotor  is  usually  a 
n 


flapping  rotor  with  low  disk  loading  (2:264),  is  the 

identity  matrix  and  L . .  =  L... 

tD  ID 

Once 


^lcj  is  known,  Eq  (3.10)  can  be  written  as 


*  WB  ?  *  -  *  F'lM 


B 


B 


iBl-  JB 


-«B 


LrJB 


(3.29) 


The  inertial  rate  of  change  of  the  main  rotor's  angular 
momentum  must  also  be  found,  noting  again  that  the  coupling 
of  each  blade's  flap  and  lag  motions  to  the  airframe  are 
ignored  for  the  present.  Assuming  the  rotor’s  angular 


velocity,  Cl  ,  is  constant,  the  following  equation  results: 

ID 


B 


where 


v).  ■ 


(3.30) 


In  terms  of  the  blade 


axis  system, 


rb  = 


0 

0 


-B/b 

co 


0 

0 

o 


(3.31) 


Consequently 


•  (vL  ■ 


h  i  n  a 

ED  ID  ID  D 


Still,  a  rotation 


from  the  blade  axes  to  the  hub-fixed  rotating  axis  system  by 
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means  of  the  flapping  angle 


n 

i ' 


and  the  lagging  angle. 


r 

s 


is  required,  yielding 


im 


m 


'  [“B/I]BLB»L.iLib 

The  rate  of  change  of  the  tail  rotor's  angular  momentum 
about  the  mass  center  may  be  derived  in  a  similar  fashion, 
resulting  in  the  following: 


0 

0 

H  I  O 

L  in  wi  mJ 


(3.32) 


0 

0 

.Vt°t 


(3.33) 


Jb 


The  external  moments  acting  on  the  helicopter  result 
from  Bj  'a  aerodynamic  forces,  as  well  as  from  the  stain  and 
tail  rotors.  For  instance,  a  moment.  Mg,  is  due  to  the 
influence  of  .  In  addition,  the  horizontal  tail's 
aerodynamic  forces  exert  moments  on  the  helicopter. 
Expressed  in  the  body  axis  system,  the  moment  arm  from  the 
helicopter's  center  of  mass  to  the  horizontal  tail  is 


(3.34) 


Hence , 

the  helicopter's 

external  moment 

due 

to  the 

aerodynamic 

forces  acting  on 

as  well  as 

the 

moments 

resulting  from  the  horizontal  tail's  aerodynamic  forces  is 


(3.35) 


The  main  rotor  also  exerta  moments  on  the  helicopter. 
In  addition  to  the  main  rotor's  pitching  moment  about  the 


hub,  M  ,  the  aerodynamic  forces  acting  on  the  main  rotor 
ym 


also  generate  external  moments  on  the  helicopter.  If  the 
moment  arm  from  the  helicopter’s  c.m.  to  the  main  rotor  is 


w,  ■  - 


=  -t  A  ,  then 
z  B 


M, 


0 

II 

♦ 

ym 

0  m 

B 

0/0 

z 


0  0  0 


(3.36) 


In  a  similar  manner,  the  tail  rotor  exerts  a  pitching 
moment  on  the  helicopter.  Further,  if  the  moment  arm  from 
the  c.m.  to  the  tail  rotor  is 


w.  ■  [  5 ' 

-  7  - 


(3.37) 


z  J  B 


the  tail  rotor’s  contribution  to  the  external  moments  acting 
on  the  helicopter  is  as  follows: 


Wa  *  °u 


0 

0 

•f 

M 

L  yJ 

B 

0  t  -t 

z  y 

tot 

z  x 

t  -t  0 

y  x 


(3.38) 


Consolidating  the  various  terms  derived  from  Eq  (3.12) 
results  in  the  helicopter's  moment  equation: 


Vm°Jb 


+LBtLtb  ° 


MtWb 


vrw  v 


irw  irw xrw  vwtui m i.-%rv"%  "v  i  -w  i n  v1  v  n  "S” 


l.’’ 


*  [I=]B‘[-[ic]B^  *K.,o}b-{‘W}b*{“tS}b 

r  „  J 


(3.39) 


SuBimar izing ,  the  equations  of  notion,  as  expressed  in 


the  body  axis  system,  include  the  kinematic  relationships 


found  in  Eq  (3.11),  as  well  as  the  force  equations,  Eq 


(3.18),  and  the  moment  equations,  Eq  (3.39) 


Aerodynamic  Forces  and  Moments 


The  aerodynamic  forces  acting  on  include  the  body's 


side  force,  lift,  and  drag,  as  well  as  the  horizontal  tail's 


lift  and  drag.  For  this  analysis,  the  side  force  term  will 


be  neglected.  If  ,  Vq ,  and  SeQ  are  the  equilibrium  values 


for  the  angle  of  attack,  velocity,  and  elevator  angle,  the 


total  lift  acting  on  the  vehicle  can  be  expressed  in 


equilibrium  and  perturbation  terms  as 


L  =  L  +  ^  Aoc  +  H  AV  +  A6e 

e  do  oV  doe 

e  e  e 


=  L@  +  L^Ac*  +  LyAV  +  L^Ade 


(3.40) 


where  A a,  AV,  and  A6e  are  the  perturbations  from  equilibrium. 


The  lift  acting  on  the  body  and  horizontal  tail  may  be 


written  directly  as  (50:23-29) 


lb  =  5ot2sb°l  «  •  lt  *  Xst(°l  “t  *  V5*)  <3-4,) 

a  v  c*r  <5e  J 


where:  and  are  the  body  and  tail  frontrl  area. 


respectively;  C  .  Cf  ,  and  C  are  the  helicopter,  tail, 
L«  \  L<Se 

and  elevator  lift  curve  slopes;  and  a^,  ®  ot  +  i  +  e  ,  with  i 
being  the  tail  incidence  and  e  the  downwash.  Noting  that  L 
=  +  Lt  ,  then  substituting  Eq  (3.41)  into  Eq  (3.40) 

results  in 


X{“„  (SBCL  *  STCL  ]  ’  St[°L  M 


♦  C  6e 

L6e  °- 


rTc  (SB°L  *  STCL 

'•a  a_ 


j2  fs  C 
o  (_  B 


(3.42) 


pVofo  SBCL;  ST°L  J  *  ST  °L 

~pv2 SC 

O  1  Jj  . 


K  H 


+  C.  6e 
L<5e 


The  drag  term  may  be  expanded  in  a  manner  similar  to  Eq 
(3.40).  Johnson  neglected  drag  due  to  the  horizontal  tail 
(51:208-209)  and  considered  only  the  aerodynamic  drag  acting 
on  the  body.  Expanding  the  drag  coefficient  (53 : 165 ; 54 : 204) 
yields  the  following  term  for  the  drag: 


D  =  ±PV*(0.00«B  ♦  f^SBC*  o2] 


(3.43) 


where  is  the  helicopter’s  wetted  body  area. 
Thus , 


-  rVo  (°  ■  0044B  *  5n®B°L  “o) 


5npVoSBCL  % 

OL 


(3.44) 


5$ 


o  T3  fin  B  La  oj 


The  moment  equation  is  found  to  be  (50 : 28 ; 5 1 : 208 ; 53 : 52) 


If  =  |pV2S  Rcfc  ♦  C  a  +  C  «5e^ 
2  B  ^  mo  “a  m6e  ^ 


(3.45) 


where  c  is  the  distance  from  the  helicopter’s  center  of 
gravity  to  the  moment  reference  center  of  the  airframe,  C 

m 

o 

is  the  helicopter’s  pitching  moment  coefficient  at  zero  lift, 

C  is  the  helicopter’s  pitching  moment  slope  and  results 
m 

a 

from  the  fact  that  B^’s  lift  and  drag  do  not  act  through  the 

helicopter’s  center  of  gravity,  and  C  is  the  pitching 

m6e 

moment  coefficient  due  to  the  elevator. 

Employing  an  equation  similar  to  Eq  (3.40)  results  in 


ipV2SDcfc  +  C  a  +  C  6e 
2  o  B  I  m  mo  m  r  o, 
v  o  a  6e  J 


— pV" S_cC 
2  o  B  m. 


pV  Snc  [c  +  C  a  +  C  <5e  I 
oBlm  mo  oj 

o  a  6e  J 


(3.40) 


?=VoSB5V 


Trim  and  Perturbation  Equations 


Mow  that  the  generalized  equations  of  motion  for  the 
coupled  rotor/ fuselage  system  have  been  derived,  they  may  be 


specialized  for  particular  problems  of  interest,  such  as 


flight  in  the  horizontal  or  vertical  plane.  For  this 
analysis,  straight  and  level  flight  will  be  considered.  This 
requires  y  =  <p  =  p  =  0 .  Further,  assuming  a  flat, 

non-rotating  earth,  since  the  heading  angle  y  is  a  constant, 
it  can  be  set  to  zero.  Parasite  drag  on  the  main  rotor  and 
side  forces  on  the  fuselage  are  neglected.  It  is  assumed 
that  the  sideslip  angle  ft,  as  well  as  the  tail  incidence,  i, 
are  zero.  The  downwash  £  is  a  small  angle  and  can  be 
neglected  (50:21).  The  acceleration  of  gravity  is  constant 
and  the  properties  of  the  atmosphere  are  known  functions  of 
altitude . 

The  equilibruim  and  perturbed  values  of  the  airframe's 
motion  can  be  expressed  as: 

V  =  V  ♦  AV 
o 

a  =  a  +  Aa  (3.47) 

o 

e  =  &  + 

o 

The  equilibruim  and  perturbed  expressions  for  both  the 
tail  and  the  main  rotors’  flapping  and  lagging  motions  were 
summarized  in  Eq  (2.44). 

Based  on  the  above  assumptions,  Eqs  (3.11),  (3.18),  and 

(3.39)  yield  the  following  equilibrium  and  perturbation 
equations,  neglecting  second  order  and  higher  perturbations: 
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FORCE  EQUATIONS 


EQUILIBRIUM 


L  aa  -D  c a  -H  -H. -Mgs$  =  0 
e  o  e  o  m  t  o 

e  e 


-L  oa  -D  sc*  -T  +Y.+Mgc$  =  0 
e  o  e  o  m  t  o 

e  e 


PERTURBATIONS 


AVca  -A aV  aa  =  — 
o  o  o  M 


ao*fi[i':<{(1'e-DJ«=‘-o*(VDe)*a„-VBtJ 


-gA&cd  -AqV  aa 
o  o  o 
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ca  -T  + Y. 
°  ma  I 


-AV(L_ca  +D_aa  +T  -Y. 

1  V  O  V  O  tv 


J-Tn<^.C)+Yt(/?.C> 


-gA &B&  +AqV  ca 
o  o  o 


MOMENT  EQUATIONS 


EQUILIBRIUM 


M  -L_  \t  aa  +J“  ca  1+f  H  +t  H.+t  Y.+M  =  0 
e  TeL  *  °  x  oj  *  «e  z  te  x  te  ym 


(3.48) 


(3.40) 


(3.50) 


(3.51) 


(3.52) 


PERTURBATIONS 
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(3.53) 


irtiere 


ipVZS_(cT  a  +  C.  <5e  \ 
2  o  T1  L  o  L,  of 
^  a^,  6e  J 
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(3.54) 


pv  S-^C-  a  +  C,  6e 
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KINEMATICS 


A$  =  Aq 


(3.55) 


Calculating  the  equilibrium  conditions  requires  an 
iterative  procedure.  Initially,  the  tail  rotor  is  neglected 
and  hover  is  assumed  in  order  to  calculate  the  maximum  value 
of  the  main  rotor's  thrust  with  Eq  (3.49): 


CT  .  Mg/jpnnV] 


(3.56) 


Using  this  value  for  the  thrust  the  main  rotor's  trim 

conditions  can  be  found  from  Eqs  (2.67-2.71).  Terms  for 

a  ,  and  6e  can  then  be  calculated: 
o  o 


V*  'V*  -fV  *  *, 


(3.57) 


D.  *  ^Vo(0  004Ab)  <3  58> 

Mote  that  although  Eq  (3.57)  possesses  tail  rotor 
components,  they  are  neglected  for  now,  as  are  higher  order 
terms.  Using  the  initial  set  of  equilibrium  conditions,  a 
new  value  for  the  main  rotor’s  thrust  coefficient  may  be 
calculated  for  non-hover  conditions  using  Eq  (3.49),  again 
neglecting  the  tail  rotor  components: 

CT  '  )/(PnnX]  <3  M> 

m  v  ey  v  J 


With  this  new  value  for  the  thrust,  equilibrium 

conditions  for  a  ,  9  ,  and  <5e  are  again  calculated,  and  this 

o  o  o 

process  is  repeated  until  CT  converges. 

Once  the  main  rotor's  thrust  coefficient  has  converged, 
the  tail  rotor  is  added  to  the  system,  and  new  equilibrium 
values  must  be  calculated.  To  accomplish  this  task,  the  tail 
rotor’s  thrust  is  calculated  using  Eq  (3.60): 
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Mew  values  for  C_  ,  a  ,  $  ,  and  £e  are  now  calculated 

Too  o 

m 

after  adding  the  tail  rotor  contributions  into  Eqs  (3.57)  and 
(3.59) .  and  the  process  is  iterated  until  both  C  and  C 

T  T. 

m  t 

converge.  These  thrust  coefficients  can  then  he  used  to 

generate  the  trim  values  for  both  blades,  as  well  as  values 

for  &  ,  a  ,  and  <5e  .  The  result  is  a  set  of  equilibrium 
o  o  o 

conditions  for  the  helicopter. 

With  equilibrium  established,  the  perturbation  equations 
can  be  formed.  The  state  vector  includes  the  flap  and  lag 
angles  and  rates  of  each  rotor  and  the  vehicle’s  forward 
velocity,  angle  of  attack,  pitch  angle  and  pitch  rate,  or 


x  (\f)  =  A/3 


u.s)  ,ACu.<0 


.A/3 


U.j) 


Av.A«2.A*.Aq 


(3.61) 


where  Act  =  V  Act,  i  is  either  the  main  or  tail  rotor,  and  j 
/.  o 

varies  from  1  through  M,  the  number  of  blades  on  each  rotor. 

It  is  desirable  to  differentiate  the  system  with  respect  to 

the  main  rotor’s  azimuth,  v  ,  and  nondimensionalize  the 

m 

helicopter’s  speed  with  respect  to  the  main  rotor.  Hence,  Eq 
(3.61)  may  be  written  as 


T 

^ .AS.*.,.** VJ  (3.62) 


where  A/j  =  AV/ (O  R)  ,  /u  =  V/(f)R)  ,  and  Act.  =  fj  Aa. 

BID  O  O  BA  ID  X  O 
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Before  the  equations  of  motion  can  be  put  into  the  form 


x  =  A  (y')x(v')  [Eq  (2.82)]  two  additional  terms  must  be 


accounted  for.  These  terms  result  from  the  tilt  of  the  main 


rotor’s  thrust  vector  with  the  tip-path  plane  (the  plane 


described  by  the  blade  tips  as  they  rotate) ,  as  well  as  the 


tilt  due  to  the  helicopter’s  angular  velocity  (2:164; 


788-789).  These  terms,  added  to  Ajj  and  A&  ,  are  generated 


by  the  following  stability  derivative: 
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Mote  that  Eqs  (3.50),  (3.51),  and  (3.53)  possess  two 


sets  of  terms,  one  including  the  vehicle’s  state,  x^y*). 


while  the  other  considers  the  coupling  terms  resulting  from 


the  interaction  of  the  blade  dynamics,  x^(v').  Bearranging 


the  equations  results  in  the  following: 
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Because  the  tail  rotor  is  a  flapping  rotor,  lag  angle 


and  rate  perturbations  can  be  neglected.  Consequently,  the 
main  and  tail  rotor  terms  in  the  D  matrix  differ.  For  the 
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(3.68) 


(3.69) 
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where,  with  i  denoting  the  summation,  from  one  to  H,  of  the 
number  of  blades  on  either  the  main  or  the  tail  rotor. 
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(3.71) 


The  flap-lag  equations  of  motion  previously  derived  for 

an  isolated  blade  still  hold  for  the  coupled  rotor/fuselage 

system,  with  a  few  exceptions.  Adding  the  fuselage 

introduces  additional  coupling  terms  resulting  from  A/j,  Aa^  , 
# 

Ad,  and  Ad  .  These  coupling  terms  must  be  determined  and 
added  to  the  existing  blade  equations.  In  addition,  recall 
that  the  moment  equations  for  were  derived  ignoring  any 
contributions  from  each  blade’s  flap  and  lag  motions. 
Coupling  the  rotors  to  the  airframe  introduces  blade  flap  and 
lag  terms  in  the  airframe's  moment  equations,  and  these  must 
be  accounted  for.  Finally,  in  an  attempt  to  expand  the 


helicopter’s  operating  envelope  reversed  flow  should  be 


considered.  Each  of  these  issues  Bust  be  addressed 
separately.  It  should  be  noted  that  the  following  notation 
is  consistent  with  the  notation  used  in  deriving  the 
equations  for  an  isolated  blade’s  flap-lag  motion. 

One  new  term  introduced  by  coupling  the  fuselage  to  the 
rotor  is  the  pitch  effect  on  the  blade’s  rotation;  i.e., 
there  is  a  ‘rocking*  motion  of  the  blade  due  to  the  body’s 
pitching  motion.  This  can  be  added  to  the  angular  velocity 
of  the  non-rotating  hub-fixed  system  (considered  to  be 
inertial  in  the  previous  chapter) .  The  resulting  angular 
velocity  becomes 

wb/I  =  QAh  -  +  C^b+  (3.72) 


Referring  to  Eq  (2.11),  the  absolute  velocity  of  a  mass 
point  on  the  blade,  expressed  in  the  hub-fixed  rotating 
system  is  then 


-r(Cc/?sC  +  /?s/?cC  +  OsC  -  #s/?cC) 
r(CcC  ♦  Ocr?cC) 
r(/?c/3cC  -  Cs/?sC  “  $c/?cO 


(3.73) 


The  Lagrangian  becomes,  for  the  main  rotor. 


L  =  |l  |^C2+/?2c2C+&2c2C*fi2  (s2C+cVc2c}-2/7£c2C 

+  2nCc^+2fV)cCsC8^-20es^cCscj+^Ic^2+|MV2-ik^2-ikcC2  (3.74) 

where  the  body  terms  are  now  included.  Vote  that  the 
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translational  kinetic  energy  of  the  airframe  is  included  in 
the  Lagrangian.  Because  ft,  and  C  do  not  appear  in  this 
term  it  will  not  contribute  any  components  to  Lagrange’s 


equations  of  motion. 

Assuming  the  angles  ft,  and  C  are  small  and  neglecting 
terms  higher  than  second  order,  Lagrange’s  equations  of 
motion  are: 


*  -  ®  *  2njft  ♦  *  V'v  ft  * 

V1] 

ft  =  V1* 

(3.75) 

C  -  20  ft  ft  ♦  Cl  &  ft  +  c  C/I  + 

m  m  c,  in 

kCC/Im  =  V*. 

(3.76) 

S  -  \[ft  -  -  V'?]/(I.*Ic) 

-  vl 

[1m+Ic] 

(3.77) 

The  trim  equations  will  not  change.  Differentiating  the 

above  equations  by  the  main  rotor’s  azimuth  angle,  ,  yields 

l& 

the  following  perturbation  equations  of  motion: 
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(3.78) 

(3.79) 

(3.80) 


Note  that,  by  including  the  pitching  terms,  the 
perturbed  equations  of  motion  are  indeed  altered.  In  order 
to  generate  the  perturbed  equations  in  the  form  of  Eq  (2.82), 
Eqs  (3.78)  and  (3.80)  must  be  simplified: 
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Eq  (3.82)  is  the  moment  equation  for  B^  and  included  the 
contributiona  of  each  blade 'a  flap  and  lag  motions  on  the 
airframe.  Thia  equation  replace?  Eq  (3.53)  when  coupled 
rotor/f uaelage  dynamics  are  considered.  However,  when  the 
uncoupled  dynamics  are  to  be  examined,  Eq  (3.53)  ia  used  as 
the  moment  equation  for  Bj. 

The  coupling  terms  associated  with  the  aerodynamic 
forces  must  also  be  calculated.  Referring  to  Figure  6, 
changes  in  the  helicopter’s  velocity  and  angle  of  attack  do 
indeed  influence  the  aerodynamic  forces  acting  on  the 
helicopter  blade  through  the  advance  ratio,  ij,  and  inflow 
ratio,  X.  These  terms  can  be  expanded  to  be: 

M  =  v  ♦  A^eosa  -  Aa.sina  (3.83) 

o  o  1  o 

X  =  X  +  A^sina  +  Aa.cosa  (3.84) 

o  o  1  o 

It  will  be  assumed  that  changes  in  the  induced  velocity 
due  to  small  changes  in  V  and  o  can  be  neglected.  Hence, 
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Au  =  AXsinot  +  A/jcosot 
o  o 


(3.85) 


Aa,  =  AAcosa  -  Ausina 
1  o  o 


(3.86) 


If  CQ  =  0,  then  perturbations  of  the  relative  normal  and 
tangential  velocity  components  can  be  expressed  as  (16:880) 


Jp  =  RQ  j^cosv'A^i  -  Axj 


AU  =  BDsinv'A^ 

T 


(3.87) 


(3.88) 


Using  the  perturbed  values  of  the  elemental  aerodynamic 
forces,  A(dF  )  and  A(dF  )  ,  found  in  Eqs  (2.76)  and  (2.77)  , 

and  integrating  over  the  length  of  the  blade,  expressions  for 

•  * 

AQ^  and  AQ^.  may  be  determined: 
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where  t_  and  t_  are  given  by  Eq  (3.71). 
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The  expression  for  AQ^,  previously  defined  in  Eq  (3.53), 
is  repeated  below  to  show  only  the  contributions  from  A/j  and 
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Combining  the  results  of  Eqs  (2 . 84) - (2 . 86)  ,  (3.70), 

(3.81),  (3.82),  and  (3 . 89) - (3 . 01)  yields  a  matrix  in  the  form 

__  0  ____ 

x(y)  =  A(v,)x(v')  .  and  includes  not  only  the  perturbations  of 


A/7  and  AC  but  also  the  coupling  terms  resulting  from  A/u ,  Aa^  , 


AO,  and  AO  .  For  a  single  blade,  this  matrix  is  as  follows: 


0  0  0 


XjtV')  = 


0  0  0  0  -  .  0  000-,. 

x.  (v>)  +  _  x_  (y/) 

C31  °32  C33  34  1  ®31  ®32  0  ®34  2 


(3.02) 


°41  C42  C43  C44 


e  e  0  e 

41  42  44 


where  (referring  to  Eqs  (2.85)  and  (3.67)  as  required) 
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If  the  helicopter's  operating  envelope  is  to  be 
expanded,  advance  ratios  greater  than  0.50  must  be 
considered.  This  means  reversed  flow  should  be  accounted 
for.  As  mas  previously  discussed,  Sissingh  (8)  described 
three  different  flow  regions  a  blade  encounters.  In  the 
first  region,  called  normal  flow,  the  air  approaches  the 
blade  from  the  leading  edge,  while  in  reversed  flow  the  air 
approaches  from  the  trailing  edge.  In  the  mixed  flow  region 
part  of  the  blade  experiences  normal  flow  while  the  other 
part  encounters  reversed  flow.  The  aerodynamic  flapping  and 
lagging  moments  describing  normal  flow  were  already  defined 
in  Eqs  (2.37)  and  (2.38);  the  moments  for  reversed  flow  are 
simply  the  negatives  of  the  normal  flow  moments.  All  that 
remains  is  to  derive  the  aerodynamic  moments  for  mixed  flow. 

According  to  Sissingh,  the  boundaries  of  the  mixed  flow 
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region  occur  from  y  =  n  to  yt  =  rt  ♦  c  and  from  y>  -  2 n  -  c  to  y> 

-t  ' 

=  2n ,  where  c  =  sin  (1/#j)  (8:57).  Hence,  the  moments  Q„ 


and  Q_.  are  modified  to  be 
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If  the  above  equations  are  solved,  the  equations  of 
motion  for  the  blade  elements  in  the  mixed  flow  region  will 
be  modified  from  Eq  (3.93)  to  be 
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Note  that  the  D  matrix  terms  for  the  tail  rotor  will  be 
modified  in  a  similar  manner. 

To  simplify  the  helicopter’s  equations  of  motion,  the 
tail  rotor  can  be  neglected.  Then  the  equations  of  motion  of 
an  N-bladed  coupled  rotor/fuselage  helicopter  may  be 
expressed  in  the  form  of  Eq  (2.82),  where  x(yO  is  defined  in 
Eq  (3.62),  and 
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where  the  f  terms  also  account  for  reversed  flow  [refer  to 
Eq  (3.95)]. 

Note  that  the  terms  in  A(v/)  have  been  previously 
defined,  and  the  dimensions  of  this  matrix  depend  upon  the 
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number  of  blades  in  the  main  rotor.  For  instance,  if  a 
four-bladed  main  rotor  is  specified,  then  the  state  matrix, 
x(v')  ,  will  be  a  20-column  vector,  where  the  first  J6 
components  are  the  flap  and  lag  angles  and  rates  of  each  of 
the  four  blades.  The  last  four  components  will  consist  of 
the  helicopter's  airframe  terms.  The  coupled  rotor/fuselage 
equations  reduce  to  the  equations  for  the  blades  alone  and 
the  helicopter’s  body  (excluding  the  rotor)  alone  when  the 
coupling  terms  are  ignored.  The  coupling  terms  are  included 
to  account  for  the  helicopter  body  pitch,  velocity,  and  angle 
of  attack  changes  being  considered  in  the  blade’s  equations 
of  motion,  as  well  as  the  rotor  blades’  flapping  and  lagging 
changes  being  included  in  the  helicopter  body’s  equations. 

Ver i f ication  of  the  Helicopter ’ s  Equations  of  Motion 

The  equations  of  motion  have  been  derived  for  a  coupled 
rotor/f uselage  system;  however,  they  still  need  to  be 
verified.  Other  studies  have  derived  and  examined  the 
equations  of  motion  for  a  coupled  helicopter,  but,  in  many 
cases,  it  is  not  feasible  to  make  a  direct  correlation 
between  those  studies  and  the  present  effort.  For  instance. 
Bo us man  (28)  investigated  the  aeromechanical  stability  of  a 
helicopter  on  the  ground  and  in  hover.  Friedmann  and 
Venkatesan  (33; 39; 40)  and  Straub  and  Warmbrodt  (31)  presented 
analytical  models  to  examine  the  aeromechanical  stability  of 


a  helicopter  in  ground  resonance.  Butko wski  (32)  and  Hodges 
(35;38)  modeled  their  systems  using  finite  element  theory. 
On  the  other  hand,  Johnson  provided  an  analytical 
examinination  of  a  helicopter's  motion  both  in  hover  and  in 
forward  flight  (2:774-843).  His  analysis  can  be  directly 
correlated  with  the  equations  derived  earlier  in  this 
chapter,  and  therefore  warrents  closer  examination. 

Because  a  helicopter's  flying  qualities  are  different  in 
hover  and  in  forward  flight,  Johnson  analyzes  the  two  regimes 
separately  (2:775).  For  hover,  Johnson  assumes  the  aircraft 
has  complete  axisymmetry  and  separates  the  vertical  and 
longitudinal-lateral  dynamics  to  simplify  his  analysis.  In 
addition,  he  assumes  that  the  longitudinal  and  lateral 
dynamics  are  also  separable  and  only  the  low  frequency 
dynamics  of  the  rotor  are  used,  since  the  motion  involved  in 
helicopter  flight  dynamics  is  slow  when  compared  to  the 
rotor’s  motion  (2:775;779).  For  the  case  of  longitudinal 
dynamics,  two  degrees  of  freedom  exist  in  hover:  pitch;  and 
longitudinal  velocity.  For  an  articulated  rotor  with  no  flap 
hinge  offset  and  no  pitch-flap  coupling,  the  helicopter’s 
pitching  moment  is  due  only  to  the  in-plane  hub  force. 
Consequently,  the  system’s  characteristic  equation  generates 
three  roots:  a  negative  real  root  due  principally  to  the 
main  rotor’s  pitch  damping;  and  a  long  period,  mildly 
unstable  oscillation  resulting  from  the  coupling  of  the  pitch 
and  longitudinal  velocity  (2:787-793). 

Johnson  also  examines  a  helicopter’s  performance  in 


forward  flight,  which  is  significantly  different  from  the 


hover  case.  As  he  states. 

Forward  speed  introduces  new  forces  acting  on  the 
helicopter:  centrifugal  forces  due  to  the  rotation  of 
the  trim  velocity  vector  by  the  angular  velocity  of  the 
body  axes;  aerodynamic  forces  on  the  fuselage  and  tail; 
and  major  rotor  forces  that  are  proportional  to  the 
advance  ratio  [2:822]. 

Again,  Johnson  assumes  the  longitudinal  and  lateral 
dynamics  can  be  analyzed  Separately.  In  forward  flight, 
three  longitudinal  degrees  of  freedom  exist:  longitudinal 
velocity;  pitch  attitude;  and  vertical  velocity.  As  in 
hover,  Johnson  obtains  the  rotor  forces  and  moments  acting  on 
the  helicopter  from  the  low  frequency  response;  the  rotor 
dynamics  do  not  add  degrees  of  freedom  to  the  system  (2:824). 
There  are  three  primary  influences  the  helicopter’s 
longitudinal  dynamics:  the  pitching  moment  due  to  vertical 


velocity;  the  vertical  acceleration  due  to  pitch  rate;  and 
the  helicopter’s  longitudinal  inertia  (2:829).  In  hover  the 
system’s  characteristic  equation  typically  yields  two  real. 


negative  roots  for  the  vertical  and  pitch  moments  and  a 


complex  conjugate  pair  in  the  right-half  plane  resulting  from 


the  longitudinal  velocity  and  pitch  coupling.  Without  a 


horizontal  tail,  the  main  rotor  produces  a  net  angle  of 


attack  instability.  As  the  velocity  increases,  the  vertical 


mode  becomes  more  unstable  while  the  pitch  mode  becomes  more 


stable.  For  the  oscillatory  mode,  the  unstable  period 


increases  and  the  damping  decreases.  However,  with  a  large 


enough  horizontal  tail,  the  helicopter,  in  forward  flight. 
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TABLE  I 

RSRA  Parameters 


) 

82.25 

21.52 

om  eg 

to  moment  center. 

m)  0.22 

m2) 

135.80 

rtia . 

(kg-m2) 

08102.41 

TAILS: 

HORIZONTAL 

VERTICAL 

area  (m  ) 

3.29 

10.584 

x*  Cm) 

13.005 

11.624 

a 

y  Cm) 

0.0 

o 

o 

z*  Cm) 

3.512 

1.758 

"distance  from  eg  to  tall  moment 

center 

BOTORS : 

MAIN 

TAIL 

number  of  blades 

4 

4 

radius  (m) 

9.450 

1.616 

chord  Cm) 

0.648 

0.306 

weight  of  a  single  blade  (kg) 

123.391 

9.954 

angular  velocity  (rad/sec) 

22.579 

132.068 

2 

mass  moment  of  inertia  (kg-m  ) 

3861.722 

4.585 

solidity  (cr) 

0.0873 

0.241 

Lock  number  (y) 

10.300 

3.500 

x*  (m) 

0.0 

11.252 

ft 

y  (m) 

0.0 

0.670 

z*  (m) 

«...  .  .  . 

2.108 

1.385 

Table  II 


RSRA  Aerodynamic  Coefficients 


LIFT 

MOMENT 

CL 

aB 

0.13751 

C 

m 

o 

0.00000 

CL 

°T 

0.27510 

C 

at 

-0.85944 

Se 

0 . 05813 

C 

“6e 

-0.04125 

data  contained  in  the  above  mentioned  souces.  The  main  rotor 
has  four  blades  and  is  articulated,  without  any  hinge  offset; 
hence,  there  is  no  hub  moment  transmitted  to  the  helicopter 
(2: 150 ; 788) . 

Using  this  helicopter  model,  the  vehicle’s  performance 
is  evaluated  for  various  flight  conditions.  To  verify  the 
equations  of  motion,  the  blade  equations  are  first  uncoupled 
from  the  body  and  the  results  compared  to  those  obtained  with 
the  isolated  blade  flap-lag  equations.  As  anticipated,  the 
results  are  identical.  In  addition,  for  the  four-bladed 
helicopter,  the  results  mirror  those  reported  by  Calico  and 
Wiesel  when  they  examined  a  two-bladed  flapping  rotor 

(43:62).  That  is,  when  the  body  and  rotor  are  uncoupled,  the 

✓ 

four  pairs  of  Poincare  exponents  associated  with  each  of  the 

blade’s  flapping  modes  are  identical.  Similarly,  the 

✓ 

Poincare  exponents  associated  with  the  lagging  modes  are 
identical.  Figures  11(a)  and  12(a)  show  the  real  parts  of 
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the  uncoupled  flap  and  lag  nodes  as  the  velocity  is  increased 
from  hover  through  300  m/s  (/u  =  1.061).  For  these  cases,  the 
blade  natural  frequencies  are  set  at  p  =  1.00  and  =  0.50, 

and  =  0.000. 

Both  figures  show  that  Initially  the  flap  and  lag  modes 

can  be  viewed  as  four  identical  pairs  of  complex  conjugate 

✓ 

Poincare  exponents.  The  flap  modes  split  off  as  two  sets  of 

✓ 

roots,  each  of  which  consists  of  four  identical  real  Poincare 
exponents.  Note  that  Figure  11(a)  shows  one  of  these  sets  of 
roots  becoming  unstable  at  very  high  speeds.  At  low  speeds, 

the  lag  modes  are  also  four  identical  pairs  of  complex 

✓ 

conjugate  Poincare  exponents.  Once  again,  the  modes  become 
real  as  the  velocity  increases,  but,  according  to  Figure 
12(a),  the  modes  become  increasingly  stable  oscillations  at 
high  speeds. 

The  uncoupled  flap  and  lag  modes  can  also  be  examined 
when  reversed  flow  is  considered.  Examination  of  Figures 
11(a)  and  12(a)  shows  that  reversed  flow  does  not 
significantly  alter  the  uncoupled  flap  and  lag  modes. 
However,  reversed  flow  does  stabilize  the  flap  siodes  at  the 

upper  end  of  the  velocity  spectrum.  Further,  the  range  of 

* 

velocities  where  the  lag  oxides  possess  real  Poincare 
exponents  (rather  than  complex  conjugate  pairs)  increases 
when  reversed  flow  is  added  to  the  system. 

Now  that  the  uncoupled  blade  flap  and  lag  siotions  are 
known,  the  uncoupled  motion  of  the  rigid  airframe  needs  to  be 
verified.  To  this  end,  the  airframe’s  roots  must  be 


Figure  12.  Uncontrolled  Lag  Modes  as  a  Function  of  Velocity 


calculated  at  varying  velocities  and  compared  with  Johnson  a 


observations 


Consequently,  the  helicopter's  uncoupled 


longitudinal  dynamics  are  examined  in  hover  and  in  forward 


flight. 


In  hover,  the  airframe’s  characteristic  equation 


generates  four  roots.  Two  real,  but  stable,  Poincare 


exponents  correspond  to  Johnson’s  roots  for  the  helicopter’s 


vertical  and  pitch  motions.  The  other  two  roots  represent  an 


unstable  oscillation  resulting  from  the  coupling  of  the  pitch 


and  longitudinal  velocity.  This  again  mirrors  Johnson’s 


observations  (2 : 793 ; 829-831 ) . 


Forward  flight  can  be  examined  by  looking  at  two 


different  scenarios:  a  helicopter  with  and  without  a 


horizontal  tail.  Johnson's  analysis  of  the  helicopter  in 


forward  flight  is  based  on  several  assumptions  which  simplify 


the  system's  characteristic  equation.  Specifically,  the 


helicopter’s  forward  velocity  stability  derivative  t®2i  from 


Eq  (3.67)]  is  considered  to  be  very  small.  In  addition,  the 


pitch  moment  stability  derivatives  are  directly  proportional 


to  the  corresponding  longitudinal  force  derivatives:  a. 


-MR  l  a.,/1  ;  a  =  -MR  l  a,_/I  ;  a..  =  -MR  l  a.. /I 

mz  11  c  42  m  z  12  c  44  m  z  14  c 


[again  referring  to  Eq  (3.67)]  (2:788;829).  Using  these 


assumptions ,  the  uncoupled  airframe's  characteristic  equation 


mirrors  Johnson’s  equation  (2:829),  and,  consequently,  the 


responses  are  similar.  Figure  13  shows  the  root  locus  of  the 


longitudinal  roots,  with  and  without  a  horizontal  tail,  as 


the  velocity  is  varied.  Without  a  tail,  an  increase  in 
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Figure  13.  Velocity  Hoot  Locus  of  the  Uncontrolled, 
Uncoupled  Body  Modes;  Ho  Reversed  Flow;  Includes  Only 
Johnson’s  Terms  (2: 788; 829) 

forward  velocity  tends  to  destabilise  both  the  helicopter’s 
vertical  motion  and  the  longitudinal  velocity.  On  the  other 
hand,  the  horizontal  tail  chosen  for  the  RSRA  does  cause  the 
vertical  and  pitch  roots  of  hover  to  transform  into  stable 
oscillatory  roots  as  the  velocity  increases.  Further,  the 
oscillatory  roots  representing  the  coupling  between  the 
helicopter’s  pitch  and  longitudinal  velocity  become  more 
stable  with  increasing  speed.  Once  again,  these  results 
mirror  Johnson’s  analysis. 

Including  all  the  airframe  contributions  in  the  system’s 
characteristic  equation  does  not  significantly  alter  the 
velocity  root  locus,  as  shown  in  Figure  14(a).  The  only 


noticeable  change  is  that,  when  the  tail  is  added,  the  modes 


representing  the  coupling  between  the  pitch  end  longitudinal 
velocity  do  not  stabilize,  even  at  very  high  velocities. 
Figure  14(b)  plots  the  real  parts  of  the  uncoupled  airframe 
modes  as  a  function  of  velocity. 

It  should  be  noted  that  adding  reversed  flow  to  the 
system  does  not  significantly  influence  the  airframe's  roots. 
However,  Figures  15(a)  and  15(b)  do  show  that,  at  very  high 
speeds,  adding  reversed  flow  to  the  system  does  stabilize  the 
modes  representing  the  coupling  between  the  pitch  and 
longitudinal  velocity. 

As  the  above  discussions  illustrate,  the  uncoupled 
airframe  does  mirror  Johnson's  analysis  as  the  velocity  is 
varied.  The  final  step,  then,  is  to  add  the  components  which 
couple  the  airframe  and  the  flap  and  lag  terms  of  each  blade 
on  the  main  rotor.  With  p  *  1.00,  *  0.50,  and  = 
0.000,  the  coupled  rotor/fuselage  system’s  roots  will  be 
examined  as  the  velocity  is  increased  from  hover  through  300 
m/s . 

The  velocity  root  locus  of  the  coupled  body  modes. 
Figure  16(a),  differs  significantly  from  the  velocity  root 
locus  of  the  uncoupled  airframe  with  reversed  flow  added  to 
the  system  [Figure  15(a)].  In  addition.  Figure  16(b)  shows 
that  the  real  parts  of  the  coupled  airframe  modes  differ  from 
the  uncoupled  cases  [Figures  14(b)  and  15(b)]  as  the  speed 
increases.  How,  recall  that,  without  any  coupling  terms  in 
the  system,  the  vertical  and  pitch  modes  generated  stable 
oscillations  as  the  forward  velocity  increased.  However, 
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Figure  16.  Dncontrol led ,  Coupled  Body  Modes;  Reversed  Floi 


this  is  not  the  case  when  the  coupling  terns  are  added.  The 
hover  roots  are  nearly  the  same,  but  now,  as  the  velocity 
increases,  the  root  representing  the  helicopter's  vertical 
■lotion  becomes  more  stable  while  the  pitch  root  becomes  less 
stable.  In  addition,  the  unstable  oscillations  in  hover 
become  more  unstable  with  increasing  velocity.  This  trend 
also  differs  from  that  exhibited  by  the  uncoupled  modes 
representing  the  coupling  between  the  body’s  pitch  and 
longitudinal  velocity,  for  those  modes  became  more  stable  as 
the  speed  increased.  Finally,  notice  that,  at  higher 
velocities,  the  pitch  mode  and  one  of  the  modes  representing 
the  coupling  between  pitch  and  longitudinal  velocity  generate 


a  slightly  unstable  oscillatory  pair  of  roots  for  a  bit 


before  becoming  unstable  real  modes  again. 

The  coupled  flap  and  lag  modes  [Figures  11(b)  and  12(b)] 
also  differ  from  the  uncoupled  modes  [shown  in  Figures  11(a) 
and  12(a)]  as  the  velocity  is  varied.  Figure  17  shows  the 
velocity  root  loci  of  the  eight  flap  modes.  The  magnitudes 
of  the  imaginary  parts  of  each  of  the  lag  modes  remain  at 
approximately  0.500  throughout  the  velocity  spectrum 
examined.  At  low  speeds.  Figures  11(b)  and  17  show  that  the 
eight  coupled  flap  modes  are  very  similar  and  mirror  the 
uncoupled  modes.  However,  as  the  speed  increases  the  coupled 
modes  diverge.  In  fact,  only  one  pair  of  modes  is  real  above 
80  m/s.  Further,  at  very  high  speeds,  two  other  modes  become 
unstable.  The  remaining  four  modes  consist  of  two  pairs  of 
stable  oscillatory  roots  through  300  m/s.  As  for  the  lag 
motions  of  the  four  blades .Figure  12(b)  indicates  that  the 
coupled  lag  modes  remain  stable  complex  conjugate  pairs.  Six 
modes  are  oscillatory  roots  similar  to  the  uncontrolled  case 
for  low  speeds.  Two  roots,  however,  differ  drastically  from 
the  other  six  roots  as  the  velocity  is  increased.  While 
stable  throughout  the  velocity  spectrum  analyzed,  these  two 

f 

Poincare  exponents  are  much  more  stable  (except  at  very  high 
speeds)  than  the  other  modes,  even  in  hover. 

Clearly,  coupling  the  flap  and  lag  motions  of  the  main 
rotor’s  four  blades  to  the  helicopter’s  airframe  (including 
the  horizontal  tail)  does  indeed  influence  the  results.  For 
this  analysis,  the  coupled  equations  of  motion  will  be  used 


i  < 


for  the  remainder  of  this  study. 
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IV.  Modal  Control  Theor' 


Once  a  set  of  equations  of  motion  have  been  derived  and 
the  stability  regions  determined,  attention  can  be  focused  on 
developing  a  control  technique  which  reduces  or  eliminates 
the  unstable  regions.  The  modal  control  theory  used  in  this 
study  is  summarized  below. 


Floquet  Theory 


The  stability  of  a  system  of  linear  equations  with 
periodic  coefficients  may  be  determined  using  Floquet  theory, 
which  only  requires  the  knowledge  of  the  state  transition 


matrix  at  the  end  of  one  period. 


Because  Floquet  theory 


involves  no  assumptions  beyond  those  used  in  deriving  the 
equations  of  motion,  the  accuracy  of  its  results  depends  only 
upon  the  computational  procedures  used  (7:26-27). 

Understanding  modal  control  theory  for  periodic  systems 
requires  a  brief  review  of  Floquet  theory.  The  ensuing 
development  follows  closely  that  found  in  Calico  and  Wiesel 
(43 ; 62) . 

Consider  a  set  of  linear  ordinary  differential  equations 
of  the  form 


x(t)  =  A ( t) x (t) 


(4. 1) 


m.:':  : 


where  the  matrix  A(t)  is  periodic  with  period  T. 


For 


convenience,  the  independent  variable  is  changed  to  v*.  and  Eq 
(4.1)  becomes 

x  (vO  =  Aty'Jxfv')  (4.2) 


where  A(y>)  =  A  ( y/  +  2 n)  The  solution  to  these  equations 

can  be  written  in  terms  of  the  state  transition  matrix, 
*(^.0),  as  follows: 


x(v')  =  *(y/,0)x(0) 

where  satisfies  the  matrix  equations 

*  (y.0)  =  A  (y/)  *  (\t/. 0) 

*(0,0)  =  I 


(4.3) 


(4.4) 


A  direct  result  of  Floquet  theory  is  that  $(y.',  0)  can  be 
written  as 

i(V'.O)  =  F(y/)eJv/F_1(0)  (4.5) 

where  the  matrix  F(y/)  is  periodic  with  the  same  period  as 
A(V')  ,  and  J  is  a  constant  matrix  which  can  be  expressed  in 
the  Jordan  normal  form. 

* 

The  diagonal  elements  of  J  are  Poincare  exponents. 

f 

Vote  that  the  system’s  stability  is  governed  by  the  Poincare 
exponents  alone,  since  F(v')  is  periodic  and  therefore 
bounded.  That  is,  if  all  the  real  parts  of  the  X A  are 
negative,  the  system  is  stable.  However,  the  system  is 


100 


1 


unstable  if  any  of  the  Poincare  exponents  have  positive  real 
parts . 

How,  because  the  knowledge  of  the  state  transition 
matrix  over  one  period  determines  the  solution  everywhere, 
solving  for  all  time  thus  requires  solving  the  matrices  J  and 
F(y).  Since  F  is  a  periodic  matrix,  F(0)  =  F(2n)  ,  and 

evaluating  Eq  (4.5)  at  the  end  of  one  period  yields 

§(2rr,0)  =  F(0)e2rTJF'4(0)  =  F(0)AF-1(0)  (4.0) 

where  A  consists  of  the  eigenvalues  of  $(2rr,0). 

Thus,  F(0)  can  be  referred  to  as  the  matrix  of 
eigenvectors  of  the  monodromy  matrix  $(2rr,0).  In  addition, 
the  eigenvalues,  or  characteristic  mul  tipi  iera ,  of  $(2rr,0) 
are  related  to  the  Poincare  exponents  by 


A.  =  exp(2nX,) 


(4.7) 


In  general,  both  A  and  X  are  complex  quantities;  thus 


A.  =  A.  +  (A 

A  *1 


♦  ^4 


(4.8) 

(4.9) 


from  which 


•  kin[(A^  *  vi 


“*  *  ktan"‘<AAI/AAB> 


(4.10) 

(4.11) 


Constructing  a  complete  solution  to  Eq  (4.5)  requires 
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the  knowledge  of  the  eigenvector  matrix  F(v')  over  one  period. 
Substituting  Floquet’s  results  from  Eq  (4.5)  into  Eq  (4.4) 
results  in 

f'  (w)  =  A(^)F(v)  -  F(w)J  (4.12) 

The  initial  conditions  for  this  differential  equation 

are  available  from  the  eigenvector  matrix  found  in  Eq  (4.6). 

Hence,  by  numerically  integrating  Eq  (4.12)  over  one  period, 

F(v')  can  be  found  throughout  the  entire  period. 

Even  though  F(^)  and  J  can  be  calculated,  the  results 

may  be  inconvenient  if  the  matrices  are  complex.  If  this  is 

the  case,  both  matrices  can  be  rearranged  to  stake  them  real. 

F(vO  should  consist  of  column  vectors  T  which  are  either  1) 

the  real-valued  eigenvectors  associated  with  the  real 

Poincare  exponents,  or  2)  two  columns  representing  the  real 

and  imaginary  parts  of  the  eigenvector,  f  and  f.  , 

real  imag 

✓ 

associated  with  a  complex  conjugate  pair  of  Poincare 
exponents.  The  J  matrix  will  consist  of  either  1)  diagonal 

r 

entries  of  the  real  Poincare  exponents,  or  2)  diagonal  blocks 
of  the  form 

‘  <4 
"“A 

* 

for  a  complex  conjugate  pair  of  Poincare  exponents  (62:672). 

How,  it  is  often  necessary  to  find  the  solution  to  the 
inverse  eigenvector  matrix,  F  1(v/).  Bather  than  numerically 
inverting  F(v'),  which  is  expensive  and  leads  to  round-off 


“A 

*4, 
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error  (62:672),  F  (yj)  can  be  calculated  by  differentiating 
FF  1  =  I  and  substituting  Eq  (4.12)  into  the  result, 
yielding 


-F_1(V)A(W)  ♦  JF_1(¥0 


(4. 14) 


This  equation  can  also  be  numerically  integrated,  and 
the  results  reduced  to  a  convenient,  usable  form  by  harmonic 
analysis  (63:108-109). 

By  introducing  a  set  of  modal  variables,  r> ,  such  that 


x(v>)  =  F  ly/)r)(v>) 


(4.15) 


the  periodic  system  of  Eq  (4.2)  can  be  written  as 


7)  (v>)  =  F 


“(VO  [ 


A(vOF(v/>  -  F 


’  (V')  j  7) 


(V>)  =  J7)(v>) 


(4. 16) 


Hence,  the  periodic  system  can  be  reduced  to  a 
constant-coefficient  system  by  using  the  eigenvector  matrix 
F (y)  as  a  periodic  transformation  (62:672). 

Control  of  the  system  can  be  accomplished  by  adding 

/ 

state  variable  feedback  to  change  the  unstable  Poincare 
exponents.  Consider  the  standard  control  problem 


x  ( V')  =  A  (v')x(y/)  +  B(vj)u(v/) 


(4.17) 


where  u(v')  is  the  control  vector  and  B(y'),  a  matrix  which  is 
periodic  with  the  same  period  as  the  fundamental  dynamical 
system,  determines  the  control  distribution. 

Assuming  full  state  feedback. 


e 


u(v)  =  K(v)x(v^)  (4.10) 

where  K(y)  is  the  gain  matrix.  Thus,  Eq  (4.17)  can  be 
written  as 

xc(V)  =  [*(,)  +  B(v')K(v/)  Jxc(V/)  (4.19) 

where  x^  is  the  closed-loop  state. 

By  introducing  the  modal  variables  into  the  feedback 
control  system  of  Eq  (4.19),  the  following  results: 

*  -i 

r)  (y)  =  Jv(y)  *  F  (y)  B(y)u(.y) 

or  (4.20) 

77  c(V/)  =  +  F_1  (v/)B(v^)K(v,)F(V/)J  Vc  (y) 

If  K(v')  is  chosen  to  be  periodic  with  the  same  period  as 
Eq  (4.2),  then  Eq  (4.20)  is  a  Floquet  problem.  To  insure 
that  the  closed- loop  system  meets  the  desired  specifications, 
K(v')  must  be  calculated.  Calico  and  Wiesel  (43; 62)  have 

developed  a  modal  control  technique  which  generates  the 

✓ 

required  K(y')  by  changing  the  unstable  Poincare  exponents  and 
leaving  the  others  unaltered.  Both  scalar  and  vector  control 
are  considered. 

Scalar  Control 

Consider  first  a  system  of  four  equations  of  motion 
where  either  one  pair  of  complex  conjugate  Poincare  exponents 
or  two  real  Poincare  exponents  are  unstable.  This  situation 
occurs  when  an  isolated  helicopter  blade  experiencing 
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flap-lag  motion  encounters  unstable  regions  at  various 
advance  ratios,  as  was  discussed  in  Chapter  II.  The  modal 


variables  may  be  numbered  so  that  T)^  and  T7^  are  the  two 
stable  modes  while  r)_  and  r)  are  the  unstable  modes.  The 

control  can  be  given  by 


u  (vO  =  k  (v')r>(V/) 


(4.21) 


where  k  (¥>)  is  a  row  matrix  of  modal  feedback  gains,  and  is 


given  by 


-T 

k1  (v) 


-  Jkj  <w)  . 


k2(y)  ,  kg(v') 


,  k4(v>)J 


(4.22) 


The  modal  feedback  control  system  in  Eq  (4.20)  now  takes 


the  form 


VV> 


J  ♦  gfv^kfv')  hf)c(v') 


(4.23) 


where  g(w),  a  periodic  modal  controllability  matrix,  is 


g(V>)  =  F~  (V')B(y/) 


(4.24) 


Modes  rj  (y)  are  controllable  if  the  corresponding 

C  •  ^ 


are  nonzero . 


Assuming  X  and  X  are  the  pairs  of  the  Poincare 
1  «  A  W  |4 

exponents  associated  with  T)  and  7?_  ,  respectively,  and 

1 v  2  3»4 

assuming,  for  the  present  case,  that  k^(w)  =  k^  ( v7)  =  0,  Eq 

(4.23)  can  be  expanded  to  be 


Vv)  - 


~2  ?2 


k3g2 


k4*2 


^3+k3*3  °°3+k4®3 

u4+^l8l  ?.+lc  g 

4  m  4  4  4  4 


ri  (v'i 

c 


(4.25) 


By  inspection,  the  Poincare  exponents  X^  ^  are  the  same 
as  those  of  the  open-loop  system,  while  the  coupled  equations 

r 

for  the  modes  7>3  and  n ^  determine  the  Poincare  exponents  for 
X-3  The  equations  defining  and  can  be  decoupled  from 

Eq  (4.25)  and  separated  to  form  the  two-dimensional  system 


,  £,+k,g,  w,+k  g 

77  (y0  =  3  3  3  3  4  3  } 

L  a)4  +  k3«4  f4+k4g4  J 


(4.26) 


It  is  desired  to  choose  gains  k3(y0  and  k^(y/)  such  that 

the  Poincare  exponents  X  are  stabilised.  First,  however, 

o ,  A 


g^(V')  has  to  be  expressed  in  a  convenient  manner. 


Clearly , 


this  function  is  periodic;  it  may  be  expanded  in  a  Fourier 
series  as  (64:81-82) 


*i<v0  =  * 


S.h, 


coa(nV)  ♦  g .  sin(n^)  (4.: 

J 


where  the  Fourier  coefficients  are  obtained  from  g^(y<)  as 
follows  (2:153-154): 


io~  2n  i 


(W)  dy> 


=  fg. 


(yj)  cos  (ny>)  dy/ 


(4.28) 


n  J 


( V')  sin  (nyO  d^ 


Consequently,  assuming  that  and  are  either 

constant,  or  can  be  expressed  in  the  form 


k^(v>)  =  k^sininv)  or  k^(v>)  =  k^cos  (nyO 


(4.29) 


Eq  (4.26)  may  be  written  as 


D  (v>)  =  Aj'Oj,  (V')  ♦  A2(v')r>c(v') 


(4.30) 


where  A^  is  a  constant  matrix,  and  is  a  purely  periodic 
matrix.  Even  though  the  gains  k^  and  k^  may  be  chosen  such 
that  A j  has  stable  eigenvalues,  the  stability  of  the  system 
is  not  assured  (62:673-674).  However,  Calico  and  Wiesel 
(43:61-62;62:674)  developed  a  technique  which  can  be  used  to 
set  the  sum  of  the  real  parts  of  the  two  new  roots  to  any 
desired  value.  This  technique,  which  determines  the  required 
values  for  k^  and  k^ ,  is  reviewed  below. 

Defining  D(v/)  =  det[$(y,0)]  and  tr(°)  to  be  the  trace, 
it  can  be  shown  that  (62:674) 


D  (y)  =  tr 


[a(vO  j 


D(y/) 


(4.31) 


This  first  order  ordinary  differential  equation  may  be 
integrated  by  means  of  an  integrating  factor  to  yield 


D(y)  =  D (0) exp 


J  tr|^A(v/)j< 


(4.32) 
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Evaluating  Eq  (4.32)  at  v/  =  2n  and  recalling  that  the 


determinant  of  a  matrix  is  the  product  of  its  eigenvalues, 

2n 

A.  =  exp]  | 

=  3 


4  2n 

det  (2rc  ,0)  j  =  n  =  exp-  Jtr  |a(v/)  jdy/ 


(4.33) 


But,  the  eigenvalues  of  the  stonodromy  matrix  are  related 
to  the  Poincare  exponents  by  Eq  (4.7)  so  that 


4  2tt 

Z  \  =  h  |tr(A^)]d^ 

i  =  3  0 


(4.34) 


Substituting  for  A(v/)  from  Eq  (4.26)  yields 

2  n 

V?4  =  Wib  J  [k3(V,)g3(V/)  +  V¥',*4(V')]dV'  (4.35) 

oJ 


The  products  kgtv'lg.^V')  and  k^  (\p)  (y/)  can  be  assumed  to 

contain  constant  terms.  For  instance,  if  k^  is  constant  and 
g3(v0  has  a  constant  term  in  its  Fourier  expansion,  k^g^(v/) 
will  have  a  constant  term  and  a  series  of  periodic  terms. 
If,  on  the  other  hand,  k  (y/)  has  the  form  of  Eq  (4.20)  ,  the 
product  k^vOg^v')  will  generate  a  constant  term  via  the 
standard  trigonometric  identities  (62:673).  Thus,  Eq  (4.35) 
becomes 


♦  z'.  =  ♦  ?.  +  [k  g  ]  +  [k.g.]  (4.36) 

3  4  3  4  330  440 

where  [k_g_]  and  [k.g.]  are  constant  values. 

u  w  O  ^  %  O 

Vote  that  if  only  one  root,  say  f j,  is  unstable.  Calico 
and  Wiesel  (2:673)  showed  that  the  root  shifts  along  the  real 


axis  as  a  linear  function  of  the  gain  k  : 

3 

^3  ”  ^3  +  *k3*3"*o  (4.37) 

For  stability  it  is  necessary  that  the  sum  of  the  real 
parts  of  the  desired  Poincare  exponents  in  Eq  (4.36)  be 
negative.  Sufficiency,  however,  requires  that  each  Poincare 
exponent  has  a  negative  real  part.  Eq  (4.36)  can  be  used  to 
define  pairs  of  and  k^  to  yield  a  specific  sum,  and,  using 
these  values,  the  closed  loop  system  can  be  solved  by  Floquet 
analysis  to  find  the  actual  values  of  £  and  £  . 

Several  methods  exist  to  choose  values  of  k,  and  k . .  In 
the  case  where  g^  and  g^  have  constant,  non-zero  terms  g^Q 

and  g.  ,  respectively,  in  their  Fourier  series,  and  k_  and  k „ 

4o  3  4 

are  chosen  as  constant  gains,  one  method  is  to  simply  chose 
values  of  k^  and  k^  such  that  Eq  (4.36)  is  satisfied. 

Values  for  k  and  k  may  also  be  found  using  the 
equations 

k3  =  (?3  "  ^3}  /®3o  ;  k4  =  [?4  ?4]/g4o  (4.38) 

A  third  method  is  to  choose  a  value  for  k^  and  solve  the 
following  equation  for  k^ : 

k4  ‘  ?4o[  *  -  k3*3o]  M'3B) 

where  =  £3  ♦  £  -  £  - 

Different,  non-constant  gains  may  be  used  to  select 


higher-order  coefficients  from  the  Fourier  expansions  of 


e 


t 


c* 


w 


4 


4 


g  (y/)  and  g  (yj)  .  For  instance,  if  k  is  chosen  to  be  a 

u  4  vi 

constant,  and  k^(v/)  =  k^sin(v')  ,  then  the  analogue  of  Eq 
(4.36)  is 

*  ■  k3*3o  *  K*«Sl  (4'40> 

where  g.  is  the  coefficient  of  the  first  sine  term  in  the 
40 1 

Fourier  series  expansion  of  g4(y>).  Similarly,  for  k^Cv7)  = 

k  cos(v') 

*  =  k3g3o  +  2k4g4c1  (4.41) 

where  g„  is  the  coefficient  of  the  first  cosine  term  in  the 
4cl 

Fourier  series  expansion  of  g  (.yj)  .  The  values  of  k  and  k 

4  o  4 

can  still  be  calculated  with  one  of  the  three  methods  used 

when  both  k  and  k  are  constant. 

u  4 

In  terms  of  the  physical  coordinates,  x(v')  and  u(v')  ,  the 
control  required  to  change  the  two  unstable  modes  is  given  by 


u(v>)  =  j^k ^  ( )  f 

where  f  1  (y/)  and  f  .1(V')  are  the  third  and  fourth  rows  of  the 
F  1(v/)  matrix,  respectively.  The  gains  k^(v/)  and  are 

based  on  the  Fourier  series  for  g3(y>)  And  g^  (y/)  ,  as 
previously  discussed. 

To  summarise,  designing  a  scalar  controller  to  shift  a 
pair  of  unstable  roots  requires  the  following  steps: 

1)  through  numerical  integration  calculate  the  state 


-l 


-l 


(y)  +  k.(v)f.  (¥<)  x(y/) 


(4.42) 
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transition  matrix  at  the  end  of  one  period  to  determine 

the  Poincare  exponents  and  the  associated  eigenvectors 

* 

2)  choose  the  desired  Poincare  exponents 

3)  form  the  controllability  matrix,  g(v').  and  the 
associated  Fourier  series  expansions 

4)  choose  the  gains  k^(y)  as  either  constants  or  in 
the  form  of  Eq  (4.29)  and  determine  their  values  with 
one  of  the  three  methods  previously  described 

5)  calculate  the  control  using  Eq  (4.42),  and  insert 
it  into  the  feedback  control  system  of  Eq  (4.17) 

6)  integrate  the  controlled  system’s  state  transition 
matrix  for  one  period,  and  calculate  the  new  Poincare 
exponents  to  determine  the  individual  root  locations. 

7)  If  the  individual  root  locations  are  not 

acceptable,  pick  another  pair  of  values  for  the  gains 
which  satisfy  Eq  (4.36)  and  repeat  step  6. 

Using  the  technique  summarized  above,  two  unstable 

Poincare  exponents  can  be  shifted  into  the  stable  region 
while  the  other  exponents  remain  the  same.  Even  though  this 
analysis  was  accomplished  using  a  fourth  order  system  of 
equations,  the  sum  of  an  arbitrary  number  of  roots  may  be 
set.  However,  as  the  number  of  roots  increases,  determining 
the  gain  is  complicated.  Consider  the  coupled  rotor/fuselage 
system  discussed  in  Chapter  III,  where  the  dimensions  of  the 
state  vector,  x(v/).  are  of  order  20.  It  is  entirely  possible 
for  this  system,  when  uncontrolled,  to  generate  more  than  two 

unstable  modes.  If  m  is  the  number  of  unstable  modes,  then 


_  rccc 


k(^)  can  now  be  expressed  as 


k  (vO 


“  [  kl  ,k2 . km’°  ,0 . 0  ] 


(4.43) 


If  the  trace  rule  is  again  applied,  Eq  (4.36)  becomes 


1  <1  ■  5  (fi  *  'Wo] 


(4.44) 


i=  1  i  =  1 


Again,  for  stability  the  sum  of  the  real  parts  of  the 

new  Poincare  exponents  must  be  negative.  Eq  (4.44)  can  be 

used  to  calculate  k  after  values  for  k  ,k  , . . . ,k  .  have 

EQ  1  jL  m- 1 

been  specified.  However,  since  only  the  sum  of  the  roots  is 
negative  and  not  the  individual  values,  stability  is  not 
assured.  A  numerical  search  procedure  is  of  possible  use 
here,  but  attempting  to  find  m  different  values  for  the 
individual  components  of  the  gain  matrix  is  quite  difficult 
for  multiple  modes.  Consequently,  scalar  control  is 
difficult  to  apply  to  systems  with  more  than  two  unstable 


modes . 


In  this  situation,  vector  control  may  be  an 


alternative  solution  to  controlling  the  system. 


Vector  Control 


As  shown  in  the  previous  section,  scalar  control  can 
eliminate  a  system's  instabilities.  How,  suppose  the  control 
u(y)  has  more  than  one  dimension.  Consider  again  the  case  of 
a  fourth  order  system  (blade  flap-lag  motion)  with  and  77  ^ 
being  the  unstable  modes.  Bestricting  the  control  to  these 


V 


a 


two  inodes  yields 


u(v/)  =  1tlv)r)lw)  = 


k33  k34 

k43  k44 


(4.45) 


Just  as  in  the  scalar  case,  the  above  expression  for 


S 

fi 


u(^)  can  be  substituted  into  Eq  (4.20),  and,  since  the  stable 
nodes  are  left  unchanged,  the  closed-loop  equations  for  the 
controlled  inodes  are  (43  :  63  ;  62  :  674-675) 


V*0  = 


k33g33  +  k43g34  +  ?3  k34g33*k44g34+a>3  . 

k  g  +k  g  +w  k  g  +k  g  +£  Vv) 
33843  43844  4  34843  44844  ^4 


where  the  gain  elements  k.  (V')  are  all  functions  of  w  and  the 

V 

controllability  matrix,  g(y'),  is  periodic. 

The  two  oscillatory  modes  in  Eq  (4.46)  can  be 
transformed  into  a  pair  of  uncoupled,  purely  damped  modes  by 
decoupling  the  two  modes  in  question.  This  requires  that 


k33(V;)g43<V')  +  k43(V/)g44(Vy)  = 


k34<V/)g33(V/)  +  k44(V/)g34(V/)  = 


(4.47) 


The  desired  value  of  the  real  parts  of  the  Poincare 


exponents,  if,  .  ,  can  now  be  chosen,  and  the  diagonal  terms  in 
Eq  (4.46)  can  be  forced  to  assume  these  values.  This  results 
in  two  more  equations: 


k33(v)g33(^)  +  k43('t'*g34*v')  *>3  ~^3 


k34  (V/)  g43  (V)  +  k44  (^)  ®44  ^  ^4  ”  ^4 


(4.48) 


v'-.'  %'  s’ 


■  Va  a.v 


I5r»aw 


Thus,  Eqs  (4.47)  and  (4.4B)  constitute  four  equations 
with  four  unknowns--the  elements  of  the  matrix  k(y').  Since 
the  values  of  the  matrix  g(v>)  nay  be  calculated  at  evenly 
spaced  intervals  throughout  the  entire  period,  a  harmonic 
analysis  algorithm  may  be  generated,  and  a  Fourier  series 
representation  of  the  elements  of  kCv7)  can  be  obtained. 
However,  the  controllability  condition 

[g33(v/)g44(v/)  "  *34(v/)g43(v/)]  *  0  (4.49) 

must  be  satisfied  at  any  point  in  the  period,  or  infinite 
values  of  k(v')  will  be  generated  (62:675). 

If  k(v)  is  infinite  at  any  point,  the  gain  functions 
k  (v-)  cannot  be  used  to  explicitly  decouple  the  system.  On 

the  other  hand,  these  functions  can  be  chosen  to  obtain  the 
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largest  possible  shift  in  the  Poincare  exponents  while 
minimizing  the  required  control  forces  (43:60). 

To  accomplish  this,  assume  k.^(y')  and  k^(v')  are  zero 
and  then  apply  the  trace  rule  to  Eq  (4.46).  The  result  is  an 
equation  similar  to  the  scalar  control  equation  [Eq  (4.36)1: 

*  =  [k33(v/)g33(w)]o  + 

where  the  bracketed  terms  are  the  constant  parts  of  the 
Fourier  series  expansion  of  the  product  k(v/)g(v/).  Mow,  the 
exponential  Fourier  series  for  both  kg.j(v')  and  k^^(v/)  may  be 
written  as  (63:60-62) 


k44(V,)g44 


(4.50) 


4 

4 


Klip)  =  K.  2  k^e 

l=-CD 


(4.51) 


where  the  constant  K  will  be  chosen  later.  The  integral 
expression  for  the  coefficients  can  be  expanded  to  be 


^  Jx(V')  |cos((v)  ~  iain Ity) Jdy 


(4.52) 


Hence  , 


ko  ■  koo  :  k<  *  Kkc/  4k.J  :  k-<  *  Kv  ik*J 


(4.53) 


The  Fourier  expansion  for  either  gggtV7)  or  g^  (y>)  can  be 
similarly  stated  to  be 


UU 

o(V')  =  J  g. 


where 


yn.  /n." 


The  product  icty/lafv')  can  then  be  written  as 


>w«H  ■  ■$  l 


;=-0D  flV=-00 


(4.54) 


*o  ■  gc  ■  =  K«c  -  tg«  )  :  *-«.  ■  Kgc„*  ig«J  ,4- 

v  ffX  ffir  v  ***- 


(4.56) 


Mow,  let  m  =  -t  ,  yielding  the  constant  tern 


coefficient 


k  ( v7)  o  (v) 


1  =  K?  k,g  , 


(4.57) 


If,  for  the  minimum  control  condition,  the  sum  of  the 

squares  of  the  coefficients  k  and  k  are  constrained  to  be 

cl 

unity,  then 


?2  K*  k*J  -  1  =  2  2  V-/ +  ko  =  1  (4  58) 

Z=-00  *  £=-00  £  =  1 

Maximizing  the  constant  term  |^K(v/)a(V/)J  while 

minimizing  the  control  k(v')  can  be  viewed  as  an  optimization 
problem,  where  the  Lagrangian,  L,  can  be  expressed  as  (65:55) 


L  ■  f  *  x[2  *  ko] 


(4.50) 


l--ao 


with  X  being  the  Lagrangian  multiplier. 


Applying  the  necessary  condition  for  stationarity 


results  in 


Se.  ■  3k.  JE  k£*-£  +  2X.f !***-*  *  -.1  = 

j  J  [£5*-oo  £=  1  J 


or,  for  j  =  0, 


k  =  -  — 
o  2X 


and,  for  j  *  0, 


(4.60) 


(4.61) 


k,  ■  -  r 


(4.62) 


Substituting  Eqa  (4.61)  and  (4.62)  into  the  constraint 


yields 


e 


=  [  Tgo  +  2  1  *z*-z  ] 


(4.63) 


If  D  is  defined  to  be  D  =  2\  ,  then 


2 

gC  + 

o 


•.to*  or 


(4.64) 


and  the  coefficients  of  k(v')  become 


k  =  -  — g  ■  k  =  -  —  ■a  •  k  =  -  — a 

c  Dbc  •  c,  D8c,  ’  s,  D8s, 

o  o  Z  Z  Z  Z 

The  constant  term  ^K(y)a(y/)j  is  then 

[K(Wo(v-)]o  =  -  £[  |  («c/  *ij] 


(4.65) 


(4.66) 


The  parameter  K  is  still  available  for  pole  placement. 
Consequently,  if  the  Fourier  coefficients  for  633(V/)  and 
g44(v/)  are  g3c/*  g3s^  and  g4c^'  g4s/  resPecti  veiy ,  then  Eq 
(4.50)  becomes 


•  ■  •  Ui.,  •  .y] 


(4.67) 


Choosing  and  so  the  new  Poincare  exponents  are 
placed  in  their  desired  locations  requires  a  numercial  search 
procedure  similar  to  that  used  with  scalar  control.  In  other 
words,  a  value  for  K  may  be  chosen,  and  then  the  following 

WO 


equations  can  be  used  to  find  K44: 


(4.68) 


Once  K  and  K  have  been  calculated,  the  resulting 

JO 

optimal  gain  functions,  k^^(v^)  and  k^(v'),  maximize  the  shift 
in  the  Poincare  exponents  with  the  smallest  possible  control 
u  (v')  ,  where 


u(y>) 


k33(V/)f31(w) 
k44(^) f~1(v) 


x(v>) 


(4.60) 


Consequently,  a  vector  controller  can  be  designed  as 
follows,  with  steps  (1)  through  (3)  being  identical  to  those 
steps  used  in  the  scalar  control  case: 

4)  determine  the  gain  matrix  k(w)t  either  by  solving 
four  linear  equations  with  four  unknowns,  or,  if  the 
controllability  condition  is  violated,  by 

a)  choosing  a  value  for 

b)  using  Eq  (4.68)  to  solve  for  the  constant  K44 

c)  calculating  the  optimal  gain  functions  k33(v') 
and  k44(v«)  with  Eqs  (4.64)  and  (4.65) 

5)  calculate  the  control  required  in  terms  of  the 
physical  coordinates  u(y>)  and  x(yO  and  insert  it  into 
the  feedback  control  system  of  Eq  (4.17) 

6)  verify  the  root  shift  by  integrating  the 
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controlled  system’s  state  transition  matrix  for  one 

/ 

period  and  obtaining  the  new  Poincare  exponents. 

/ 

Iterate  to  obtain  the  desired  Poincare  exponents. 

Following  the  above  procedure  will  result  in  developing 
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a  vector  controller  which  will  shift  two  unstable  Poincare 
exponents  while  leaving  the  others  unaltered.  However,  it 
may  be  necessary  to  stabilize  more  than  two  modes  in  a 
system.  Again  consider  using  the  modal  control  technique  to 
stabilize  a  system,  such  as  the  coupled  rotor/fuselage 
system,  with  more  than  two  unstable  modes.  Unfortunately  the 
same  situation  which  occurs  with  the  scalar  controller  also 
arises  when  the  vector  controller  is  used.  That  is,  it 
becomes  extremely  difficult  to  determine  the  correct  gains 
required  to  shift  more  than  two  unstable  modes  to  their 
desired  locations. 

To  summarize,  both  scalar  and  vector  controllers  have 
been  developed  to  shift  two  unstable  nodes  into  the  stable 
region.  It  appears,  though,  that  the  only  way  to  shift  more 
than  two  modes  is  by  guessing  the  values  of  the  individual 
terms  in  the  gain  matrix.  However,  another  possible  solution 
is  to  apply  the  modal  control  technique  to  shift  a  single 
pair  of  unstable  roots  and  then,  using  that  controlled 
system,  design  a  control  which  shifts  two  other  unstable 
roots.  This  type  of  control  could  then  be  used  to  stabilize 
any  number  of  modes.  All  that  remains  is  to  verify  these 
observations  by  attempting  to  control  a  system  of  equations 
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V.  Besul ts  from  Controlling  Blade  Flap-Lag  Motion 


Chapter  IV  outlined  the  procedure  used  to  develop  modal 
control  theory.  The  results  of  this  theory  can  now  be 
implemented  in  a  controller  to  reduce  or  eliminate  a 
helicopter  blade's  unstable  regions. 

The  first  step  in  constructing  a  controller  is  to  choose 
the  control  system.  Calico  and  March  (42)  applied  the  time 
periodic  modal  control  technique  to  the  problem  of 
controlling  the  flapping  instabilities  of  a  helicopter  blade 
using  a  flap  torque  actuator  situated  at  the  blade  root. 
Calico  and  Wiesel  (43)  designed  a  time  periodic  modal  control 
system  which  used  existing  collective  and  cyclic  pitch 
mechanisms  on  a  conventional  swashplate.  Stabilizing 
helicopter  blades  with  this  type  of  active  control  could 
possibly  eliminate  the  need  for  mechanical  lead-lag  dampers 
and  expand  a  helicopter's  operating  envelope  (31:13-14).  It 
remains  to  be  seen  what  computational  difficulties  arise  when 
the  modal  control  technique,  as  implemented  by  Calico  and 
Wiesel,  is  used  to  control  more  complex  systems,  such  as 
those  described  in  Chapters  II  and  III.  Therefore,  this 
technique  will  first  be  used  to  control  the  flap-lag  motion 
of  a  single  rotor  blade. 

The  control  vector  u(v')  for  the  control  of  a  helicopter 
blade’s  motion  typically  has  components  which  represent 
collective  pitch,  9^ ,  and  cyclic  pitch,  6^,  6*: 
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The  control  matrix  B(v')  is  then  composed  of  periodic 
pitch  control  functions  : 
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The  above  expressions  are  generated  by  extracting  the 
collective  and  cyclic  pitch  terms  from  the  generalized 
aerodynamic  forces  of  Kqs  (2.72)  and  (2.73).  Terms  Second 
order  and  higher  are  neglected. 


Results  Using  Scalar  Control 


For  scalar  control,  u(v>)  can  come  from  either  the 
collective,  cyclic  cosine,  or  cyclic  sine  pitch  controls. 
Designing  a  scalar  control  system  requires  determining  the 
gain  matrix,  k(V'),  needed  to  move  the  unstable  modes  to  their 
desired  locations.  To  obtain  the  gains,  Eq  (4.36)  must  be 
solved.  The  modal  controllability  matrix  is  given  by 


gj (VO  = 
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where  j  =  o ,  c ,  or  s  for  collective,  cyclic  cosine,  or  cyclic 
sine  pitch  control,  respectively.  The  elements  of  g^  (v>)  can 
be  written  in  a  Fourier  series  form  using  harmonic  analysis. 
The  gains  k  and  k  can  then  be  calculated  by  one  of  the 

u  4 

three  methods  described  in  Chapter  IV. 

These  gains  can  be  used  in  Eq  (4.42)  to  determine  the 
control  required  to  shift  the  roots  to  their  desired 
locations.  Once  this  control  is  inserted  into  the  system, 
the  Floquet  solution  confirms  the  new  Poincare  exponents. 

To  demonstrate  the  modal  control  technique,  five 
different  scalar  controllers  are  designed  to  stabilize  a 
point  in  the  unstable  region  at  fu  =  0.40  (Figure  0)  .  This 
design  point  is  defined  by  p  =  1.15  and  =  1.40.  With  = 
Cc  =  0.000,  the  real  parts  of  the  uncontrolled  Poincare 

exponents  are  ? 3  4  *  0.00189894.  The  real  parts  of  the 
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desired  pole  locations  are  chosen  to  be  f  =  -0.025.  The 

3,4 

five  scalar  controllers  are  as  follows: 

1)  collective  pitch,  using  the  constant  terms  of  the 

Fourier  series  expansion  of  g(y')  corresponding  to  the 

unstable  modes,  denoted  as  g_  and  g. 

O  O  40 

2)  cyclic  cosine  pitch,  using  the  constant  terms  of 

g(V').  denoted  as  c_  and  c. 

30  40 

3)  cyclic  sine  pitch,  again  using  the  constant  terms 
as  described  in  the  first  controller  and  denoted  as 


s_  and  8. 
3o  4o 


4) 


collective 


constant/cosine  pitch. 


using  the 


constant  term  of  the  Fourier  expansion  of  g(y>) 
corresponding  to  the  first  unstable  mode  (g_  )  and  the 

JO 

coefficient  of  the  first  cosine  term  in  the  Fourier 

aeries  expansion  of  g(y')  corresponding  to  the  second 

unstable  mode  (g.  ) 

4C1 

5)  collective  cons tant/sine  pitch,  using  the  terms 

described  in  the  fourth  controller,  except  ,  the 

4s  i 

term  corresponding  to  the  second  unstable  mode,  is  the 
coefficient  of  the  first  sine  term  of  g(v')- 
The  values  for  and  are  calculated  by  using  all 

three  methods  described  in  Chapter  IV.  The  first  and  third 
methods  of  choosing  the  gains  require  iteration  procedures  to 
find  values  for  k  and  k  such  that  the  control  succeeds  in 

ij  4 

shifting  the  unstable  roots  to  their  desired  locations.  The 
accuracies  of  the  new  pole  locations,  when  compared  to  the 
desired  values,  depend  on  how  many  iterations  are  used  to 
find  the  required  gains.  Only  a  few  iterations  are  necessary 
if  a  low  degree  of  accuracy  is  desired.  However,  placing  the 
poles  with  increasing  accuracy  requires,  in  general,  an 
increasing  number  of  iterations. 

On  the  other  hand,  the  second  method  requires  only  one 
iteration  to  arrive  at  values  of  k^  and  k^ ,  and  the  accuracy 

of  the  resulting  root  locations  is  quite  impressive.  Four  of 
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the  five  different  controllers  tested  produce  new  Poincare 
exponents  accurate  to  at  least  four  decimal  places  when 
compared  to  their  desired  locations.  This  corresponds  to  a 
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TABLE  III 


Controller  Poincare  Exponents,  (C. 


0.000) 


TYPE  OF 
CONTBOLLEB 

GAIN  VALUES 

COEFFICIENTS 
FOB  g(v>) 

PO INC ABE 
EXPONENTS 

COLLECTIVE 

PITCH 

k3=-0. 1 17589 
k4=  1.008789 

g,  =  0.228753 

30 

gA  =-0.266650 
4o 

-0.025062 

CYCLIC 

COSINE 

PITCH 

k  =-0. 142515 

3 

k4=  0.141415 

c„  =  0.188744 

dO 

c,  =-0.190212 

4o 

-0.025000 

CYCLIC 

SIHE 

PITCH 

k3=-0. 116426 
k4=-0. 187789 

s_  =  0.231038 

dO 

s,  =  0.143241 

4o 

-0.024805 

COLLECTIVE 
CONSTANT/COS I HE 
PITCH 

k  =-0. 117589 

3 

k4=  0.135280 

g„  =  0.228753 
30 

g,  =-0.397678 
4°1 

-0.025010 

COLLECTIVE 

CONSTANT/SINE 

PITCH 

k3=-0. 1 17589 
k4=-0. 173313 

g_  =  0.228753 

dO 

g^  =  0.310409 
48  1 

_ i _ 

-0.024994 

99.75  percent  accuracy.  The  fifth  controller,  cyclic  sine 

pitch  control,  produces  results  accurate  to  within  99.20 

percent  of  the  desired  pole  locations.  These  results  are 

consistent  for  a  wide  variety  of  flight  conditions  tested. 

As  an  example  (refer  to  Table  III)  ,  if  the  real  parts  of  the 

desired  pole  locations  are  .  =  -0.025,  the  cyclic  sine 

3,4 

f 

pitch  control  yields  Poincare  exponents  with  real  parts  at 
£ ,  .  =  -0.024805,  whereas  the  other  controllers  produce 

3,4 

results  better  than  =  -0.025062.  It  should  be  noted 

3 , 4 

that,  in  all  the  results  obtained,  the  stable  modes  are  left 
unchanged . 

How,  each  feedback  controller  developed  is  designed  to 
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stabilize  a  pair  of  Poincare  exponents  at  a  specific  flight 
condition.  Similar  controllers  could  be  produced  for  any 
flight  condition.  This,  however,  would  require  changing  the 
feedback  control  gains  as  the  parameters  are  varied. 
Therefore,  it  would  be  advantageous  for  a  flight  controller, 
designed  to  control  the  system  at  a  single  point,  to  be  able 
to  reduce  or  eliminate  blade  instabilities  at  off-design 
conditions.  However,  it  should  be  emphasized  that  no 
guarantee  of  stability  exists  at  equilibrium  conditions  other 
than  those  the  controller  is  designed  for. 

Consequently,  the  nominal  design  point  for  each  of  the 
five  scalar  controllers  is  chosen  to  be  at  ^  =  0.40,  with  p  = 
1.15,  w^  =  1.40,  and  =  0.000.  As  before,  the 

controllers  are  designed  to  shift  the  real  parts  of  the 

/  / 

unstable  Poincare  exponents  to  £  =  -0.025.  With  each 

3 , 4 

controller  set  to  shift  the  unstable  roots  at  the  chosen 
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design  point,  flight  conditions  are  varied  and  new  Poincare 
exponents  are  generated. 

Figure  18  illustrates  the  five  flight  controllers' 
performances  when ,  with  the  advance  ratio  held  constant  at  /li 
=  0.40,  p  and  are  varied.  Although  convergence  is  not 
guaranteed  for  off-design  conditions,  it  can  be  seen  from  the 
figures  that,  using  this  point  design,  each  controller 
eliminates  the  unstable  regions  over  a  wide  range  of  p  and 
combinations.  In  fact,  two  of  the  controllers--collective 
constant/cosine  and  collective  constant/sine--el iminate  the 
entire  unstable  region  for  ij  =  0.40.  The  other  three 
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(a)  (b) 

Figure  18.  Stability  Regions  at  p  =  0.40 

control lers--col lective ,  cyclic  sine,  and  cyclic 

cosine--generate  small  bands  of  unstable  regions  near  the 
lower  end  of  co^.  .  The  bands  produced  by  the  cyclic 

controllers  extend  past  p  =  1.75  [Figure  18(a)],  while  the 
collective  pitch  controller  eliminates  the  unstable  region 
for  p  ^  1.30  [Figure  18(b)].  Note  that  these  bands  of 
unstable  regions  are  much  smaller  in  area  than  the 
uncontrol led  case.  Indeed,  the  controllers,  using  this 

single  point  design,  significantly  reduce  the  unstable  region 
for  v  =  0.40.  It  should  be  mentioned  that  other  point 

designs  were  not  examined  in  off-design  conditions. 

With  the  same  five  scalar  controllers  in  place, 

stability  in  hover  is  considered  as  a  function  of  p  and  . 


Figure  10  compares  the  controlled  and  uncontrolled  stability 
boundaries  for  hover.  Figure  10(a)  shows  that  the  collective 
constant/cosine  and  constant/sine  pitch  controllers  now 
produce  an  unstable  region  extending  past  p  =  1.75.  This 
narrow  band  is  similar  to  the  region  generated  by  both  the 
cyclic  controllers  for  the  M  =  0.40  case  [Figure  16(a)). 
Figure  10(b)  illustrates  the  results  of  using  the  cyclic 
cosine  and  cyclic  sine  pitch  controllers.  Once  again,  narrow 
bands  of  unstable  regions  are  generated,  but  the  cyclic 
cosine  controller  eliminates  this  unstable  region  for  p  £ 
1.245  and  the  cyclic  sine  controller  eliminates  the  band  for 
p  >  1.285.  It  is  interesting  to  note  that  these  results  are 
similar  to  those  obtained  when  the  collective  controller  is 
used  at  fj  =  0.40  [Figure  18(b)). 

On  the  other  hand,  in  hover  the  collective  pitch 
controller  produces  results  very  much  different  from  those 
previously  summarized.  With  the  first  four  controllers  the 
unstable  regions  are  greatly  diminished  when  compared  to  the 
uncoupled  case.  But  now  Figure  19(c)  shows  that  when  the 
collective  controller,  designed  at  p  =  0.40,  is  used  in 
hover,  a  large  unstable  region  is  generated.  In  fact,  this 
region  is  larger  than  that  generated  by  the  uncontrolled 
system. 

Another  way  of  observing  how  a  point-designed  controller 
behaves  in  off -design  cases  is  to  hold  p  and  constant  at 
1.15  and  1.40,  respectively,  and  vary  the  advance  ratio  from 
)j  =  0.00  to  )J  =  0.50.  Figure  20  summarizes  the  results  for  a 
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system  with  £  ^  =  (;  =  0.000.  As  can  be  seen  from  the  figure, 
the  uncontrolled  blade  is  unstable  throughout  the  entire 
flight  regime  examined.  Each  controller,  designed  at  = 
0.40  to  shift  the  real  parts  of  the  unstable  roots  from 


=  0.00189894  to 


-0.025,  accomplishes  this  task. 


Further,  four  of  the  five  controllers  generate  stable  systems 
from  t-i  =  0.00  through  tu  -  0.50.  As  the  advance  ratio 
increases,  the  cyclic  sine  and  cosine  pitch  controllers 
[Figure  20(a)]  produce  fairly  constant,  but  increasingly 


stable,  values  for  the  new  Poincare  exponents. 


The  same 


observations  can  be  made  for  the  collective  constant/cosine 
and  constant/sine  pitch  controllers  [Figure  20(b)],  although 

x 

at  lower  advance  ratios  the  new  Poincare  exponents  are  less 
stable  than  those  generated  by  the  cyclic  pitch  controllers. 

On  the  other  hand,  the  collective  pitch  controller 
[Figure  20(b)]  produces  results  which  differ  from  the  other 
four  controllers.  Even  though  the  unstable  roots  are  shifted 
to  their  desired  locations  at  (J  =  0.40,  the  controller 

generates  unstable  modes  below  advance  ratios  of  0.06.  In 
fact,  below  fj  -  0.05,  the  controlled  system  is  more  unstable 
than  the  uncontrolled  system.  Further,  above  )J  =  0.40  the 

x 

pair  of  complex  conjugate  Poincare  exponents  crosses  the  real 
axis  and  splits  off  as  two  real  roots,  although  they  remain 
stable . 

It  is  also  interesting  to  examine  the  effect  of  adding 
passive  damping  to  the  system’s  lag  hinge.  To  accomplish 
this,  Cr  is  set  to  0.001,  and,  for  a  design  point  of  p  =  1.15 
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TABLE  IV 


Controller  Poincare  Exponents  =  0.001) 


TYPE  OF 
CONTROLLER 

GAIN  VALUES 

COEFFICIENTS 
FOB  g(^) 

POINCARE 

EXPONENTS 

COLLECTIVE 

PITCH 

k3=-0. 109599 
k4=  0.523310 

g3o=  0.232754 

g.  =-0.048746 

40 

-0.024980 

CYCLIC 

COSIME 

PITCH 

k3=-0. 144719 
k4=  0.119706 

c  =  0.176269 
oo 

c.  =-0.213101 

4o 

-0.024999 

CYCLIC 

SINE 

PITCH 

k3=-0. 101668 
k4=-0. 202978 

8-  =  0.250909 
oo 

s.  =  0.125676 

4o 

-0.024810 

COLLECTIVE 
CONST ANT /COSINE 
PITCH 

k3=-0. 109599 
k4=  0.115905 

g_  =  0.232754 
Oo 

g  =-0.440178 
*C1 

-0.024995 

COLLECTIVE 
CONSTANT/ SINE 
PITCH 

k3=-0. 109599 
k4=-0. 183796 

g3o=  0.232754 

g,  =  0.277584 
48 1 

-0.025005 

controllers  is  evaluated  at  off-design  conditions;  for  this 
case,  at  varying  values  of  the  advance  ratio.  These  results 
are  tabulated  in  Figure  21  and  compared  to  the  uncontrolled 
system.  With  -  0.001,  the  uncontrolled  system  generates 
stable  Poincare  exponents  for  y  5  0.175.  Above  this  point, 
however,  the  uncontrolled  system  generates  unstable  modes. 
The  controllers  produce  results  and  trends  almost  identical 
to  the  case  where  =  0.000.  Moreover,  even  the  collective 
pitch  controller  generates  stable  modes  throughout  the 

spectrum  of  flight  conditions  examined.  But,  below  ij  =  0.20 

✓ 

the  Poincare  exponents  are  not  as  stable  as  the  uncontrol led 
system,  and  the  collective  pitch  controller  generates  complex 
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conjugate  pairs  of  stable  Poincare  exponents  through  p  = 
0.50. 

The  transient  responses  of  the  system  resulting  from  an 

initial  disturbance  can  be  determined  assuming  unit  initial 

#  • 

values  for  each  of  the  four  states,  /?,  C.  0  ,  and  C  A 

collective  pitch  controller,  designed  at  p  :  0.40,  with  C  ;•  = 

0.001,  p  =  1.15,  and  =  1.40,  is  used.  Figure  22 

illustrates  the  transient  responses  after  ten  periods. 
Figure  22(a)  indicates  that  both  the  uncontrolled  and 
controlled  responses  of  ft  to  a  step  input  are  similar.  The 
responses  damp  out  rapidly  in  a  periodic  manner.  As  for  the 
lag  modes’  responses.  Figure  22(b)  shows  that  the 
uncontrol led  system  produces  a  very  slightly  unstable 
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response  while  the  controlled  system  damps  out  the  transients 
fairly  quickly.  Figure  22(c)  plots  the  magnitude  of  the 
control  response  as  a  function  of  the  azimuth  angle,  and 
shows  that  the  required  control  is  decreasing.  It  should  be 
noted  that  if  a  rotor's  nominal  rotation  rate  is  730 
revolutions  per  minute,  a  blade  completes  ten  revolutions  in 
0.822  seconds.  Thus,  the  scalar  collective  pitch  controller 
stabilizes  the  system  in  a  very  short  period  of  time. 

To  verify  the  above  observations,  the  other  four 
controllers  can  also  be  used  to  stabilize  the  system  at  the 
same  design  point  for  ten  periods,  and  the  flap,  lag,  and 
control  transient  responses  recorded.  The  responses  of  each 
controller  are  all  very  similar.  Figure  23  shows  the 
transient  responses  for  the  flap  and  lag  modes  after  ten 
periods  for  two  of  the  control lers--col lective  constant/sine 
pitch  and  cyclic  cosine  pitch.  As  can  be  seen  by  comparing 
these  figures  with  the  responses  generated  by  the  collective 
pitch  controller  (Figure  22) ,  the  controllers  do  not 
significantly  alter  the  stable  roots’  responses  (i.e.,  ft), 
but  they  effectively  stabilize  the  lag  modes’  responses  in  a 
short  period  of  time. 

Figure  24  compares  the  control  required  during  ten  of 
the  blade *8  revolutions  for  three  of  the  controllers: 
collective;  collective  constant/sine;  and  cyclic  cosine. 
Clearly,  the  required  control  dramatically  decreases  as  the 
number  of  revolutions  increases. 

It  should  be  noted  that  the  trends  shown  by  the 


Figure  24.  Transient  Responses  for  Collective,  Cyclic 
Cosine,  and  Collective  Constant/Sine  Controls 

collective  constant/cosine  and  cyclic  sine  pitch  controllers 
mirror  those  recorded  by  the  collective  constant/sine  and 
cyclic  cosine  controllers,  respectively.  In  addition,  the 
trends  observed  are  repeated  when  is  changed  from  0.001  to 
0.000.  Therefore,  each  of  the  five  scalar  controllers 
stabilizes  the  blade  flap-lag  motion. 

The  above  results  clearly  demonstrate  the  success  of 
controlling  blade  flap-lag  motion  using  scalar  control.  The 
next  step  is  to  see  if  a  vector  controller  can  be  as 
effective . 
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Reaul is  Using  Vector  Control 

A  vector  controller  can  be  designed  in  much  the  same 
manner  as  a  scalar  controller.  However,  in  this  case  the 
multiple  elements  of  the  control  vector,  u(vO  ,  will  be  used. 
For  the  current  study,  both  cyclic  sine  and  cyclic  cosine 
pitch  control  will  be  vised.  Hence, 


u(y) 


v*> 

V*> 


(5.5) 


Again,  determining  the  gain  matrix  k(y/)  requires 
knowledge  of  the  modal  controllability  matrix  g(y/).  Eq 
(4.24)  shows  this  matrix  to  be 
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Using  the  harmonic  analysis  technique,  values  for  g.  .(yO 
may  be  calculated  at  any  given  azimuth  angle.  If  and 

are  the  unstable  modes,  only  the  last  two  rows  of  the  g(v/) 
matrix  are  of  interest.  Given  that  the  controllability 
condition  is  satisfied  throughout  the  period,  the  four  linear 
equations  described  in  Eqs  (4.47)  and  (4.48)  can  be  solved  at 
evenly  spaced  intervals  to  obtain  the  Fourier  series 
expansions  for  the  four  components  of  the  gain  matrix  k(v'). 
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With  =  C £  =  0.000,  a  design  point  of  p  =  1.15  and  co^. 
=  1.40  at  p  =  0.40  is  chosen.  As  can  be  seen  from  Figure  9, 

this  point  is  in  the  unstable  region,  and  the  real  parts  of 

✓ 

the  unstable  Poincare  exponents  are  £  =  0.001898936.  The 

3,4 

real  parts  of  the  desired  pole  locations  are  again  chosen  to 
be  =  -0.025. 

In  calculating  the  gain  matrix,  k(v0,  it  mas  discovered 
that  the  controllability  condition  is  violated  twice  in  the 
interval  between  0  and  2 n .  Hence,  infinite  values  of  k(y/) 
are  obtained,  and  the  vector  controller  fails  for  this 
situation.  It  should  be  noted  that  Calico  and  Wiesel  (43:63) 
documented  similar  behavior  in  the  investigations  of  an 
isolated  rotor  blade’s  flapping  motion. 

A  second  vector  controller  was  designed  so  that  the  gain 
matrix  maximizes  the  constant  terms  in  Eq  (4.50).  The 

Fourier  expansion  of  the  modal  controllability  matrix,  g(v>)  , 
provides  the  Fourier  coefficients  needed  to  calculate  the 
elements  of  [k  (v>)  g,,  (v>)  J  and  [k .  .  (^)  g  .  .  (v>)  )  ,  as  defined 
in  Eqs  (4.64)  through  (4.66).  The  values  for  K33  and 

✓ 

required  to  move  the  real  parts  of  the  unstable  Poincare 

* 

exponents  to  =  -0.025  are  obtained  by  varying  K  and 

u  i  4  uj 

using  Eq  (4.68)  to  calculate  K44.  The  Fourier  coefficients 
of  833^)  and  £44  (V0  are  again  used  to  construct  the  Fourier 
series  representation  of  the  gain  matrix  k(y»)  ,  and,  after 

calculating  and  ,  this  matrix  is  used  in  the  feedback 
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control  system  to  move  the  new  Poincare  exponents  to  their 
desired  locations.  With  K33  =  0.21275545  and  J<44 
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r 


r 


V 


f 


* 

0.06116429,  the  real  parts  of  the  new  Poincare  exponents  are 
shifted  to  =  -0.025000  when  vector  control  is  applied. 
Agai  n,  all  other  modal  frequencies  are  left  unchanged. 

The  efficiency  of  the  numerical  search  procedure  used  to 
find  the  values  for  and  depends  upon  the  accuracy 
desired.  Obtaining  results  with  a  high  degree  of  accuracy 
requires  a  large  number  of  iterations  which,  of  course,  can 
be  a  time-consuming  endeavor.  Lowering  the  required  accuracy 
results  in  fewer  iterations  needed  to  obtain  K  and  K 

u  j  44  . 

In  this  case,  the  results  obtained  using  the  vector 
controller  are  much  more  accurate  than  those  results  obtained 
with  the  five  scalar  controllers  (refer  to  Table  III) ,  since 

K  and  K  are  obtained  after  14  iterations.  However, 

44 

recall  that  the  scalar  control  results  were  obtained  without 
performing  any  iterations.  Chapter  IV  outlined  a  numerical 
search  procedure  for  scalar  control  where  one  gain  is  chosen 
and  the  other  calculated,  with  the  process  being  iterated 
until  the  results  are  satisfactory.  With  this  procedure  the 
scalar  controllers  can  be  just  as  accurate  as  the  vector 
controller,  if  so  desired. 

The  performance  of  the  vector  controller  in  off-design 
cases  is  also  considered.  The  design  point  is  chosen  to  be 
identical  to  the  one  used  for  the  scalar  controllers  (p  = 
1.15,  =  1.40,  and  =  0.000  at  /j  =  0.40),  and,  with 
the  advance  ratio  fixed  at  (J  =  0.40,  p  and  are  varied. 
For  this  case,  the  vector-controlled  system  produces  stable 
results  with  every  possible  p  and  co^.  combination.  If  the 
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are 


same  controller  is  used  to  generate  results  as  p  and 
varied  in  hover,  stable  responses  are  again  produced 
throughout  the  flight  regime.  In  other  words,  the  vector 
controller,  designed  for  a  single  point,  completely 
stabilizes  the  system  for  all  values  of  p  and  in  hover  and 
at  m  =  0.40.  Comparing  these  results  to  those  generated  by 
the  scalar  controllers  (Figures  18  and  19),  leads  to  the 
conclusion  that  the  vector  controller  eliminates  a  wider 
range  of  unstable  p  verses  to  regions  than  any  of  the  scalar 


control lers 


The  vector  controller’s  performance  can  also  be 
evaluated  by  varying  the  rotor  blade’s  advance  ratio.  Thus, 
with  the  values  for  p  and  fixed  at  the  design  point  (p  = 
1.15,  =  1.40.  K p  =  C ^  =  0.000,  with  n  =  0.40),  the  advance 
ratio  is  varied  from  (J  =  0.00  to  fJ  =  0.50.  Figure  25(a) 
shows  that  the  vector  controller  generates  stable  modes 
throughout  the  flight  regime  inspected.  Just  as  with  the 


scalar  controllers,  the  vector  controller  produces  modes 
which  become  more  stable  as  the  advance  ratio  increases. 

With  =  0.001,  results,  similar  to  those  summarized 
above,  are  obtained.  For  K  =  0.21855978  and  K  = 
0.04988617,  the  vector  controller  produces  roots  with  real 
components  at  .  a  -0.025000.  In  other  words,  the  desired 
root  shift  is  achieved,  and  the  stable  roots  are  unaltered. 
With  the  design  point  set  at  p  =  1.15  and  =  1.40  for  tu  = 
0.40,  the  advance  ratio  is  varied  from  n  =  0.00  to  (J  =  0.50. 


The  results,  shown  in  Figure  25(b) ,  are  almost  identical  to 
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(a)  =  0.000  (b)  =  0.001 
Figure  25.  Controlled  Lag  Modes  as  a  Function  of  ij 

the  case  where  =0.000.  The  controller  generates  stable 
inodes  throughout  the  spectrum  of  flight  conditions,  and  these 
inodes  become  more  stable  as  the  advance  ratio  increases. 

Finally,  the  controlled  system's  transient  responses  are 
generated  for  ten  periods.  Once  again,  the  design  parameters 
are  as  follows:  p  =  1.15;  =  1.40;  /l<  =  0.40;  and  = 
0.001.  Figure  26(a)  is  the  flap  modes’  responses  after  ten 
periods,  while  Figure  26(b)  shows  the  lag  modes'  responses 
after  ten  periods.  Both  curves  are  similar  to  those 
generated  when  the  five  scalar  controllers  are  used.  The 
trends  show  that  the  stable  flap  responses  are  not 
significantly  altered,  and  the  vector  controller  stabilizes 
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the  lag  modes’  responses,  rapidly  damping  out  any  transients. 
Finally,  the  physical  control  required  by  the  vector 
controller  after  ten  periods  is  summarized  in  Figure  26(c). 

It  is  interesting  to  compare  the  two  figures  summarizing 
the  required  scalar  and  vector  control.  Figures  24  and  26(c). 
The  vector  controller  is  designed  to  maximize  the  shift  in 
the  unstable  roots  while  minimizing  the  required  control 
forces.  Indeed,  as  can  be  seen  from  the  figures,  the  control 
required  by  the  vector  controller  is  less  than  that  needed  by 
the  scalar  controllers,  though  not  significantly. 

In  the  same  vein.  Figure  27  compares  the  magnitude  of 
the  gain  matrix,  k(v'),  as  a  function  of  azimuth  angle  for 
several  of  the  controllers  examined.  Again,  the  collective 
constant/cosine  and  cyclic  sine  pitch  controllers  mirror 
their  respective  counterparts--col lect ive  constant/sine  and 
cyclic  cosine  controllers.  The  gains  for  the  collective  and 
cyclic  cosine  pitch  controllers  are  constants  while  the  gains 
for  the  collective  constant/sine  pitch  and  vector  controllers 
are  periodic  functions  of  the  azimuth  angle.  This  is  to  be 
expected  since,  for  the  constant/sine  controller,  k^fv7)  is  a 
function  of  the  sine  of  the  azimuth  angle,  according  to  Eq 
(4.29),  whereas  the  collective  and  cyclic  cosine  controllers 
use  constant  gains.  As  for  the  vector  controller,  the  gain 
k(w)  is  calculated  by  a  Fourier  series  expansion  as  a 
function  of  V',  and  it,  too,  is  periodic.  Note  that  the  gain 
required  by  the  collective  pitch  controller  is  greater  than 
any  of  the  other  controllers  examined.  This  observation 
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Figure  27.  Transient  Responses  of  the  Gain  Magnitudes 


mirrors  that  found  when  examining  Figure  24:  the  control 
required  by  the  collective  controller  is  greater  than  that 
needed  by  the  other  controllers.  These  observations  should 
be  the  same,  since  the  required  control  is  generated  by  the 
gain  k  (v')  . 

Finally,  to  insure  that  the  observations  summarized 
above  are  not  unique  to  a  single  design  point  or  flight 
condition,  the  scalar  and  vector  controllers  were  designed 
for  a  variety  of  desired  pole  locations  and  flight 
conditions.  In  each  case,  the  results  and  observations 
mirrored  those  described  above. 


Coupled  Botor/Fuselage  Motion 


Chapter  III  presented  and  verified  the  equations  of 
motion  for  a  coupled  rotor/fuselage  system.  Without  control 
the  uncoupled  system  generated  unstable  roots  for  both  the 
rigid  body  motions  and  the  flapping  motions  of  the  individual 
blades  on  the  main  rotor.  When  coupling  and  reversed  flow 
were  added,  similar  results  were  obtained.  Consequently,  to 
successfully  stabilize  the  coupled  rotor/fuselage  system  two 


different  types  of  control  are  required: 


one  which 


stabilizes  the  constant  coefficient  body  modes;  and  a  second 
which  stabilizes  the  periodic  blade  modes.  Control  of  the 
rigid  body  motion  will  be  accomplished  first. 


Controlling  The  Airframe 


Because  the  helicopter's  airframe  (including  the 
horizontal  tail)  is  modeled  as  a  constant  coefficient  system, 
standard  feedback  control  may  be  used  (50:357-384).  Johnson 
(2:794-800)  examines  three  different  feedback  control  systems 
for  a  hovering  helicopter  with  an  articulated  rotor: 
longitudinal  velocity  feedback;  pitch  feedback;  and  lagged 
pitch  feedback.  Further,  he  states. 

In  order  to  achieve  stable  flight,  the  longitudinal 
dynamics  of  the  hovering  helicopter  require  feedback 
control,  either  from  the  pilot  or  from  an  automatic 


control  system  (perhaps  a  mechanical  system,  often  using 
a  gyro) .  The  longitudinal  velocity  and  pitch  attitude 
must  be  sensed  and,  after  appropriate  compensation,  fed 
back  to  the  longitudinal  cyclic  pitch  [2:794]. 

The  pitch  feedback  control  system  Johnson  examines  is 
adopted  for  the  present  analysis  (2:796-798).  The  feedback 
used  is  either  the  helicopter's  pitch  attitude,  pitch  rate. 


or  a  combination  of  both. 


Pitch  attitude  feedback  can 


stabilize,  with  positive  gain,  the  oscillatory  modes  which 
result  from  the  coupling  of  the  pitch  and  the  longitudinal 
velocity.  However,  positive  gain  decreases  the  damping  of 


the  real  root  representing  the  pitch  mode. 


Pitch  rate 


feedback,  also  using  positive  gain,  increases  the  real  root 
damping  and  increases  the  period  and  time  to  double  amplitude 
of  the  oscillatory  roots.  Unfortunately,  this  oscillatory 
response  remains  unstable.  Johnson  suggests  using  a 
combination  of  pitch  attitude  feedback,  which  stabilizes  the 
oscillatory  roots,  and  pitch  rate  feedback  to  keep  the  pitch 
damping  high. 

Consequently,  the  pitch  feedback  to  the  main  rotor’s 
longitudinal  cyclic  pitch  can  be  described  as  follows 
(2:797) : 


■-«[ 


TS  +  18 


(6.1) 


where  the  gain  K  is  positive.  The  lead  T  must  be  large 
enough  so  that  the  zero  is  to  the  right  of  the  open  loop  pole 
(which,  in  hover,  is  principally  due  to  the  main  rotor’s 


pitch  damping).  Otherwise,  the  damping  of  the  real  root  will 


decrease  with  positive  gain  (2:798) . 

Another  type  of  control  system  uses  the  elevator. 

Though  ineffective  at  low  speeds,  deflecting  the  elevator 

through  an  angle  6e  produces  increments  in  both  the  C  and  C. 

in  JL 

of  the  airframe,  which  helps  maintain  helicopter  stability 
throughout  the  flight  envelope  (50:28). 

A  control  system  utilizing  the  elevator  control  surfaces 
as  well  as  pitch  feedback  control  to  the  main  rotor's 
longitudinal  cyclic  pitch  can  be  designed  specifically  for 
the  coupled  rotor/fuselage  system.  A  pitch  feedback  control 
system  will  first  be  designed  for  the  helicopter's  uncoupled 
longitudinal  dynamics  and  compared  with  Johnson's  analysis. 
For  simplicity  this  control  system  will  be  referred  to  as  a 
pitch  attitude  feedback  control  system.  Bemember ,  though, 
that  pitch  rate  feedback  is  also  used  in  conjunction  with 
pitch  attitude  feedback  for  control  of  the  rigid  body  inodes. 

Since  the  pitch  attitude  feedback  control  system  is 
implemented  using  the  main  rotor's  longitudinal  cyclic  pitch, 
a  control  vector  B1  can  be  generated  by  collecting  the  terms 
which  explicitly  use  6g  in  the  airframe’s  uncoupled 
perturbation  equations  [Eqs  (3.50),  (3.51),  and  (3.53)]: 
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Letting  =  K  and  K  ,  =  rK  ,  a  closed- loop  feedback 

control  system  of  the  form 


x2(y/)  =  A(y/)  ♦  B  x2<V') 


(6.3) 


can  be  generated  such  that 
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Setting  T  =  100.0,  the  longitudinal  dynamics  of  the 
uncoupled  airframe  in  hover  can  be  generated  for  various 
gains.  Figure  28  shows  the  root  loci  for  the  pitch  attitude 
feedback  control  to  the  longitudinal  cyclic  for  both  positive 
and  negative  gains,  and  mirrors  Johnson’s  results  (2:797). 
Clearly,  increasing  the  gain  not  only  stabilizes  the 
oscillatory  hover  roots,  but  also  increases  the  real  root’s 


£ 


damping . 

Pitch  attitude  feedback  control  can  be  added  to  the 
uncoupled  body  modes  to  stabilize  the  helicopter  in  hover 
and,  then,  with  the  control  in  place,  the  velocity  can  be 
varied.  Becall  that,  without  any  control.  Figure  14  showed 
that  the  unstable  oscillatory  modes  became  less  unstable  for 
increasing  velocities,  although  they  never  crossed  over  into 
the  stable  region.  In  addition,  the  vertical  and  pitch  roots 
transformed  into  stable  oscillatory  modes.  Figure  29  shows 
the  results  of  letting  K  =  0.60  and  adding  pitch  attitude 
feedback  control  to  the  uncoupled  system  in  hover.  In  hover, 
the  modes  representing  the  coupling  between  the  helicopter’s 
pitch  and  longitudinal  velocity  are  indeed  stabilized.  In 
addition,  the  pitch  root  is  Significantly  more  stable  with 
feedback  control.  However,  as  the  velocity  increases,  the 
oscillatory  modes  become  unstable  and  eventually  become  two 
real  modes ,  one  becoming  more  unstable  and  the  other 
approaching  the  stable  region.  The  vertical  and  pitch 
motions  remain  as  real,  but  stable  roots  with  increases  in 
the  speed . 

It  appears  that  pitch  attitude  feedback  control  is 
effective  in  limited  off-design  cases.  Achieving  stable 
flight  throughout  the  flight  envelope  requires  a  stability 
augmentation  system  which  adjusts  the  pitch  feedback  at 
various  speeds.  At  higher  speeds  the  elevator  can  be 
deflected  to  help  maintain  stable  flight.  In  this  case,  the 
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Figure  29.  Uncoupled  Body  Modes  with  Pitch  Attitude  Feedback 
Control;  No  Reversed  Flow 


control  is  6e  and  the  control  vector,  B1 ,  is  as  follows: 
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Using  combinations  of  pitch  attitude  feedback  through 
the  rotor’s  longitudinal  cyclic  pitch  and  elevator 
adjustments,  it  is  possible  to  generate  stable  roots  for  the 
helicopter’s  uncoupled  longitudinal  dynamics  through  300  m/s. 


For  the  remainder  of  this  analysis,  only  the  pitch  attitude 
feedback  control  system  will  be  used;  the  elevator  will  not 
be  implemented  for  control. 

Stabilization  of  the  coupled  rotor/fuselage  system 
requires  modification  of  the  feedback  control  system.  The 
elevator  control  matrix  remains  the  same,  but  the  matrix 
defining  the  pitch  attitude  feedback  control  changes  because 
now  the  individual  blades  on  the  main  rotor  provide  the 
longitudinal  cyclic  input.  Thus,  for  each  of  the  rotor 
blades,  the  pitch  feedback  vector  B1  is,  without  reversed 
f  low, 
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where,  with  i  denoting  the  blade  and  varying  from  one  through 
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In  the  mixed  flow  region 
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Using  Eq  (6.6)  as  the  control  vector,  integrating  Eq 


(6.3)  over  one  period,  and  summing  over  the  number  of  blades 


on  the  main  rotor  reduces  Eq  (6.6)  to  Eq  (6.2).  Likewise,  a 


root  locus,  obtained  by  adding  pitch  attitude  feedback 


control  [Eqs  (6.4)  and  (6.6)]  to  the  coupled  rotor/fuselage 


system  in  hover  and  varying  K,  airrors  Figure  28.  With  t  = 


100.0  and  K  =  0.60,  the  pitch  feedback  control  system 


stabilizes  the  oscillatory  modes  while  generating  a  much  more 


stable  pitch  root. 


Using  this  control,  the  velocity  is  varied  and  the 


results  summarized  in  Figure  30.  This  figure  shows  that  the 


hover  roots  are  the  same  as  the  hover  roots  of  Figure  29. 


However,  now  the  pitch  attitude  feedback  control  system 


stabilizes  the  body  modes  up  to  122  m/s.  Above  this 
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(a)  Velocity  Root  Locus  (b)  Modes  vs.  Velocity 

Figure  30.  Coupled  Body  Modes  with  Pitch  Attitude  Feedback 
Control;  Reversed  Flow  Added 

velocity,  one  mode  is  slightly  unstable.  The  other  three 
roots  remain  stable  through  300  m/s.  In  addition,  the  pitch 
and  vertical  modes  become  more  stable  as  the  velocity  is 
increased.  Thus,  applying  pitch  feedback  control  to  the 
coupled  system  in  hover  generates  a  wider  range  of  stable 
off -design  airframe  roots  than  the  uncoupled  case.  Again, 
changing  the  pitch  attitude  feedback  control  or  adjusting  the 
elevator's  control  surface  extends  the  airframe’s  stable 
flight  evelope  through  300  m/s. 

Clearly,  pitch  attitude  ftodback  control,  using  a 
combination  of  pitch  attitude  and  pitch  rate  feedback,  is 


beneficial  to  the  airframe  s  coupled  longitudinal  dynamics. 
Figures  31  and  32  show  the  variations  in  the  main  rotor 
blades'  flapping  and  lagging  roots  as  the  velocity  increases. 
In  both  cases  the  system  includes  the  pitch  attitude  feedback 
control  designed  for  hover.  Note  that  Figure  32  shows  the 
velocity  root  loci  of  the  controlled  flap  modes.  The  general 
trends  of  the  controlled  flap  modes  shown  in  Figures  31(a) 
and  32  are  similar  to  the  uncontrolled  roots  summarized  in 
Figures  11(b)  and  17.  At  low  speeds  two  flap  modes  transform 
from  a  pair  of  oscillatory  roots  to  two  real  modes  and  back 
to  a  pair  of  oscillatory  roots.  Two  more  roots  are  real,  but 
stable,  at  speeds  above  79  m/s  while  the  other  four  flap 


modes  generate  stable  oscillations  through  300  m/s, 


Mote 


that,  when  pitch  attitude  feedback  control  is  applied  to  the 
coupled  rotor/fuselage  system,  all  of  the  flap  modes  remain 
stable  through  300  m/s. 

Figure  3l(b)  shows  the  lag  modes  of  the  four  blades  when 
pitch  attitude  feedback  control  is  added  to  the  coupled 
rotor/fuselage  system.  Once  again,  the  magnitudes  of  the 
imaginary  parts  of  the  lag  modes  are  approximately  0.500  from 
hover  through  300  m/s,  so  velocity  root  loci  of  these  modes 
are  not  plotted.  A  comparison  between  Figure  31(b)  and  the 
plot  of  the  real  parts  of  the  uncontrolled  lag  modes  as  a 
function  of  velocity  [Figure  12(b)]  indicates  that,  at  low 
speeds,  the  modes  are  similar.  With  control  added  to  the 
system,  all  eight  lag  modes  are  stable  through  300  m/s. 
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Three  pairs  of  lag  modes  are  nearly  identical  when  compared 
to  the  corresponding  uncontrolled  modes;  the  fourth  pair’s 
roots  are  much  more  stable. 

The  above  discussions  indicate  that  pitch  attitude 
feedback  control  applied  to  the  coupled  rotor  fuselage  system 
in  hover  does  not  significantly  alter  the  individual  blade 
flapping  and  lagging  motions,  except  at  high  speeds  where 
there  are  stabilizing  effects  on  the  two  previously  unstable 
flap  modes.  As  for  the  airframe’s  modes,  the  pitch  attitude 
feedback  control  system  stabilizes  the  oscillatory  roots  not 
only  in  hover  but  also  for  a  fairly  wide  range  of  off-design 
cases  (i.e.,  different  velocities).  Further,  one  of  the 
modes,  after  it  crosses  into  the  unstable  region,  remains 
only  slightly  unstable,  even  at  very  high  velocities. 
Finally,  pitch  feedback  has  a  stabilizing  effect  on  the  pitch 
root  as  the  velocity  increases. 

For  the  coupled  rotor/fuselage  system,  the  net  result  of 
designing  pitch  attitude  feedback  control  through  the  main 
rotor’s  longitudinal  cyclic  pitch  for  hover  is  that  only  one 
body  mode  goes  unstable  as  the  velocity  is  increased  through 
300  m/s.  As  was  mentioned  earlier,  this  mode  can  be 
stabilized  either  by  using  pitch  attitude  feedback  control  or 
by  adjusting  the  elevator.  Consequently,  the  helicopter  can 
achieve  stable  flight  using  standard  feedback  control 
mechanisms  even  at  very  high  speeds. 

Vote  that,  with  this  particular  helicopter 


configuration,  the  blades’  periodic  flapping  and  lagging 
motions  are  stable  throughout  the  flight  envelope  inspected. 
It  is  quite  possible  that  different  helicopter  configurations 
can  generate  unstable  flap  or  lag  modes  with  variations  in 
the  helicopter’s  forward  velocity.  It  is  desirable  to  alter 
the  helicopter’s  parameters  in  order  to  produce  unstable 
responses  in  either  the  flap  or  lag  modes,  and  then  control 
those  modes  using  the  modal  control  technique. 

Therefore,  for  case  ii,  R  =  8.00,  c  =  0.300,  and  = 

mm  (, 

0.10.  This  results  in  unstable  lag  modes  from  hover  through 
16  m/s  for  the  uncoupled  rotor  blades.  Figure  33(a)  shows 
that,  above  this  speed,  all  eight  lag  modes  are  stable.  A 
plot  of  the  real  parts  of  the  uncoupled  flap  modes  as  a 
function  of  helicopter  speed,  shown  in  Figure  34(a), 
indicates  that  each  of  the  four  pairs  of  flap  roots  produce 
stable  oscillations  at  lower  speeds.  As  the  velocity 
increases,  though,  real  roots  are  generated  and  four  of  these 
roots  are  unstable  at  high  speeds. 

Adding  reversed  flow  to  the  uncoupled  blades’  equations 
of  motion  does  not  alter  the  lag  modes  except  at  high  speeds; 
the  modes  are  more  stable  when  reversed  flow  is  included 
[refer  to  Figure  33(a)],  Figure  34(a)  shows  that  adding 
reversed  flow  stabilizes  the  uncoupled  flap  modes  at  higher 
velocities.  Mote  that  these  results  mirror  those  obtained 
using  the  original  rotor  parameters. 

Figure  35  shows  the  uncoupled  body  modes  (without 


300.00 


reversed  flow  being  included)  for  case  ii.  Figure  35(a)  is  a 
velocity  root  locus  of  the  modes  while  Figure  35(b)  plots  the 
real  parts  of  the  body  modes  as  a  function  of  velocity. 
Because  the  airframe’s  equations  are  normalized  with  respect 
to  R  and  O  ,  changing  either  of  these  two  parameters  alters 

ID  ID 

the  uncoupled  body  roots. 

The  coupled  body  modes  (with  reversed  flow  added  to  the 
system)  can  also  be  generated  for  case  ii,  as  shown  in  Figure 
36.  Note  again  that  the  coupling  terms  do  alter  the  modes 
with  increases  in  speed  when  compared  to  the  uncoupled  modes 
(Figure  35).  The  four  blades’  coupled  flap  modes  are  also 
changed  with  case  ii.  Comparing  Figures  11(b)  and  34(b) 
reveals  that  the  magnitudes  of  the  real  parts  of  the  Poincare 
exponents  are  smaller  in  case  ii.  In  addition,  one  pair  of 
flap  modes  is  unstable  at  very  high  speeds.  Figure  37  shows 
the  velocity  root  loci  of  the  four  pairs  of  coupled  flap 
modes  for  case  ii.  The  dashed  and  solid  lines  in  the  figure 
indicates  corresponding  pairs  of  roots. 

Finally,  Figure  33(b)  shows  that,  for  case  ii,  three 
pairs  of  the  coupled  lag  modes  are  unstable  at  velocities  up 
to  approximately  1B6  m/s.  Note  also  that,  at  very  high 
speeds,  a  pair  of  lag  modes  is  unstable.  It  is  interesting 
to  compare  these  results  to  the  uncoupled  lag  modes  generated 
for  case  ii  [Figure  33(a)].  It  appears  that  coupling  the 
airframe  longitudinal  dynamics  to  the  blades’  lag  modes  is 
destabilizing.  It  should  be  noted  that  the  imaginary  part  of 
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Figure  35.  Uncontrolled,  Uncoupled  Body  Modes;  Ho  Reversed 
Flow;  Case  ii 


(a)  Velocity  Root  Locus  (b)  Modes  vs.  Velocity 

Figure  36.  Uncontrolled,  Coupled  Body  Modes;  Reversed  Flow; 
Case  ii 


Figure  37.  Velocity  Hoot  Loci  of  Uncontrolled,  Coupled  Flap 
Modes;  Case  ii 

each  of  the  lag  modes  has  a  magnitude  of  about  0.100  from 
hover  through  300  m/s . 

For  case  ii,  pitch  attitude  feedback  control  can  be 
designed,  in  hover,  to  control  the  airframe's  unstable 
oscillatory  roots.  With  t  =  100.0  and  K  =  0.60,  the  coupled 
rotor/fuselage  sytem’s  roots  are  shown  in  Figures  36  through 
40  as  functions  of  velocity.  Figure  38  indicates  that,  for 
case  ii,  pitch  attitude  feedback  control  stabilizes  the 
coupled  airframe 's  oscillatory  roots  through  300  m/s.  The 
root  representing  the  helicopter’s  pitch  motion  does, 
however,  become  unstable  at  speeds  above  175  m/s.  Figures 


Figure  40.  Velocity  Boot  Loci  of  Coupled  Flap  Modes  with 
Pitch  Attitude  Feedback  Control;  Case  ii 


39(a)  and  40  show  all  eight  flap  modes  to  be  stable  through 
300  m/s  when  pitch  attitude  feedback  control  to  the  main 
rotor's  longitudinal  cyclic  pitch  is  added  to  the  coupled 
rotor/fuselage  system.  Note  again  that  the  dashed  and  solid 
lines  represent  corresponding  pairs  of  flap  roots. 

The  coupled  lag  modes,  shown  in  Figure  39(b),  are 
altered  when  pitch  attitude  feedback  control  is  used. 
Without  any  control  [Figure  33(b)],  six  lag  modes  are 
unstable  up  to  186  m/s,  and  the  other  two  lag  modes  become 
unstable  at  very  high  speeds.  When  control  is  applied  to  the 


coupled  rotor/f uselage  system  for  case  ii,  six  of  the  lag 
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modes  are  unstable  at  speeds  as  high  as  218  m/s.  However, 
the  two  remaining  modes  are  now  stable  throughout  the  flight 
regime  inspected. 

Thus,  for  case  ii,  six  unstable  lag  modes  exist  after 
pitch  attitude  feedback  control  is  added  to  the  coupled 
rotor/fuselage  system.  The  modal  control  technique  can  now 
be  used  to  control  these  unstable  modes. 


Control 1 ing  The 


Blades 


C 


* 


First,  however,  a  control  system  which  uses  the  modal 
control  technique  must  be  developed.  As  was  the  case  in 
controlling  an  isolated  blade’s  flap-lag  motion,  the 
controller  chosen  for  the  coupled  rotor/fuselage  system  uses 
a  conventional  swashplate  mechansm: 


u(y') 


(6. 10) 


The  control  matrix  B(v')  is  generated  by  extracting  the 
collective  and  cyclic  pitch  terms  from  the  A(v')  matrix  found 
in  Eq  (3. 101) : 
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where  4  =  4  ( <-  -  1)  ,  +  steps  from  1  through  H,  and  ^ 
Incorporating  the  B(v')  matrix  terms  derived  for  an 
blade  in  normal  flow  [Eqs  (5.2)  and  (5.3)], 
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Note  that  t_  and  t_  were  defined  in  Eq  (3.71). 
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In  reveraed  flow,  the  Signs  of  the  terua  in  Eq  (6.12) 
are  reveraed.  In  mixed  flow  theae  terms  are  altered  such 
that 


'*♦3.1  *  "e'*')  *  eu 

OO 


A+3 ,2 


4+3,3 


4+4, 1 


4+4,2 


me'£>  *  ei2 


»e<W  *  .13 

u8 


+  ®21 

4o 


mej*')  *  e22 

4c 


4+4,3 


-e'*"  *  *23 

4a 


v>  ■  >  *>  \>\av>  ■_»  •v'-i,  ‘  »  wN  A  w' 


IWWV 


(6. 14) 


Vl.l  =  -r»{[Pll  +  Pll]Cao+[P21+P2l]flao) 

Vl,l  ‘  -r«{K^P12]':“o'[P22*['22)*ao} 

Vl.3  '  -r»{[P13*P’l3]OC,o*(P23*P23)aao} 

V2.I  *  r.{(Pll+Pil)'V(P21*Pil]“o} 

b,,2.2  =  r»{(P12*P'l2)*“o-(P22*P22)=ao} 

Va.3  -  r»{(P13+Pi3]a<:,o-(P234pycao} 
V4.I  =  n»{(P31*P3l]"'z(Pll*Pn)/B»} 
V.,2  =  nm{[P32*P32H(P12*p12]/E»} 
*>,♦4.3  =  nm{(P33*P33]*,z(P13*P13]/B»} 


where 


63o 


(^) 

e3c 


v>) 
e3  8 


m^(w) 


“e,v)  -  Tr»Hvi]*4V,°’v'‘} 

3o  v  v  ' 

"e<v)  *  iT.{(V1)0V2,,v'‘} 

WW 

3s  b 


ma(v) 

®4c 


'  I"'', [l^o (* ‘^o] ■ 1  }**'t *5? [* ^o] cvi 

*kK-|j!)"/r>'('*<!„)  K]] 


(V) 

e4s 


(V)  + 


(* -n)  - [«x2 -^)] 


"ii  - 


P12  2^ 


r2(xflccvl] 


p!3  m  -K{3X'?.*V§K-1)} 


'  *  2 
*21  "  2^ 


(Is  vov,4 [* -3,,o)} 


p22  =  l«[“2K-3)*v{12XX^2(1-,1',o)}<”,'i]*rj2x 

p23  *  ^{K*5K(v,)4*ia(1wH 

P31  '  |^2«*'<-(vi)OVi} 

*32  ‘ 

p33  =  -|,(p‘2X',oOV'i) 


(6. 15) 


11  =  ell+l6^L'^olPlSao  +  C<:,oJ^lCao",,0,oJ|CV>i 

[ca0-*a„)] 

12  =  el2*ia''«[{SXHo,COo]*>J(3'5o*2](CC'o-0ao)}c*'i 


(ac,o*c0,o)  *x  (c“o-s“o)}] 


e  13  "  e 


ca  -sot  lay . 
o  ol  t 


'21  *  e2l‘5-X'»{t'["°'o',CC,o}*3X(OC'o-Sao)}an 


*22  *22 


~T-^r  ft  (cot  -sot  Is  y . 
12  ID  o[  O  O  J  4. 


*23  "  * 23+ 


(aao  +  cao]+9X  [ca0-aao)-4^0  (Gao'8CXo)CV/  J 


The  modal  control  technique  uses  the  control  matrix 
defined  in  Eq  (6.11)  to  control  the  unstable  blade  nodes  in 
the  coupled  rotor/fuselage  system.  Scalar  control  is 
examined  first. 

Determining  the  gain  k(w)  required  to  shift  the  real 
parts  of  the  unstable  blade  modes  of  a  coupled  roter/ fuselage 
system  to  their  desired  locations  necessitates  calculating 
the  modal  controllability  matrix,  g(vO-  Recall  that  Eq  (5.4) 
described  the  g(vO  matrix  used  in  controlling  an  isolated 
blade's  unstable  lag  modes.  If  the  helicopter  is  modeled 


with  a  4-bladed  rotor,  the  control  matrix,  B(y.'),  is  of  order 


20  when  scalar  control  is  to  be  used.  Since  g(y/)  = 

F1(v')B(v/)  [Eq  (4.24)],  g(v')  is  also  a  20x1  vector.  Each 

row  of  g(v0  (where  t  denotes  the  row)  can  be  expressed  as 


g.  .(V')  =  f-/,  B_  .+  f"4  B.  ,  A 
+  *  ^  1B>+ 1  ,  A  +  *  t ,  >+2B/*2  ,  A+  * /,<*  +  4B/»4  ,A 


(0.16) 


where  <•  steps  from  1  through  N,  ^  =  4xt  ,  and  A  is  the 
desired  controller:  collective  pitch;  cyclic  sine  pitch;  or 
cyclic  cosine  pitch. 

The  scalar  collective  pitch  controller  is  first 
examined.  Recall  that,  with  pitch  attitude  feedback  control 
stabilizing  the  body  modes  in  hover,  the  coupled  body  modes 
remain  stable  through  175  m/s  (Figure  38) .  The  addition  of 
collective  pitch  control  to  the  pitch  attitude  feedback 
control  signal  in  hover  shifts  the  most  unstable  pair  of  lag 
modes  to  their  desired  locations.  The  results  are  summarized 
in  Figures  41  through  43. 

Figure  41  plots  the  eight  roots  representing  the  coupled 
lag  modes  when  both  pitch  attitude  feedback  and  scalar 
collective  pitch  controls  are  applied  to  the  coupled 
rotor/fuselage  system.  Figure  41(a)  shows  the  real  parts  of 
the  three  coupled  conjugate  pairs  of  lag  roots,  while  Figure 
41(b)  shows  the  pair  of  real  roots.  This  figure  may  be 
compared  to  Figure  30(b),  which  shows  six  lag  modes  to  be 
unstable  at  speeds  up  to  218  m/s  when  only  pitch  attitude 
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Figure  41.  Coupled  Lag  Modes  as  a  Function  of  Velocity 
Controls;  Case  ii 


feedback  control  is  used  to  stabilize  the  helicopter’s 

coupled  body  modes.  In  hover  and  without  any  modal  control, 

✓ 

the  two  most  unstable  lag  modes  have  Poincare  exponents  at 
5.19259x10  *  ±  1.05750x10  4i.  The  scalar  collective  pitch 
controller  is  designed  to  shift  these  modes  so  that  the  real 
parts  of  the  new  Poincare  exponents  are  at  -0.001.  Figure 
41(a)  shows  this  to  be  the  case.  In  addition,  the  other 
modes  are  unaltered.  Thus,  for  the  coupled  rotor/fuselage 
system,  scalar  collective  pitch  control  functions  as 
predicted  in  Chapter  IV. 

However,  when  scalar  collective  pitch  control  is 


applied  to  an  unstable  pair  of  lag  modes  at  a  single  design 
point,  the  other  three  pairs  of  lag  modes  are  significantly 
altered  in  off -design  cases.  This  is  readily  apparent  when 


comparing  Figure  39(b)  to  Figure  41.  The  latter  figure 
indicates  that  the  lag  modes  for  which  the  control  was 
designed  remains  stable  through  300  m/s.  Another  pair  of 
modes  is  stabilized  at  lower  speeds;  however,  they  are 
unstable  at  higher  speeds.  A  third  pair  of  lag  modes  remains 


unstable 

throughout  most  of 

the 

velocity 

regime 

inspected ; 

they  are 

slightly  stable  at 

low 

speeds 

and 

at 

very  high 

speeds . 

The  final  pair  of 

modes 

,  shown 

in 

Figure  41 (b)  , 

✓ 

splits  off  as  two  real  Poincare  exponents  just  above  the 

✓ 

design  point.  One  Poincare  exponent  becomes  unstable  as  the 
velocity  increases  while  the  other  becomes  more  stable  with 
increasing  velocity.  Except  for  the  two  lag  roots  which 
become  real ,  the  imaginary  parts  of  the  lag  roots  do  not  vary 
much  as  the  velocity  is  increased;  hence,  velocity  root  loci 
are  not  plotted. 

The  collective  pitch  controller,  when  designed  to 
stabilize  a  pair  of  lag  modes  in  hover,  also  alters  the  flap 
modes  in  off -design  cases,  as  shown  in  Figure  42.  Vote  that 
Figure  42(a)  is  a  plot  of  the  real  parts  of  the  flap  modes  as 


a  function 

of 

velocity  while  Figures 

42(b)  and 

42(c) 

plot 

the 

magnitudes 

of 

the  imaginary  parts  of 

the  flap 

modes 

versus 

velocity. 

Comparing  these  figures  to 

Figures 

39(a) 

and 

40 

indicates 

that 

although  the  flap  modes  are  unaltered 

at 

the 

176 


Figure  42.  Coupled  Flap  Modes  as  a  Function  of  Velocity  with 
Controls:  Case  ii 


design  point,  they  do  change  as  the  velocity  increases. 

✓ 

Still,  all  four  pairs  of  Poincare  exponents  remain  stable 
through  300  m/8.  Note  that  one  pair  splits  off  as  two  real 
roots  at  the  upper  end  of  the  velocity  spectrum. 

Finally,  Figure  43  shows  the  body  modes  as  a  function  of 
velocity  when  scalar  collective  pitch  control  is  used  to 
stabilize  two  lag  modes  in  hover.  Figure  43(a)  plots  the 
real  parts  of  the  body  modes  as  the  velocity  is  increased, 
and  Figure  43(b)  plots  the  magnitude  of  the  imaginary  part  of 
the  oscillatory  body  roots  as  a  function  of  velocity.  This 
figure  can  also  be  compared  to  Figure  36,  which  shows  the 
body  modes  as  a  function  of  velocity  when  pitch  feedback  is 
the  only  control  applied  to  the  system.  It  appears  that  the 
body  modes  are  altered  when  scalar  control  is  used  to  control 
one  of  the  lag  modes.  The  oscillatory  roots  shown  in  Figure 
38  become  more  stable  with  increasing  speeds;  at  higher 
velocities  the  roots  become  real,  but  remain  stable.  On  the 
other  hand.  Figure  43  shows  that,  with  Scalar  collective 
pitch  control  added  to  the  coupled  rotor/fuselage  system, 
these  two  roots  become  less  stable  as  the  speed  is  increased. 
The  pitch  root  is  not  too  different  from  the  case  where  only 
pitch  feedback  control  is  used.  The  vertical  mode,  though, 
is  different.  With  just  pitch  attitude  feedback  control  this 
mode  becomes  unstable  at  speeds  above  175  m/s;  when  both 
pitch  attitude  feedback  control  and  collective  pitch  control 
are  applied  to  the  system  in  hover,  the  vertical  mode  becomes 
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Figure  43.  Coupled  Body  Modes  as  a  Function  of  Velocity  with 
Controls;  Case  ii 

more  stable  with  increasing  speed  (see  Figure  43).  Thus, 
modal  control  does  influence  the  airframe’s  modes  in 
off-design  cases. 

To  confirm  the  observations  summarized  above,  various 


design  points  and  off -design  cases  were  examined.  In 
addition,  the  other  four  scalar  controllers,  as  well  as  the 
vector  controller,  were  tested.  In  all  cases,  the  general 
trends  previously  observed  were  repeated.  In  other  words, 
modal  control  does  shift  the  unstable  modes  to  their  desired 
pole  locations  while  leaving  the  other  modes  unaltered. 
Further,  the  controlled  modes  are,  in  general,  more  stable 
than  their  uncontrol led  counterparts  in  off-design  cases. 


However,  the  modal  controllers  do  influence  the  uncontrol led 
inodes  in  off -design  cases.  It  must  be  emphasized,  though, 
that  the  modal  control  technique  does  not  guarantee  a 
system’s  stability  in  off-design  cases.  Enhancing  a  system's 
performance  over  a  wide  range  of  operating  conditions  while 
using  a  scalar  controller  designed  at  a  single  point  is  a 
bonus . 

It  remains  to  be  demonstrated  that  the  modal  control 
technique  can  be  used  to  stablize  more  than  two  modes  at  one 
time.  Since  six  unstable  lag  modes  are  generated  in  case  ii, 
this  question  can  be  addressed. 

Scalar  collective  pitch  control  of  four  unstable  lag 
modes  was  first  attempted  at  various  design  points  using  the 
scalar  control  to  set  the  Sum  of  the  real  parts  of  the  four 
modes.  Unfortunately,  gain  values  could  not  be  found  which 
shifted  all  four  lag  modes  to  stable  locations.  The  same 
results  were  repeated  when  the  other  scalar  and  vector 
controllers  were  used.  Even  though  the  trace  rule  is  still 
valid,  finding  the  required  values  for  the  individual 
components  of  the  gain  matrix  is  quite  difficult  once  there 
are  more  than  two  unstable  modes  to  control. 

On  the  other  hand,  it  is  possible  to  stabilize  more  than 
two  lag  modes  using  the  technique  suggested  in  Chapter  IV. 
That  is,  a  modal  control  system  which  uses  either  a  scalar  or 
vector  controller  can  be  used  to  shift  a  pair  of  unstable  lag 
modes  to  their  desired  locations  at  a  specific  design  point. 


The  resulting  closed  loop  system  is  a  new  linear  system  with 
periodic  coefficients.  Hence,  another  modal  controller  can 


r 

fc 


«? 


\ 


be  designed  for  this  new  system  to  Shift  a  second  pair  of  lag 
modes  to  their  desired  locations.  This  process  stay  be 
repeated  until  all  the  instabilities  are  eliminated  at  the 
design  point. 

The  technique  described  above  is  used  to  control  the  six 
unstable  lag  modes  of  the  coupled  rotor/fuselage  system,  and 
the  results  summarized  in  Table  V  and  VI.  Again,  a  scalar 
collective  pitch  controller  is  used  to  stabilize  all  six 
unstable  roots.  The  design  point  is  chosen  to  be  at  hover, 
and  the  real  parts  of  the  desired  Poincare  exponents  are 
chosen  to  be  at  -0.001.  The  Poincare  exponents  of  the 
uncontrolled  system  (in  hoverlare  listed  in  the  first  column 
of  Table  V  in  the  following  order:  the  first  four  pairs  of 
values  represent  the  coupled  lag  roots  (notice  that  three 
pair  are  unstable) ;  the  next  four  pairs  correspond  to  the 
coupled  flap  roots;  and  the  remaining  four  Poincare  exponents 
(two  real,  one  pair  of  complex  conjugate  roots)  represent  the 
coupled  body  roots. 

The  second  column  of  Table  V  shows  that  a  scalar 
collective  pitch  controller  does  shift  the  first  pair  of 
unstable  lag  modes  to  their  desired  locations,  and  the  other 
modes  remain  unaltered  through  six  decimal  places.  Using 
this  controlled  system,  another  scalar  collective  pitch 
controller  is  designed  to  shift  the  Second  pair  of  unstable 


TABLE  V 


Coupled  Rotor/Fuaelage  Poincare  Exponents 


UNCONTROLLED  SYSTEM 

1  CONTROLLER 

5 . 192593e”5±  1 . O57490e‘  1  i. 

-9.995077e"*±  3.026290e_1v 

3 . O59940e” 5±  1.057554e~‘i 

3 . 059855e" 5±  1 . 057554e“ 1 t 

9 . 0427 15e~  tf±  1 . 057536e' 1 «. 

9.057906e_<S±  1 . O57530e_  1  v 

-3.408090e_3±  1 . 09 1 555e~ 1 i 

-3.408137e-9±  1 . 091553e‘ 1  >. 

-3.795581e_2±  2 . O99087e~ 2 «. 

-3 . 795597e” 2 ±  2.099653e~2C 

-3.805239e“2±  2 . 389089e~ 2 i 

-3 . 805223e~ 2 ±  2.389114e_2t 

-3.806932e_2±  2 . 390504e~  2  i. 

-3.806912e‘2±  2.390474e_S\ 

-4 . 14992 le” 2 ±  2.555194e'2v 

-4 . 14999 le~  2  ±  2 . 555 188e'  2  i 

-1  . 105481e_3±  1 . 248439e~  2 i 

-1 . 104960e“3±  1 . 248497e~  2 i 

-9 . 980934 e  9 ; -6 . 22O205e~  2 

-9.980960e'3 ;-6.220364e“2 

2  CONTROLLERS 

3  CONTROLLERS 

-9.950731e"*±  3 . 025793e" 1 i 

-9.927525e_4±  3.025079e~1i 

-1 -004005e"9±  4.394540e~2i 

-1 .378458e_a±  0.525711e_2i 

9. 234463e_<*±  1 . O57530e~ 1 i 

-1 .000030e~9±  1 . 204 130e~  2 i 

-3 . 438676e~  3  ±  1.091954e_iv 

-3.438994e"S±  1.091939e_1v 

-3.802735e~Z±  2.111320e"Zi 

-3.801327e_2±  2 . 103759e" 2 t 

-3 . 803062e~2±  2.387222e"Zv 

-3 . 804800e_ 2 ±  2 . 389742e' 2  .. 

-3 . 807493e~ 2 ±  2.389073e_Zi 

-3.8O0028e“2±  2 . 389808e' 2 i 

-4 . 108726e~ 2 ±  2.5570O9e‘2v 

-4 . 158455e” 2 ±  2 . 544931e"2  <■ 

-5 . 53 1076e~  4±  1.321801e'Zv 

7.40453Oe_i ; -7 . 432805e" ‘ 

-1 .000727e“2 ;-0. 237134e'2 

-9.979118e'3 ; -0 . 22442 le“ 2 

Poincare  exponents  to  -0.001. 


The  third  column  in  Table  V 


indicates  that  indeed  this  is  the  case.  Unfortunately, 
numerical  inaccuracies  can  be  noticed  in  the  uncontrolled 
modes.  The  roots  representing  the  coupling  between  the 
body’s  pitch  and  longitudinal  velocity  (the  next  to  the  last 
pair  of  modes  shown  in  the  column)  are  different.  In 
addition,  the  other  roots  are  accurate  to  only  four  or  five 
decimal  places.  The  last  column  in  Table  V  shows  the  results 
of  controlling  the  final  pair  of  unstable  Poincare  exponents. 
Again,  the  scalar  collective  pitch  controller  shifts  the 
unstable  roots  to  their  desired  locations.  This  time, 
however,  the  numerical  inaccuracies  are  significant.  The 
roots  representing  the  coupling  between  the  body’s  pitch  and 
longitudinal  velocity  have  now  become  real;  in  fact,  one  is 
unstable.  The  other  modes  are  accurate  to  only  three  or  four 
decimal  places. 

Table  VI  shows  a  comparison  between  the  gains  required 
for  each  of  the  three  scalar  collective  pitch  controllers. 
Notice  that  the  gains  are  all  on  the  same  order  of  magnitude. 

The  results  summarized  in  Table  V  indicate  that  the 
technique  used  to  control  multiple  blade  instabilities  is 
valid.  However,  numerical  inaccuracies  are  a  factor  when 
generating  the  uncontrol led  modes  and  become  noticeable  when 
the  modal  control  technique  is  used  more  than  once  at  a 
particular  design  point. 

As  a  final  note,  the  modal  control  technique  can  also  be 
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TABLE  VI 


Controller  Gains 


CONTROLLER 

2  CONTROLLERS 

3  CONTROLLERS 

-20. 1367 

-20. 1367 

-20. 1367 

9.6640 

9.6640 

9.6640 

0.0000 

25.9670 

25.9870 

0.0000 

21.0596 

21.0596 

0.0000 

0.0000 

22.2350 

0.0000 

0.0000 

5.4816 

used  to  stabilize  the  main  rotor's  blades  with  individual 
blade  controllers.  As  Calico  and  Wiesel  observed,  any  number 
of  unstable  blade  modes  can  be  controlled  with  simple 
one-blade  scalar  controllers  (43:64). 

Therefore,  since  a  combination  of  pitch  attitude 
feedback  to  the  main  rotor’s  longitudinal  cyclic  pitch  and 
elevator  deflections  can  be  used  to  stabilize  the  body 
modes ,  a  coupled  rotor/fuselage  system  with  multiple 
instabilities  can  be  successfully  controlled  throughout  a 
wide  range  of  flight  conditions. 


VI I .  Conclusions 


The  flap-lag  equations  of  motion  of  an  isolated  rotor 
blade’s  flap-lag  motion  and  those  for  a  rigid  helicopter 
containing  four  blades  free  to  flap  and  lag  were  derived.  In 
addition,  control  techniques  were  developed  to  stabilize  both 
systems  for  a  variety  of  flight  conditions.  Some  conclusions 
concerning  the  results  generated  during  this  study  are 
summarized  below. 

Control  of  Blade  Flap-Lag  Motion 

Modal  control  theory,  implemented  through  the  collective 
and  cyclic  pitch  control  mechanisms,  reduced  or  eliminated 
blade  flap-lag  instabilities.  The  pole  placement  technique 
accurately  shifted  unstable  roots  to  their  desired  locations 
while  leaving  the  stable  modes  unaltered. 

The  five  scalar  controllers  considered  all  shifted  the 
unstable  roots  to  desired  locations.  When  used  at  off-design 
points  one  controller,  using  collective  pitch,  failed  to 
stabilize  the  blade  over  a  wide  range  of  p  and  parameters 
up  to  fj  =  0.06.  The  other  scalar  controllers  reduced  the 
blade’s  unstable  regions  from  hover  through  /u  =  0.50. 

Similar  results  for  the  vector  controller  were  obtained. 
However,  this  controller  eliminated  the  unstable  regions  in 
the  p  versus  parameter  space  for  both  hover  and  ij  -  0.40. 


! 

i; 


In  fact,  the  vector  controller  produced  results  which  were 


better  than  those  generated  by  any  of  the  scalar  controllers. 


In  summary,  both  the  scalar  and  vector  controllers 


shifted  the  unstable  roots  to  desired  locations,  and,  using  a 


single  design  point,  reduced  or  eliminated  the  unstable 


regions  of  a  blade’s  flap-lag  motion  throughout  a  wide 


variety  of  flight  conditions. 


Control  of  Coupled  Botor/Fuselage  Motion 


Coupling  the  flap-lag  equations  of  motion  of  four  rotor 


blades  to  a  rigid  airframe  altered  the  flap,  lag,  and 


airframe  roots.  A  pitch  attitude  feedback  controller,  which 


used  a  combination  of  pitch  attitude  and  pitch  rate  feedback 


to  the  rotor’s  longitudinal  cyclic  pitch,  was  designed  to 


stabilize  the  uncoupled  body  modes  in  hover.  The  same 


feedback  control  mechanism  also  stabilized  the  body  modes  of 


the  coupled  rotor/fuselage  system,  not  only  at  the  design 


point  but  also  over  a  wide  range  of  off-design  velocities. 


Although  not  specifically  considered  in  this  research 


effort,  stabilizing  the  rigid  body  modes  throughout  the 


flight  envelope  can  be  accomplished  by  scheduling  the  gains 


in  the  pitch  attitude  controller. 


The  coupled  rotor/fuselage  system  was  also  used  to 


further  analyze  the  modal  control  technique  developed  in 


Chapter  IV.  Even  though  the  system  was  of  order  20,  there 


were  no  numerical  difficulties  in  designing  a  modal 


controller  to  place  a  single  pair  of  unstable  roots.  It  was 

shown  that,  at  a  specific  design  point,  scalar  and  vector 

/ 

control  shifted  a  pair  of  unstable  Poincare  exponents  to 
desired  locations  without  altering  the  other  modes  in  the 


system.  In  addition,  the  controllers,  designed  at  a  single 
point,  stabilized  the  controlled  lag  modes  over  a  wide  range 


of  off-design  cases.  However,  the  controllers  altered  the 
other  modes  in  off-design  cases. 

Finding  gains  which  stabilized  more  than  two  roots 
proved  to  be  very  difficult.  However,  the  modal  control 
technique  was  used  to  eliminate  multiple  blade  instabilities 
by  first  controlling  a  pair  of  unstable  roots  at  a  specific 
design  point.  The  resulting  closed  loop  system  was  a  new 
linear  system  with  periodic  coefficients.  Another  modal 
controller  was  designed  for  this  new  system  to  shift  a  second 
pair  of  unstable  roots  to  desired  locations.  This  process 
was  repeated  until  all  instabilities  were  eliminated. 
Numerical  inaccuracies,  unfortunately,  became  noticeable  when 
modal  control  was  used  more  than  once. 
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and  vector  control  are  successfully  used  to  stabilize  the  blade's  motion. 

Coupling  the  flap-lag  equations  of  motion  of  four  rotor  blades  to  a  rigid 
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